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CHAPTER  1 
INTRODUCTION 


The  need  for  multidimensional  signal  processing  in  manifold  areas  of 
applications  is  now  well  recognized.  The  important  role  of  passive  digital 
filtering  and  related  multidimensional  modelling  schemes  in  the  general  area  of 
multidimensional  signal  processing,  and  in  particular,  in  the  areas  of  image 
processing,  target  tracking  etc.,  has  already  been  documented  in  an  earlier 
report  [1]  and  thus,  its  discussion  will  not  be  undertaken  here.  A  brief 
description  of  the  results  obtained  during  the  course  of  the  present 
investigation  subsequent  to  those  reported  in  [1]  follows. 

In  chapter  1  properties  of  various  multidimensional  polynomials  arising  in 
studies  of  passive  (lossless)  discrete  multidimensional  systems  are 
investigated.  Reactance  Schur  polynomials  and  immittance  Schur  polynomials 
occurring  respectively  as  the  denominators  (and  numerators)  of  discrete 
reactance  functions  and  discrete  positive  functions  are  introduced  and  their 
properties  studied.  Role  of  these  polynomials  in  scattering  or  immittance 
descriptions  of  passive  discrete  time  domain  multiports  are  brought  out.  The 
interrelation  between  classes  of  multidimensional  polynomials  arising  in 
discrete  systems  and  the  corresponding  classes  of  polynomials  in  the  context  of 
continuous  systems  is  also  studied  via  the  mechanics  of  bilinear 
transformation. 

In  chapter  2  the  problem  of  structurally  passive  synthesis  of  multidimensional 
digital  filters  of  the  quarter-plane  causal  type  as  an  interconnection  of  more 
elementary  building  blocks  directly  in  the  discrete  domain  has  been  addressed 
via  the  factorization  of  the  chain  matrix,  the  hybrid  matrix  and  the  transfer 
function  matrix  associated  with  a  prescribed  multidimensional  lossless 
two-port.  By  exploiting  recent  results  on  the  discrete  domain  representation  of 
such  matrices  a  generalized  lossless  two-port  matrix  has  been  introduced  to 
present  all  three  factorizations  in  an  unified  setting.  Necessary  and 
sufficient  conditions  for  factorability  as  well  as  algorithm  for  computing 
these  factors  when  they  exist  are  obtained.  In  particular,  it  is  shown  that  in 
one-dimension  the  factorizations  can  always  be  performed.  Thus,  in  1-D, 
discrete  domain  algorithms  for  synthesizing  previously  unpublished  internally 


passive  structures  as  well  as  alternative  methods  of  synthesis  for  more 
conventional  structures  such  as  the  cascade  structure  are  also  obtained  as  a 
byproduct  of  our  discussion.  Since  most  multidimensional  applications  dictate 
that  the  filter  be  either  symmetric  or  (quasi)  antimetric,  special  attention  is 
paid  to  the  problem  of  synthesis  of  these  subclasses  of  multidimensional 
lossless  two-ports. 

Due  to  difficulties  inherent  to  the  mathematics  of  quarter  plane 
multidimensional  filtering,  (e.g.,  polynomial  nonfactorability  etc.)  the 
alternative  approach  of  studying  passive  filters  of  the  fully  recursive 
half-plane  causal  type  is  undertaken  in  chapter  3  for  the  first  time  in  the 
literature.  Apart  from  the  tractability  of  analylis  and  design,  such  filters 
have  the  potential  to  maximally  exploit  the  currently  emerging  parallel  (VLSI 
and/or  optical)  architechture,  when  implementation  is  called  for.  Thus,  passive 
and  lossless  two-dimensional  digital  one-ports  as  well  as  two-ports  of  the 
fully  recursive  half-plane  type  are  introduced  and  are  characterized  in  terms 
of  their  transfer  function  descriptions.  An  algorithm  for  the  structurally 
passive  synthesis  of  filters  having  such  recursive  structure  is  then  derived 
from  these  representation  results  as  an  extension  of  a  recent  Schur  type 
algorithm  for  the  synthesis  of  discrete  lossless  two-ports.  Design  methods  for 
important  practical  cases  when  frequency  response  of  the  filter  is  required  to 
have  specific  symmetries  are  also  presented. 

Finally,  results  are  sumnarized,  conclusions  are  drawn  and  reconmendations  for 
further  work  are  made  in  chapter  4. 

Each  of  the  following  chapters  are  self  contained  and  can  be  read 
independently.  For  similar  discussions  in  the  open  literature  we  refer  to  the 
publications  (2],  (3]  and  ( 4 ]  in  the  following. 


2 


References 


[1]  S.  Basu,  Multidimensional  Filtering  Investigations,  Technical  report, 

RADC-TR-86-231 .  A178  745. 

[2]  S.  Basu,  A.  Fettweis,  New  results  on  stable  multidimensional 
polynomials,  Part  II:  Discrete  case,  IEEE  trans.  on  CAS,  November,  1987. 

[3]  S.  Basu  and  A.  Tam,  On  the  synthes izability  of  multidimensional  lossless 
two-ports,  Midwest  Symp.  on  Circuits  and  Systems,  Syracuse,  1987.  Also, 
submitted  for  journal  publication. 

[4]  S.  Basu,  Synthesis  and  design  of  structurally  passive  fully  recursive 
half-plane  2-D  digital  filters,  Allerton  Conference,  University  of 
Illinois,  September  1987. 


3 


2.1.  INTRODUCTION: 


CHAPTER  2 

SOME  NEW  RESULTS  CN  STABLE 
HULTIDIKENS IGNAL  DIGITAL  FILTERS 


The  scattering  Hurwitz  polynomials  have  been  introduced  recently  as  the 
denominators  of  (rational)  bounded  functions  arising  in  studies  on  passive 
multidimensional  (k-D)  networks  (1].  Subsequently,  the  denominator  polynomials 
of  (rational)  reactance  functions  and  (rational)  positive  functions  have  also 
been  characterized  as  the  reactance  Hurwitz  and  the  immittance  Hurwitz 
polynomials  respectively  (2].  We  also  refer  to  { 3 ]  for  a  related  discussion.  In 
view  of  recent  interest  in  synthesis  of  passive  digital  filter  networks 
directly  in  the  discrete  domain,  a  study  of  properties  of  the  corresponding 
polynomials  arising  in  discrete  time  systems  seems  highly  relevant. 
Furthermore,  it  is  now  known  to  researchers  in  the  field  that  unlike  the  1-D 
case,  in  multidimensions  discrete  counterparts  of  certain  continuous  domain 
results  need  not  always  be  true  (4-6).  A  separate  treatment  for  discrete  time 
systems  is  thus  needed  for  a  more  careful  analysis.  The  first  step  in  this 
direction  has  been  taken  in  (7)  by  studying  the  properties  of  discrete 
scattering  Hurwitz  polynomials.  In  the  present  paper  further  properties  of  this 
latter  class  of  polynomilas  are  studied,  discrete  reactance  Hurwitz  polynomials 
and  discrete  immittance  Hurwitz  polynomials,  ocurring  as  the  denominators  (and 
consequently  the  numerators)  of  discrete  reactance  functions  and  discrete 
positive  functions,  are  introduced  arid  their  properties,  as  they  relate  to  both 
single-port  and  multi-port  passive  digital  networks,  are  studied  in  detail.  It 
is  possible,  at  least  in  principle,  to  undertake  the  abovementioned  discussion 
via  the  utilization  of  multiple  bilinear  transform  and  analogous  results 
already  existing  [2]  for  continuous  systems.  However,  it  turns  out,  due  to 
difficulties  of  the  type  elaborated  in  (4-6),  that  such  an  indirect  approach  is 
neither  mathematically  elegant  nor  is  it  desirable  from  the  standpoint  of  a 
self-contained  theory  of  passive  multidimensional  discrete  systems.  Therefore, 
the  present  paper  has  been  organized  in  such  a  way  that  it  can  be  read 
indepedently  of  all  existing  publications  on  continuous  systems. 

In  what  follows  we  will  consider  polynomials  a-a(z)  in  k-variables  z  - 
(z^,*2, .. •Z]') •  All  polynomials  and  rational  functions  are  assumed  to  be 
functions  of  k-variables  unless  otherwise  specified.  A  k-variable  polynomial  a 
will  said  to  involve  a  variable  z^  if  the  indeterminate  z^  actually  exists  in 
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at  least  one  of  the  monomials  comprising  the  polynomial  a.  A  polynomial  will  be 
called  non-trivial  if  it  involves  at  least  one  of  the  variables  z.,  i  »  1  to  k. 
More  generally,  the  above  comments  apply  if  a  is  a  rational  matrix  function  in 
z  -(z^,Z2» . • .z^) .  If  a  -  a(z)  is  written  as  a  polnomial  in  z^  as: 


ni  v 
a  -  E  a  z. 

v-0 


(2.1.1) 


where  the  coefficients  av's  are  polynomials  in  the  remaining  variables  with  an 
?  0,  then  n^  is  called  the  partial  degree  of  a  in  the  variable  z.  and  is  to  M 
denoted  by  deg^a.  Occassionally,  we  will  write  in  a  compact  notation  zr- 
(Zj^Zj.  • .  .z^) .  Two  polynomials  will  be  said  to  be  relatively  prime  if  they  do 
not  have  a  proper  (i.e.,  a  non-cons tant )  factor  in  common.  A  rational  function 
will  be  said  to  be  in  irreducible  form  if  the  numerator  and  denominator 
polynomials  are  relatively  prime.  The  discrete  paraconjugate  a  of  a  polynomial 
a  -  afz^Zj/ . .  .z^)  in  k-variables  z  -  (z^,z-. . .  ,z^)  is  defined  as: 


),  and  the  notation  z  denotes  the 


**  n.  *  ^  1  „ — 1  1  n 

a  -  a.z  ,  jjherj*  a  -  a  (Zj^  ,z2  , ...,zk  ),  and  the  notation  z-3  denotes  the 

monomial  Zj^  ^  ni  "  deg^a.  The  superscript  *  denotes  complex 

conjugation.  Note  that  the  discrete  paraconjugate  of  a  is  not  necessarily  equal 

to  a.  The  polynomial  a  will  be  called  discrete  selfparaconjuqate  if  a  -  ya, 

where  y  is  a  (necessarily  unimodular)  constant.  In  addition,  a  will  be  called 

discrete  paraeven  or  discrete  paraodd  according  as  r*l  or  y»-l.  The  notation 

|z(  <  1^  will  mean  |z^|  <1  for  all  i  ■  1  to  k.  Obvious  variations  of  this 

notation  with  the  symbol  >  replaced  by  -,  <,  <,  >  etc.  will  be  used. 

The  tr*:ms  sequentially  almost  complete  set  0  or  sequentially  infinite  set  0,  as 
defined  in  (2]  will  be  used.  In  the  present  context,  however,  the  elements  of 
the  m- tuple  (0^,©2, . . ,0ffl)  belonging  to  ©will,  unless  otherwise  specified,  be 
chosen  from  the  field  of  real  numbers  modulo  2n. 

We  will  have  the  occasion  to  consider  the  polynomial  a  associated  with  the 

polynomial  a,  as  defined  in  the  following.  Let  a-a(z)-a(j+a1z+...anz  ,  where  the 

coefficients  a^'s  are  polynomials  in  several  variables  over  the  field  of 

complex  numbers.  Then  associated  with  a  we  define  a  polynomial  a  via  the 

formal  algebraic  operation  as:  a_»a  (z)*bn+b.z+. . .b  z  ,  where  b.«(n-2i)a.  for 

z  z  •”  u  i  n  x  i 

i  -  0,1,... n.  Further  elementary  properties  of  the  polynomial  a  are  derived 

z 
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in  appendix  2 «c. More  important  definitions  will  be  introduced  as  they  become 
necesary  in  the  main  text. 

Various  classes  of  polynomials  devoid  of  zeros  in  the  unit  polydisc  will,  in 
general,  be  referred  to  as  the  multidimensional  Schur  polynomials.  This  is  a 
departure  from  the  previously  adopted  terminology  in  (7),  18],  but  is  consistent 
with  the  terminology  used  for  the  corresponding  class  of  one-variable 
polynomials  in  the  literature.  Properties  of  the  widest  sense  Schur 
polynomials,  selfparacon jugate  Schur  polynomials  and  scattering  Schur 
polynomials  are  discussed  in  section  2.2.  Section  2.3  contains  discussion  on 
elementary  properties  of  multidimensional  discrete  positive  functions.  In 
section  2.4  reactance  Schur  and  immittance  Schur  polynomials  respectively 
occuring  as  the  numerators  (and  consequently  denominators)  of  discrete 
reactance  functions  and  discrete  positive  functions  are  introduced.  Extensions 
of  some  of  these  results  to  passive  multiports  are  discussed  in  section  2.5. 
Finally,  interrelationship  between  the  various  classes  of  multidimensional 
Schur  polynomilas  and  the  corresponding  classes  of  multidimensional  Hurwitz 
polynomials  is  examined  via  the  artifice  of  bilinear  transformation  in  section 
2.6,  and  conclusions  are  drawn  in  section  2.7.  An  appendix  is  included,  where 
some  general  properties  of  multidimensional  polynomials  proved  elsewhere  are 
stated,  and  a  few  basic  formulae  occurring  in  studies  of  discrete  systems  are 
derived.  For  consistency  in  the  logical  sequence  in  which  the  proofs  of  results 
are  arranged  in  the  present  paper,  theorems  2. 2. 2. 7  and  2. 2. 2. 8  should  be  read 
after  2.3.6  and  before  corollary  2. 2. 3. 7.  However,  theorems  2. 2. 2. 7  and  2. 2. 2. 8 
are  incorporated  earlier  in  section  2.2.2  rather  than  in  section  2.2.3,  beacuse 
this  allows  for  a  more  systematic  categorization  of  the  properties  of  various 
classes  of  multidimensional  Schur  polynomials. 
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2.2.  VARIOUS  CLASSES  OF  SCHUR  POLYNOMIALS 
2.2.1.  Widest  sense  Schur  polynomials: 

Definition  2. 2. 1.1:  A  polynomial  a  is  called  widest  sense  Schur  if  a/0  for  |z| 

<  A- 

Theorem  2. 2 .1.2:  If  a(z)  is  a  widest  sense  Schur  polynomial  in  k-vfcriables 

z-fz^^, . .  .z^),  then  the  polynomial  a'(z')  in  (k-1)  variables  z'  obtained  by 
freezing  z^  in  |z.  |£l  is  either  widest  sense  Schur  or  is  identically  equal  to 
zero.  Furthermore,  the  latter  instance  can  occur  only  for  finitely  many  values 
of  z ^  on  \z^\  -1. 

Proof:  Obviously  a'(z')  cannot  be  zero  in  \z'\<l  for  any  |z-|<l,  because  if 

for  some  z^  -  z1Q  in  |z^q|<1,  a'  (z' )  has  a  zero  z'  -  z*q  with  |z^|<l,  then 
a(Zi0,z^)  -  0,  which  contradicts  the  widest  sense  Schur  property  of  a(z). 
Furthermore,  if  a'(z')  ■  0  for  some  z'  »  z£  in  |z£|<l,  and  for  z^  -  z1Q  on 

IZj^l-1  then  the  following  two  possibilities  arise.  Firstly,  if  a'(z')  *  0 
i.e.,  a(z10,z')  *  0  for  any  z'  then  (z-z1Q)  is  a  factor  of  a(z) .  Obviously, 

then  there  are  finitely  many  values  of  z^q  on  |z^|-l,  such  that  (zi-zig)  is  a 
factor  of  a(  z) .  This  proves  the  second  part  of  assertion  of  the  theorem. 
Next,  to  prove  the  first  part,  assume  that  a'(z')  ji  0,  i.e.,  (z^-z^q)  is  not  a 
factor  of  aiz).  Then  by  moving  the  variable  z^  by  sufficiently  small  amount 
inside  |z^|<l  from  zj“ziq»  and  by  invoking  the  continuity  property  of  zeros  of 
a  polynomial  we  would  be  able  to  construct  a  zero  of  a(z)  with  |  z^  |  <1  and 
1*0 1 <i*  v'hich  contradicts  the  widest  sense  Schur  property  of  a(z)  and  is  thus 
impossible. 

A  repeated  application  of  the  above  theorem,  along  with  the  definition  of 
sequentially  almost  complete  [cf.  Appendix  A]  set  yield  the  following  result. 

Theorem  2 . 2 . 1 . 3 :  let  a(z)  be  a  widest  sense  Schur  polynomial  in  k-variazles  z. 
Then  there  exists  a  sequentially  almost  complete  set  of  m-tuples  ©m,  of  order 
m<k,  such  that  the  (k-m)  variable  polynomial  obtained  by  freezing  m  of  the  k 
variables  z^  at  z^  -  z^q  with  z^  ■  expfj©^,  0  <  9^  <  2n  for,  say  i  -  1  to  m. 


7 


is  widest  sense  Schur  if  (©^©^ . .  .©m)  £  ©m  and  is  identically  equal  to  zero  if 

(e1,e2,...em)  %  ©m- 

2.2.2  Selfparacon jugate  Schur  Polynomials: 

Definition  2. 2. 2.1:  A  polynomial  a  is  called  self-paraconjugate  Schur  if  a  t  0 
for  |z|<l^  and  3  -  ya,  where  y  is  a  complex  number. 

~  k  ni 

Since  3  -  a  n  zi  we  have  for  z^-exp(j©^),  i-1  to  k,  |3|  -  ja|  »  |a  |  -  |a|. 
Thus  the  coAsAant  y  in  definition  2. 2. 2.1  is  necessarily  unimodular  i.e., 
lTl-1. 

Lemma  2. 2. 2. 2:  Let  a  be  a  selfparaconjugate  Schur  polynomial  in  k-variables  z, 
and  be  the  (k-1)  variable  polynomial  obtained  by  freezing  in  a  one  of  the 

variables,  say  z^,  on  |z^|  *  1.  Then  deg.a  -  deg^a^  for  the  remaining 

variables  i  -  1  to  k  but  i/X- 


Proof:  Assume  X/l .  Let  us  write  the  polynomial  a  as  in  (2R1.1^  with  i-1.  Then 

it  follows  from  property  2.B8  in  the  appendix  that  yan  -aQ  II  z,  l,  Furthermore, 

since  deg.an<n. ,  it  follows  from  the  last  equation1 that*" ya  -  3n(  n  z.1), 

where  p^^  are  some  nonnegative  integers.  If  an  *  0  for  some  on  |z^|-1 

irrespective  of  the  other  variables  then  (z-z.i)  must  be  a  factor  of  a  ,  and 

au  n1 

thus  a  factor  of  Sg.  This  later  conclusion,  however  imply  that  (z^-z^gj  is  a 
factor  of  a^,  and  therefore  a-0  for  z^-z^g  and  z^-0  irrespective  of  other 
variables,  and  consequently,  either  a^«0  or  is  not  widest  sense  Schur. 
However,  since  a^O  by  hypothesis,  this  is  in  contradiction  with  theorem  2.1.1. 
Therefore,  an  jtO  for  any  z^  on  |z^|-l.  The  validity  of  the  present  lemma  foe 
other  values  o£  i  can  similarly  be  demonstated  by  writting  a  as  a  polynomial  in 
z^  with  coefficients  as  polynomials  in  remaining  variables. 


Theorem  2. 2. 2. 3:  If  a  is  a  self-paraconjugate  Schur  polynomial  in  k-variables 
then  the  polynomial  a'  in  (k-1)  variables  obtained  by  freezing  any  one  of  the 
variables,  say  z^,  is  also  self-paraconjugate  Schur  for  an  almost  conplete  set 
of  values  on  |z^|  -  1 


Proof:  From  theorem  2.1.1,  a'  is  widest  sense  Schur  for  an  almost  conplete  set 
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of  values  of  9^^  in  0£©1<2r/  where  Zj*exp(j©^)  i.e.,  for  an  almost  complete  set 

z ^  on  IZjJ  *  1.  It  remains  to  show  that  a'  -  y' .a',  where  y’is  a  constant.  By 

substituting  z1-exp(j91)  in  S-ya  it  follows  that 

*  *1  k  n. 

a  (exp(je1),z2  ,...zk  )(exp(  j©^))  n  zi  -  ya(  (exp(  j©1 )  ,z1# . .  .zk) ,  where 

n^-deg^a.  lhe  desired  result  then  follows  by  noting  that  the  right  hand  side  of 
the  last  equality  is  ya',  whereas  the  left  hand  side  due  to  the  fact  that 
deg^a-deg^a'-n^  (lemma  2.2.2.2)is  equal  to  (exp(  jn^)  )S' . 

Repeated  use  of  the  above  result  yields  the  following  theorem. 

Theorem  2 . 2 . 2 . 4 :  Let  a  be  a  self-paraconjugate  Schur  polynomial  and  a  be  the 
polynomial  obtained  by  freezing  m,  m<k,  of  the  k-variables,  say 
z^,  i  ■  1  to  »,  on  Iz^  -  1  i.e.,  z^  -  exp(j©^)  ;  0£9^<2ji.  Then  there  exists  a 
sequentially  almost  complete  set  ©m  of  m-tuples  of  order  m  such  that  for  any 
(©1,©2, . ..©m)€0m,  the  polynomial  is  selfparacon jugate  Schur. 

Theorem  2. 2. 2. 5;  Let  a  be  a  self-paraconjugate  Schur  polynomial.  There  exists 
a  sequentially  almost  complete  set  ©  of  k-tuples  of  order  (k-1)  such  that  a  -  0 
for  any  z ^  -  exp(j©i);  0<©^<2fi,  i  -  1  to  k,  and  (©^©^ . .  .©^c©. 

Proof;  From  theorem  2. 2. 2. 4,  the  polynomial  obtained  by  freezing  (k-1)  of 
the  k  variables,  say  z^,  i-1  to  (k-1)  on  z^  -  exp(je^)  ;  (k©^<2n  is 
selfparacon  jugate  Schur  for  any  (©^,©2, . .  .©^le©^,  where  is  a 

sequentially  almost  complete  set  of  order  (k-1).  Since  the  polynomial  - 
ak-l*zk*  is  self-Paraconju9ate  Schur,  for  |zk|<l  as  well  as  for  | | >1, 

and  thus  has  zeros  on  |zk|-l  only.  Therefore,  there  is  a  sequentially 
almost  complete  set  ©  of  k-tuples  of  order  (k-1)  such  that  a  -  0  for 

( ©^ » ©j  1  •  •  •  ©j^ )  c®» 

we  refer  ahead  to  definition  2. 2. 3.1  for  the  statement,  but  not  the  proof  of 
the  following  theorem. 

Theorem  2. 2. 2. 6;  A  widest  sense  Schur  polynomial  a  can  be  expressed  as  a 
product  of  a  selfparacon jugate  Schur  factor  and  a  scattering  Schur  factor. 

Proof;  Let  a-de,  S-df,  where  d  is  the  gcd  between  a  and  a,  and  thus  e  and  f 
are  prime  polynimials.  From  property  2.B6  in  the  appendix  then  d  -  ^d.  i.e., 
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d  is  discrete  self-paracon jugate.  Since  a  is  widest  sense  Schur  so  are  both  d 
and  e.  Thus  d  is  selfparacon jugate  Schur.  We  further  claim  that  e  is 
scattering  Schur.  Clearly,  d.f  -a  -  d.e  -  y.d.e  and  hence  f  -  ye,  where  the 
second  equality  follows  from  property  2.B1  in  the  appendix.  Therefore,  if  e 
and  e  had  a  common  factor  then  e  and  f  would  not  be  relatively  prime. 

The  following  theorem  characterizes  self par aeon jugate  Schur  polynomials. 

Theorem  2. 2. 2. 7;  A  polynomial  a  is  self-paraconjugate  Schur  if  and  only  if  a/0 
for  |z|<l.  as  well  as  for  |z|>l. 

Proof:  Necessity:  aj^O  in  |z|<l  by  definition.  Furthermore,  since  a  is 
discrete  self-paraconjugate,  a  zero  of  a  in  (z|>l  would  imply  existence  of  a 
zero  of  a  in  |z|<l,  and  is  hence  excluded. 

Sufficiency:  If  a  /  0  for  |z|<^L,  then  by  theorem  2. 2. 2. 6  a  is  product 
of  a  self-paraconjugate  Schur  factor  and  a  scattering  Schur  factor.  However,  a 
factor  of  latter  type,  due  to  theorem  2. 2. 3. 6  must  have  a  zero  in  |z|>l,  thus 
implying  a  zero  of  a  in  |z|>l,  which  is  excluded. 

Theorem  2. 2. 2. 8:  The  irreducible  factors  of  a  selfparaconjugate  polynomial  a 
are  also  selfparaconjugate  Schur. 

Proof:  Obviously,  the  irreducible  factors  of  a  are  widest  sense  Schur.  Thus, 
due  to  theorem  2. 2. 2. 6,  these  factors  are  either  selfparaconjugate  Schur  or 
scattering  Schur.  However,  presence  of  a  factor  of  the  latter  type,  due  to 
theorem  2. 2. 3. 6,  would  indicate  that  a  has  a  zero  in  |z|>l,  which  is  ruled  out 
in  theorem  2. 2. 2. 7. 

Lemma  2. 2. 2. 9:  Let  a  be  a  selfparaconjugate  Schur  polynomial  involving  the 
variable  z^.  Then  deg^  a  -  deg^a2  for  X  »  1  to  k. 

Proof:  If  X-i  the  proof  is  obvious  from  the  definition  of  az  .  Next,  let  X*i, 
and  a  be  written  as  in  (2.1.1).  Since  a  is  selfparaconjugate  Schur,  and  in 
particular  widest  sence  Schur,  cannot  have  a  factor  z^  for  any  X,  beacuse 
otherwise  a  would  be  zero  for  z.«0,  z.-O,  and  arbitrary  values  of  the  remaining 
variables.  Therefore,  Sq  must  contain  a  monomial  not  involving  z^  .  Since  from 
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property  2.B8  in  the  appendix,  an  -ra0  n  ,  the  ^ast  stated  property  of  aQ 

yields  that  deg.a„  -nr  Furthermore*  sin£l*in  a_  -  I1  b  z?,  b  -  -n.  .a  we 
-  a  n.  a  z.  vi  n.  1  ni 

have  degxbz  -  deg^an  -  n^.  Consequently,  deg^az  -n^. 

2.2.3  Scattering  Schur  Polynomials: 

Definition  2.2. 3.1:  A  polynomial  a  is  called  scattering  Schur  if  a  f  0  for 
|z|<l,(i.e.,  a  is  widest  sense  Schur)  and  if  a  and  a  do  not  have  any  common 
factor. 

The  term  discrete  scattering  Hurwitz  has  earlier  been  used  for  the  above  class 
of  polynomials  in  (7], (8). 

Theorem  2. 2 .3. 2:  If  a  is  a  scattering  Schur  polynomial,  then  a  cannot  have  a 
discrete  selfparaconjugate  factor. 

A 

Proof:  If  a  had  a  factor  d  i.e.,  a  -  de,  with  d  -  od,  where  a  is  a  constant, 

H.  AAA  A 

then  a  -  de  -  doe,  and  thus  a  and  a  would  not  be  relatively  prime. 

Corollary  2. 2. 3. 3:  A  scattering  Schur  polynomial  a-a(z)  in  one  variable  is  a 
strict  sense  Schur  polynomial  i.e.,  a/0  for  |z|<l. 

Theorem  2. 2. 3. 4:  Let  a  be  a  scattering  Schur  polynomial  in  k-variables  then  the 
polynomial  a'  in  (k-1)  variables  obtained  by  freezing  any  one  of  the  variables, 
say  Zy  with  z^  -  exp(j0^);  (k0^<2«,  is  also  scattering  Schur  for  an  almost 
complete  set  of  values  of  0^. 

Proof:  Let  a'(z' )  -  a(z1Q,z')  with  z1Q  -  exp(j01Q),  where  01Q  is  a  fixed  value 
of  0^  in  O<0^<2ji.  Due  to  theorem  2. 2. 1.2,  a'  is  either  identically  equal  to 
zero  or  is  a  widest  sense  Schur  polynomial.  However,  in  the  former  case 
(z-Zj^q),  which  is  selfparaconjugate  Schur,  is  a  factor  of  a,  and  thus  excluded 
due  to  theorem  2.2. 3.2. 

A 

Next,  the  polynomials  a  and  a  are  relatively  prime  by  definition  2. 2. 3.1, 
and  thus  the  (k-1)  variable  polynomials  a'  -  a(z' )  -  (a)z  and  [aJz 
with  z^Q-exp(  j0^Q) ,  due  to  property  2.A4  in  the  appendix,  Jre^elatively  ^rii§ 
for  an  almost  complete  set  ©  of  values  of  0^g.  However,  since  from  property 
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«  n.o  k  p. 

2.B7,  [aj_  _  -  z.ia'  H  z,  ,  p4>0,  we  have  that  a'  and  a'  are  relatively 

2l“z10  10  i-2  1  1 

prime  polynomials,  and  thus  a'  is  scattering  Schur  for  all  e^c©. 

Extending  the  above  result  via  a  repeated  application  we  obtain  the  following: 
Theorem  2.2. 3.5:  Let  a  be  a  scattering  Schur  polynomial  and  let  am  be  the 
polynomial  obtained  from  a  by  freezing  m,  m<k  of  the  k-variables,  say  z^,  i  -  1 
to  m  on  z^  -  exp(je^),  (ke^<2n.  Then  an  is  also  scattering  Schur  for  any 
choice  of  the  m-tuple  (©^©j, . .  .©m)e©m,  w*iere  *s  a  sequentially  almost 
complete  set  of  m-tuples  of  order  m. 

Theorem  2. 2. 3. 6:  A  scattering  Schur  polynomial  a  must  have  zeros  in  |z|>l. 

Proof:  Let  a,  ^  be  the  polynomial  obtained  from  a  by  freezing  (k-1)  of  the 

k-variables,  say  z^,  for  i  -  1  to  (k-1)  on  ■  exp(j9i>  0<9i<2ii.  Then  due  to 
theorem  2. 2. 3. 5,  there  exists  a  sequentially  almost  complete  set  of  order 
(k-1)  such  that  for  any  * '^-l^^-l  ^  PolYno,,lial  aj(_i  is  scattering 

Schur,  and  is  thus,  in  view  of  corollary  2. 2. 3. 3,  a  strict  sense  Schur 
polynomial  in  the  variable  zfc  only.  Therefore,  ak  l  has  zeros  in  |zkl>l. 
Consequently,  a  has  zeros  for  |zi|  «  1,  i  «  1  to  (k-1)  and  | | >1 .  Next,  by 
continuously  moving  the  variables  z^^  in  the  regions  |zi|>l  for  i  -  1  to  (k-1) 
by  sufficiently  small  amounts,  and  invoking  the  continuity  property  of  zeros  of 
a  polynomial  as  a  function  of  its  coefficients,  it  follows  that  a  has  zeros  in 
1*1 >1. 

Theorem  2. 2. 3. 6  can,  in  fact,  be  strengthened  in  the  following  form. 

Corollary  2. 2.3. 7:  An  irreducible  polynomial  a  is  scattering  Schur  if  and  only 
if  a  #  0  for  |z|<l  and  a  has  at  least  one  zero  in  |z|>l. 

Proof:  Necessity  of  the  theorem  is  obvious  in  view  of  theorem  2. 2. 3. 6.  To 

prove  sufficiency,  let  us  note  that  a  is  widest  sense  Schur,  and  therefore,  by 

virtue  of  its  irreducibility  and  by  theorem  2. 2. 2. 6  is  either  a 

\ 

selfparaconjugate  Schur  or  a  scattering  Schur  polynomial.  However,  due  to 
theorem  2. 2. 2. 7,  a  self-paraconjugate  Schur  polynomial  cannot  have  zeros  in 
|z|>l.  Thus,  a  is  a  scattering  Schur  polynomial. 
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Theorem  2 . 2 . 3 . 8 ;  Let  a  be  a  widest  sense  Schur  polynomial.  Then  a  is  also 
scattering  Schur  if  and  only  if  the  zeros  of  a  on  \z\-l  does  not  form  a 
sequentially  almost  complete  set  of  order  (k-1). 

Proof:  Sufficiency:  Assume  for  contradiction  that  d  is  the  nonconstant 
greatest  common  factor  between  a  and  a.  By  property  2.B6  in  the  appendix  d  is 
a  discrete  self-paracon jugate  polynomial.  Invoking  theorem  2. 2. 2. 5,  it  follows 
that  a  -  0  for  zi»exp(j0i),  0<9^<2ii,  for  i  -  1  to  k  with  (©^,©2, ..  .©^e©,  where 
6  is  a  sequentially  almost  complete  set  of  order  (k-1).  However,  this  latter 
conclusion  contradicts  the  fact  that  the  zeros  of  a  on  |z|-l  cannot  form  a 
sequentially  almost  complete  set  of  order  (k-1).  The  polynomial  d  is  therefore 
a  constant ,  thus  proving  the  scattering  Schur  property  of  a. 

Necessity:  If  a-0  for  z ^  -  exp(j©^),  0<©i<2n,  i  -  1  to  k  with  (e^,©^, . .  .©^e©, 
where  ©  is  a  sequentially  infinite  set  of  ordgr  (k-1),  then  a(exp(j©1), 
exp(j©2),...  exp(j^))  -  a  (exp(j©1),...exp( j©k)  n  exp( jni©i)-0,  where  n.  - 
deg^a,  thus  implying  that  a  and  a  would  have  Sequentially  infinitely  many 
common  zeros  of  order  (k-1)  on  (zj  -  1.  Property  2. A3  in  the  appendix  then 

A  “ “ 

implies  that  a  and  a  would  have  a  common  factor,  which  is  impossible  if  a  is 
scattering  Schur. 

Theorem  2.2. 3.9:  (a)  Factors  of  scattering  Schur  polynomials  are  scattering 
Schur.  (b)  Conversely,  products  of  scattering  Schur  polynomials  are  also 
scattering  Schur. 

Proof:  (a)  Let  a-bc  be  a  scattering  Schur  polynomial.  Obviously,  then  b  and  c 
are  widest  sense  Schur,  beacuse  a  is  so.  Furthermore,  since  due  to  property 

A  A  A 

2.B1  in  the  appenix,  a-bc,  there  cannot  be  a  nontrivial  common  factor  between  b 
and  b  or  between  c  and  c,  beacuse  otherwiswe  a  and  a  would  not  be  relatively 
prime,  and  the  scattering  Schur  property  of  a  would  thus  be  violated. 
Therefore,  b  and  c  are  both  scattering  Schur  polynomias. 

(b)  Conversely,  If  a-bc  with  b  and  c  scattering  Schur,  then  clearly  a  is  widest 
sense  Schur.  Due  to  theorem  2. 2. 2. 6  a  is  the  product  of  a  scattering  Schur 
factor  and  a  selfparaconjugate  Schur  factor.  However,  irreducible  factors  of 
the  latter  type  are  also  sel  f  par  aeon  jugete  Schur  due  to  theorem  2. 2. 2. 8,  and 
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would  thus  be  contained  in  either  b  or  c,  which  in  view  of  theorem  2. 2. 3. 2, 
violates  the  scattering  Schur  property  of  b  or  c.  Therefore,  a  is  a  Scattering 
Schur  polynomial. 

Theorem  2.2.3.10:  Let  a  be  a  scattering  Schur  polynomial  in  k  variables,  then 
the  (k-1)  variable  polynomial  a'  obtained  by  freezing  any  one  of  the  variables, 
say  z^,  in  |zil<l  is  also  scattering  Schur. 

Proof:  Assume  i  «  1,  and  z^  to  be  frozen  at  z^  »  z^.  Obviously,  a'  -  a'(z') 
-  a(z^Q,z')  is  widest  sense  Schur.  We  only  need  to  show  that  a'  and  a'  are 
relatively  prime  polynomials.  If  a'  and  a'  are  not  relatively  prime  then  the 
greatest  common  factor  d  between  them  is  by  property  2.B6  discrete 
self-paracon jugate,  and  thus  is  a  sel f pa racon jugate  Schur  polynomial.  Invoking 
theorem  2. 2. 2. 5  it  then  follows  that  d*»0,  and  thus  a  -  0,  for  z^  -  exp(j©^);  0 
<  0^<2it  for  i  -  1  to  k,  with  0'  *  (O-,©^, . . .0^)e©' ,  where  0'  is  a  sequentially 
almost  complete  set  of  order  ( k— 2 ) .  Since  a  »  0  for  z ^  •  exp(j@i)  O<0i<2n,  i  - 
2  to  k  with  0'e  ©  and  for  z^  ■  z^q,  theorem  2.2.1. 3  along  with  the  fact  |z1Q|<l 
yields  that  a  -  0  for  z^  «  exptj©^,  O<0i<2n,  with  0'e  0'  and  arbitrary  z^. 
Therefore,  by  restricting  z^  on  z ^  «  exp(j01)  it  is  possible  to  assert  that  a  - 
0  for  zi  -  exp(j0i),  0<®£<2r,  for  all  <M0lf02, . .  .0^)  e  0,  where  ©  is  a 
sequentially  almost  complete  set,  and  thus  a  sequentially  infinite  set  of 
k-tuples  of  order  (k-1).  This  latter  conclusion,  in  view  of  theorem  2 .2.3.8 
violates  the  scattering  Schur  property  of  a.  The  polynomials  a  and  a'  are  thus 
relatively  prime. 

Lemma  2.2.3.11:  If  a  is  a  widest  sense  Schur  polynomial  and  b  a  further 
polynomial  such  that  |b/a|<l  for  z  on  |z|»l  whenever  a/0,  then  |b/a|<l  for 
|z|<l.  Furthermore,  if  |b/a|»l  for  some  z  in  |z|-l  then  it  is  impossible  to 
have  |b/a|«l  for  some  z  in  |z|<l  unless  b/a  is  independent  of  z^  for  all  i 
i.e.,  a  constant. 

Proof:  Assume  for  the  purpose  of  induction  that  the  theorem  is  true  for 
polynomial  and  rational  functions  in  (k-1)  variables.  Then  due  to  theorem 
2.2.2. 3,  the  polynomial  a'  in  (k-1)  variables  obtained  by  freezing  z^  in  a  at 
z^Q  on  |z|»l  is  widest  sense  Schur  for  almost  all  (i.e.,  except  finitely  many) 
choices  of  z^q  on  |z|«l.  Also,  let  us  define  the  (k-1)  variable  polynomial 
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d'-d'(z')-d(z10#z).  Now  since  a'-a'(z'M)  implies  a(z10,z')j*0,  if  | b/a | <1  for 
jz|-l,  whenever  a^O  then  we  have  that  |b'/a'l<l  for  |z'  |«1,  whenever  a'7<0. 
Consequently,  by  invoking  induction  hypothesis  we  may  assert  that  |b'/a'|<l  for 

|z'l<l. 

Consider  next  the  polynomial  and  b^-b(zlfz^),  where  z^  is 
considered  frozen  in  |z'  |<1.  Obviously,  a^O  in  |z^|<l.  Furthermore,  from  the 
conclusion  of  the  last  paragraph  it  follows  that  ja^/b^l^l  for  almost  all  z ^  on 
|z^|-l.  Thus,  via  an  application  of  the  known  maximum  modulus  principle  for 
rational  functions  of  one  variable,  we  conclude  that  |bj/a^|<l  for  | | <1^. 
Since  z£  is  chosen  arbitrarily  in  |z'|<l,  the  latter  conclusion  yields  |b/a|<l 
for  |z|<_l. 

Finally,  the  last  statement  of  the  present  theorem  follows  from  the  the  known 
result  [91  that  a  function  of  several  complex  variables  cannot  reach  a  maximum 
at  a  point  interior  to  the  domain  of  holomorphy  unless  it  is  a  constant. 

Theorem  2.2.3.12:  (a)  If  a  is  a  nontrivial  scattering  Schur  polynomial  then 
"  1  "  ‘ 

|a/a|<l  for  |z|<l.  (b)  Conversely,  if  b  and  a  are  relatively  prime  polynomials 
such  that  | b/a | <1  for  |z|<l  then  a  is  either  a  scattering  Schur  polynomial  or 
a  nonzero  constant. 

Proof;  (a)  Since  |a/a|-l  for  |z|»l,  whenever  a^O,  the  result  follows 
ininediately  from  lemma  2.2.3.11. 

(b)  Clearly,  a^O  for  |z|<l,  beacuse  if  a(Zg)-0  with  |Zq|<1  then  in  order  for 
|b/a|<l  to  be  satisfied  we  would  necessarily  have  b(z^)-l  i.e.,  the  rational 
function  b/a  would  have  a  nonessential  singularity  of  second  kind  at  z^. 
Consequently,  b/a  would  also  have  a  singularity  of  first  kind  in  an  arbitrarily 
small  neighbourhood  of  z^  lying  entirely  within  |_z | <_1  (10],  thus  violating  the 
condition  |b/a|<l  for  |z|<l. 

Next,  we  claim  that  if  a  is  not  a  constant,  and  if  a(Zg)«0  for  some  z^  with 
1^0 1"—'  t^ien  b(*g)“0.  To  substantiate  this  claim  consider  the  open  connected 
set  2  of  points  z  lying  in  a  neighbourhood  of  z^  as  well  as  inside  |z|<l. 
Obviously,  |b/a|<l  and  a^O  for  zc2.  It  then  follows  by  invoking  the  continuity 


of  the  function  b/a  in  2  that  if  a(Zg)-0  then  also  have  b(Zg)-0.  Now,  since  the 
polynomials  a  and  b  are  relatively  prime  the  set  of  zeros  conmon  to  them 
cannot,  due  to  propertty  A3,  form  a  sequentially  almost  complete  set  of  order 
( k— 1 ) .  Thus,  the  set  of  zeros  of  a  on  |zj-l  does  not  form  a  sequentially  almost 
complete  set  of  order  (k-1).  It  then  follows  by  invoking  theorem  2. 2. 3. 8  that  a 
is  a  scattering  Schur  polynomial . 

Corollary  2.2.3.13;  Let  the  rational  function  p-b/a  be  such  that  a  is  widest 
sense  Schur  and  |p|<l  for  all  those  z  on  |z|-l  for  which  a^O.  Then,  if  p-b^/ag 
in  irreducible  form  then  Sq  is  scattering  Schur. 

Proof:  Follows  immediately  from  lemma  2.2.3.11  and  theorem  2.2.3.12. 


2.3.  ELQttNTARY  PROPERTIES  OP  (k-D)  RATIONAL  DISCRETE  POSITIVE  FUNCTIONS: 


Definition  2.3.1:  A  rational  function  C  will  be  called  discrete  positive  if  Re 
C  >  0  for  |^|  <1,. 

A  discrete  positive  function  C  -  jC,  where  C  is  a  real  constant,  is  said  to  be 
trivial.  All  other  discrete  positive  functions  are  said  to  be  non-trivial. 

Definition  2.3.2:  A  discrete  positive  function  C-b/a  with  a,  b  polynomials, 
will  be  said  to  be  a  discrete  reactance  function,  if  in  addition,  it  satisfies: 
(b/a)  +(b/a)  -  0. 

L*»mnw  2.3.3:  let  C  -  C'(z')  be  the  well  defined  rational  function  in  ( k— 1 ) 

variables  obtained  by  freezing  in  the  k-variable  rational  function  C  *  C(z)  the 
variable  z ^  on  |z^|-l.  Then  if  C  is  discrete  positive  real  so  is  C' . 

Proof;  From  definition  2.3.1  via  the  use  of  arguments  similar  to  those  used  in 
the  proof  of  lerana  2.5.11. 

Theorem  2.3.4:  A  discrete  positive  function  C  is  non-trivial  if  and  only  if 

ReC  >  0  for  |z|<l. 

Proof;  If  C  is  a  trivial  positive  function,  then  ReC  -  0  for  all  z,  and  in 
particular  for  |z|<l.  On  the  otherhand,  if  C  involves  one  of  the  variables, 
say  Zy  then  let  us  freeze  all  other  variable  z'  -  (z2,z3, . .  .z^)  at  z'  -  z^  in 
\z' |<1  in  C  and  let  the  resulting  function  be  called  C' .  Then  C' ,  due  to  lemma 
2.3.3,  is  a  non-trivial  discrete  positive  function  of  one  variable  z ^  only,  for 
which  it  is  well  known  that  ReC'>0  for  | z^ 1 <1 .  Since  z£  is  chosen  arbitrarily 
in  |z|<l,  the  present  theorem  is  established. 

Theorem  2,3.5:  If  C«b/a,  with  a  and  b  polynomials,  is  a  discret  positive 
function  in  irreducible  form  then  both  a  and  b  are  widest  sense  Schur 
polynomials.  Furthermore,  neither  a  nor  b  may  have  selfparaconjugate  Schur 
factors  of  multiplicity  larger  than  one. 

Proof:  Consider  the  nonconstant  rational  function  p  -  (1-C)/(1+C)  - 
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(a-b)/(a+b),  which  is  necessarily  in  irreducible  form.  We  have  ReC  >  0,  and 
consequently  |p|<l  for  |z|<l.  Therefore,  due  to  theorem  2.2.3.12(b),  (a+b)  is 
scattering  Schur.  lhus,  p  is  holomorphic  in  |z|<l^»  and  consequently,  the 
maximum  modulus  theorem  for  functions  of  several  complex  variables  [9]  implies 
that  | p | <1  in  |z|<l.  Now,  if  b(Zg)*0  but  a(Zg)-0  or  b(Zg)-0  but  a(Zg)#0  with 
IZglcl  then  we  would  obviously  have  |p(Zg)|»l,  which  is  impossible. 
Furthermore,  if  for  some  Zg  with  |Zgj<l,  a(Zg)  -  b(zQ)-0,  then  p  would  have  a 
nonessential  singularity  of  the  second  kind  at  z  -  zQ.  However,  p  would  then 
also  have  a  singularity  of  first  kind  in  an  arbitrarily  small  neighbourhood  of 
z»Zg  contained  in  |z|<l  (10],  which  is  ruled  out  by  |p|<1. 

To  prove  the  second  part,  let  b  -  e.fn  where  e,f  are  polynomials  such  that 
f  is  discrete  selfparacon jugate  and  does  not  divide  e.  Then  due  to  theorem 
2. 2. 1.3  there  exists  a  sequentially  almost  complete  set  ©'  of  ( k— 1 )  tuples  of 
order  ( k— 1 )  such  that  the  one  variable  polynomial  b'-b'fz^)  in  the  variable  z^ 
obtained  by  freezing  in  f  the  remaining  variables  at  z^  -  exptj^),  i  -  2  to  k 
with  e/-(e2,e3, . . . 

almost  complete  set  of  order  (k-1).  Let  a'  and  f'  be  the  polynomials  obtained 
by  freezing  respectively  in  a  and  f  the  variables  z^  at  z^^  -  exp(j0i)  for  i  -  2 
to  k.  Then,  since  the  pair  (a,f)  is  coprime,  due  to  property  2.A4  in  the 
appendix,  there  exists  a  sequentially  almost  complete  set  0"  of  (k-1)  tuples  of 
order  (k-1)  such  that  for  any  choice  of  0'-(  02,03,  •  •  *\) c0"  •  P3*1  (a',fr) 

is  also  coprime.  Therefore,  if  0'e(0'n  0")  then  the  one  variable  rational 
function  C'(z^)  »  b'  (z^)/a'  (z^)  is  well  defined,  and  in  its  irreducible  form 

the  numerator  polynomial  b'  contains  the  factor  fr-f'(z^)  of  multiplicity  n. 
Note  that  the  polynomial  V ,  due  to  theorem  2.2.2. 3,  is  selfparaconjugate 
Schur.  Furthermore,  due  to  lemma  2.3.4,  C'  is  a  discrete  positive  function  in 
one  variable,  and  is  thus  devoid  of  multiple  zeros  on  the  unit  circle  |z2|-l. 
Consequently  n-1. 

Assertion  2.3.6:  Let  C  •  b/a  be  a  discrete  positive  function.  Then  C  is  a 

reactance  function  if  and  only  if  deg.b  »  deg. a  for  all  i  -  1  to  k  and  b/a  -  - 

A  4*  1  1 

b/a. 

Proof:  let  t  -  d/c  in  irreducible  form,  where  b  -  ed  and  a  -  ec,  e  being  the 

g.c.d  between  a  and  b.  Then  due  to  property  2.B1  in  the  appendix,  b/a  -  b/a 


©l^e©'  is  also  widest  sense  Schur,  where  ©'  is  a  sequentially 


~  -  -  -  *  11  ■ -llj  • 

implies  d/c  -  -  d/c  -  -(d/c)  D  z  ,  where  nci  -  deg.c,^ndi  “  deg^  for  i  - 
1  to  k.  Also,  due  to  property^. B2  in  the  appendix  fi,  and  d  does  not  contain  a 
factor  2^  for  any  i.  Furthermore,  c  and  d,  being  numerator  and  denominator  of 
irreducible  discrete  positive  function,  cannot  have  a  zero  in  |z|<l  due  to 
theorem  2.3.5,  and  thus  cannot  have  a  factor  z^  for  any  i.  Consequently,  nci  - 
ndi  for  all  i,  therefore,  degib  «  deg.e  +  n^  -  deg^  +  nci  -  deg^a  for  all  i  - 
1  to  k.  Furthermore,  b/a  *  d/c  -  -  d/c  -  -  (d.e)/(c.e)  -  -  b/a,  where  the  last 
equality  follows  from  property  2.B1  in  the  appendix.  Conversely,  if  deg.b  - 
deg^a  for  i  -  1  to  k,  then  obviously  b/a  -  b/a.  Thus,  if  C  -  b/a  is  a  discrete 
positive  function  with  b/a  -  -  b/a,  we  would  also  have  b/a  +  b/a  -  0  implying 
that  t  is  a  reactance  function. 

Notice,  however  that  a  discrete  positive  function  can  have  a  non-zero 
difference  in  partial  degree  between  its  numerator  and  denominator.  Consider, 
t  -  l+((l-z1z2)/(l+z1z2)  ]  -  ll+z^]-1 ,  which  is  a  discrete  positive  function. 

Theorem  2.3.7:  If  a  rational  function  t  -  b/a  satisfies:  (i)  a  is  scattering 
Schur  (ii)  Re  C>0  for  |z|*l»  whenever  a*0  then  ReC  >  0  for  |z|  <  1  i.e.,  C  is  a 
discrete  positive  function. 

Proof :  Let  us  assume  for  the  purpose  of  induction  that  the  theorem  be  true  for 
(k-1)  variables.  Also,  let  the  polynomial  a'-a' [z' )-a(ZjQ,z' )  and  the  rational 
function  C'-C' (z' )-C(z1Q,z' )  be  defined  by  freezing  respectively  in  a  and  in  i 
the  variable  z^  at  z^q  on  |z^|-l.  Since  a  is  scattering  Schur,  due  to  theorem 
2. 2. 3. 4,  a'  is  scattering  Schur  for  an  almost  complete  set  of  values  z^q  on 
| z^ |*1.  For  those  z'  on  |£' l"i  for  which  a'*0,  we  have  a(z1Q,z')*0,  and 
consequently,  due  to  condition  (b)  of  the  present  theorem,  ReC'*ReC(z10,z' )>0. 
Then  due  to  te  induction  hypothesis  it  follows  that  Ret'>0  for  )z'l<l  for  an 
almost  complete  choice  of  z^Q  on  |z^|*l. 

Consider  next  a^-a^z^-aU^,^)  and  Cj-Cj^z^-CfZj/Z^) ,  where  z^  is  frozen  in 
Iz^^l.  Then  from  the  conclusions  of  the  last  paragraph  it  follows  that  the 
one-variable  rational  function  satisfies  ReC^>0  for  an  almost  complete  set 
of  values  on  |z^|*l.  Furthermore,  a^,  due  to  theorem  2.2.3.10  and  corollary 
2. 2. 3. 3,  is  a  strict  sense  Schur  polynomial  in  z^  only.  Therefore,  from  a 
classical  one-variable  result  it  follows  that  ReC^>0  for  |z^|<l.  Since  z£  is 
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arbitrarily  chosen  in  |z'  |<1,  the  last  condition  yields  that  ReC>0  for  |z|<l, 
thus  completing  the  proof  of  the  present  theorem. 


Remark  2.3.8;  The  above  theorem,  in  the  one  variable  case,  can  be  obtained  from 
an  application  of  lemma  2.2.3.11  with  the  nonrational  function  exp(-t). 
However,  since  a  multidimensional  version  of  lemma  2.2.3.11  valid  for 
nonrational  functions  is  not  known  to  us,  a  separate  proof  of  theorem  2.3.8  is 
needed. 


Theorem  2.3.9;  If  C»b/(ac)  is  a  discrete  positive  function  in  k-variables  z, 
where  a  is  a  scattering  Schur  polynomial  and  c  is  a  reactance  Schur  polynomial 
in  one  variable  only  then  C  can  be  decomposed  as  t-(c'/c)+(a'/a) *  where  c'/c  is 
a  one  variable  discrete  reactance  function  and  a '/a  is  a  discrete  positive 
function  in  k-variables  z. 


Proof:  Let 
distinct.  Cons 


n 

c-  £  c. , 


iair1 


where 


vw 


with 


the  rational  function  K  (z')  -  [ (z.-a  )b/acl 

V  —  XV  z 


1  and  o^'s  are  necessarily 
,  which  is 


v  -  *  v  ‘i  “v 

obviously  holomorphic  in  |z'|<l.  Furthermore,  for  all  z£  with  tylo*  is 
the  residue  of  the  one-variable  discrete  positive  function  Cj-C^z^-CIz^Zq) 
at  z, ■<*..,  and  is  thus  real  and  positive.  Since  a  function  of  several  conplex 


'l'V 
variables, 


constant 


which  is 
ye  have 


h«h(z)«b-a  E  Kic/ci 
each  v  . , 
straightforward 
c'-  E  l^c/c^. 
thefceiore,  I 


real  and  holomorphic  in  a  domain  must  be  equal  to  a 
that  K^( z' )-K^-constant.  Also,  since  the  polynomial 
satisfies  h(ov,z')»0,  (z^-a^)  must  be  a  factor  of  h  for 

h«ca'  for  some  polynomial  a'«a'U).  It  then  follows  by 
algebraic  manipulation  that  t«(c'/c)+(a'/a),  where 
Clearly,  c'/c  is  a  one-variable  discrete  reactance  function,  and 
Re(c'/c)-0  for  |z^|-l,  whenever  c*0.  Since  C  is  a  discrete  positive 
function,  ReC>0  for  |z|-l,  whenever  C  is  regular.  Therefore,  Re(a'/a)>0  for 
|z|-l,  whenever  a*0.  It  then  follows  via  the  use  of  theorem  2.3.8  that  a'/a  is 
a  discrete  positive  function. 


Theorem  2.3.10:  Let  t-b/a,  where  a  and  b  are  relatively  prime  polynomials.  Then 
C  is  a  discrete  positive  function  if  and  only  if  (a+b)  is  scattering  Schur,  and 
ReC>0  for  |z|-l,  wherever  C  is  regular. 
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Proof:  Consider  the  rational  function  p-(l-C)/(l+CMa-b)/(a+b)  in  irreducible 
form  and  notice  that  ReC>0  if  and  only  if  | p | <1 . 

Necessity:  Since  ReC>0  for  jz|<l  we  have  | p | <1  for  |z|<l.  Therefore,  due  to 
theorem  2.2.3.12(b),  (a+b)  is  scattering  Schur.  Finally,  the  fact  that  Re£>0 

for  all  regular  points  of  C  on  |z|«l  follows  from  lemma  2.3.3. 

Sufficiency:  If  a+b*0  but  a*0  for  some  z  then  |p|-l.  On  the  other  hand,  if 
a+b*0  but  a*0  for  some  z-z^  with  |Zq|-1^  then  p  is  regular  at  z^,  and 
consequently,  ReC(ZQ)>0,  thus  implying  |p(Zg)|<l.  Therefore,  | p ) <1  for  |z|-l, 
whenever  a+b*0  i.e.,  p  is  holomophic.  The  result  then  follows  from  lemma 

2.2.3.11  via  the  use  of  the  fact  that  (a+b)  is  scattering  Schur. 

We  have  the  following  important  result. 

Theorem  2.3.11:  If  a  is  a  widest  sense  Schur  polynomial  involving  one  of  the  k 

variables,  say  z.,  then  a  /a  is  a  discrete  positive  function. 

1  z. 

Proof:  Assume  i  -  1.  Since  a  is  widest  sense  Schur  the  polynomial  a  in  one 

variable  z1  obtained  by  freezing  in  a  the  variables  z^i-2  to  k  inside  Iz^  1, 

is  widest  sense  Schur.  The  result  then  follows  from  property  2.C3  in  the 

appendix  by  noting  that  Re(a  /a)  -  Re(a  /a)>0  for  | z | <1 . 

Z1  Z1  ~ 

Theorem  2.3.12:  If  a  is  a  selfparaconjugate  Schur  polynomial  involving  one  of 

the  k  variables,  say  z. ,  then  a  /a  is  a  discrete  reactance  function  (not 

x  z. 

necessarily  in  irreducible  form) . 


Proof:  Assume  i  »  1,  and  a  -  ra.  The  fact  that  a  /a  is  a  discrete  positive 

-  z 

function  has  been  proved  in  theorem  2.3.11.  If  a  is  expressed  as  a  polynomial 


in  z^  with  coefficients  as  polynomials  in  z' 

1,  then  az  -  E^z?  ,  where  b^  -  (n1-2v)av 
Z1  ”  x  v 


i.e.,  in  the  form  (2.1.1)  with  i  - 
for  \^0,1, . . .n^.  Furthermore,  from 


lenma  2. 2. 2. 9  we  have  that  deg. a  -  deg. a  -  n.  for  i  ■  1  to  k.  Therefore, 

1  Zm  1  1 

.  n.-  _  k  n.  n.  k  n. 

a„  -(  Exb  z7v)  11  Z;  -  (  E  -  (n.-  2v)a„  ,  z  /(  Hz.)  (via  the  substitution  b 

Z1  v-0  v  1  i-1  1  v-0  1  nl~v  i-2  1  „ 

«  ni  n. 

■  (n.-2v)a  )  -  E1  -(n1-2\>)ya  zv  (via  the  use  of  equation  (2.a1)*-y  xbz  « 

1  v  \*«0  1  v-0  v 
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-  ya  .  Therefore,  (a  /a)  +  (a  )/a  -  0.  The  fact  that  a  /a  is  a  discrete 
21  Z1  Z1  Z1 

reactance  function  then  following  by  invoking  assertion  2.3.6. 

Corollary  2.3.13:  If  a  is  an  irreducible  widest  sense  (discrete 

selfparaconjugate)  Schur  polynomial  involving  the  variable  z.  then  a  /a  is  a 

12. 

discrete  positive  (reactance)  function  in  irreducible  form. 

Proof:  Immediately  follows  from  theorems  2.3.11,  3.12  and  property  2.C2  in  the 
appendix. 

Theorem  2.3.14:  If  a  is  a  widest  sense  Schur  polynomial  involving  the  variable 
z^  then  the  polynomial  az  is  also  widest  sense  Schur. 

Proof:  let  a_  /a  -  c/b,  a  ■  bd  and  a„  -  cd,  d  being  the  g.c.d  between  a„  and 

a.  Then  due  to  theorem  2.3.11,  c/b  it  discrete  positive  and  consequently*  due 

to  theorem  2.3.5,  b  and  c  are  widest  sense  Schur  polynomials.  Furthermore,  the 

polynomial  d,  being  a  factor  of  a,  which  is  widest  sense  Schur,  is  also  widest 

sense  Schur.  Thus  a  -  c.d  is  widest  sense  Schur. 

zi 

Theorem  2.3.15:  If  C  -  b/a  is  a  nontrivial  discrete  positive  function  in 

irreducible  form  then  f  *  xa  +  jyb  is  widest  sense  Schur  for  all  real  x  and  y 
not  simultaneously  zero.  Conversely,  if  the  polynomial  f  -  xa+jyb  is  widest 
sense  Schur  for  all  real  x  and  y  not  simultaneously  zero  then  either  b/a  or 
(-b/a)  is  a  discrete  positive  function. 

Proof:  If  b/a  is  a  nontrivial  discrete  positive  function  in  irreducible  form 

then  due  to  theorem  2.3.5,  a*0,  b*0  for  |zj<l.  Consequently,  if  f(Zg)  ■  0  for 
some  Zq  in  |£|<1.  then  b/a  -  j(x/y)-,  implying  Re(b/a)  -  0  for  z  -  z^  in  |z|<l. 
The  latter  conclusion,  in  view  of  theorem  2.3.4  contradicts  with  the  fact  that 

(b/a)  is  a  nontrivial  discrete  positive  function.  Conversely,  if  f  is  widest 

sense  Schur  for  all  real  x  and  y  except  x  -  y  -  0,  then  a  and  b  are  also  widest 
sense  Schur,  because  if  a(z^)  -  0  (or  b(Zg)-O)  with  |  z^  | <1,  then  f(z^)  •  0  for 
y  ■  0(or  x  -  0)  and  for  arbitrary  x(or  y).  Thus,  Re(b/a)  is  a  continuous 
function  of  z  in  |z|<l.  Furthermore,  Re(b/a)  *  0  for  any  z  in  |z|<l,  because 
otherwise  b/a  -  -  jy/x. 


i.e.,  xa+jby  ■  0  for  some  real  x  and  y  not  simultaneously  zero.  Thus,  by 
continuity  of  b/a  in  |z|<l,  if  Re(b/a)>0  (or  <0)  for  some  z  in  |z|<l,  then 
Re(b/a)>0  (or  <  0)  for  all  |z|<l. 

Corollary  2.3.16;  The  rational  function  C  -  b/a,  where  b  and  a  are  relatively 

prime  polynomials,  is  a  nontrivial  discrete  positive  function  if  and  only  if 
2  2  2  2 

(x  a  +y  b  )  is  a  widest  sense  Schur  polynomial  for  all  real  x  and  y  except  x  - 

y  -  o. 

2  2  2  2 

Proof:  Since  (x  a  +y  b  )«(xa+jyb) (xa-jyb) ,  the  above  corollary  clearly  follows 
from  theorem  2.3.15  and  the  obvious  fact  that  a  polynomial  is  widest  sense 
Schur  if  and  only  if  its  factors  are  widest  sense  Schur. 

We  note  that  the  above  result  is  a  discrete  k-D  counterpart  of  a  result 
originaly  proved  by  Brockett  (11)  in  the  1-D  context.  The  k-D  continuous 
counterpart  of  corollary  2.3.16  is  obtained  in  (5],  however,  via  a  proof 
technique  entirely  different  from  the  one  presented  here. 
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2.4.  THE  REACTANCE  SCHUR  AND  IMUTTANCE  SCHUR  POLYNOMIALS: 

Definition  2.4.1:  A  selfparacon jugate  Schur  polynomial  a  is  said  to  be 
reactance  Schur  if  the  irreducible  factors  of  a  do  not  occur  with  multiplicity 
larger  than  one. 

The  following  results  immediately  follows  from  the  above  definition. 

Theorem  2.4.2:  (a)  Products  of  reactance  Schur  polynomials  that  are  pairwise 
relatively  prime  are  reactance  Schur  polynomials.  (b)  Any  factor  of  a 
reactance  Schur  polynomial  a  is  also  a  reactance  Schur  polynomial. 

Proof:  (a)  Follows  from  the  fact  that  the  product  of  selfparacon  jugate  Schur 
polynomials  is  also  selfparacon  jugate  Schur.  (b)  We  note  that  definition  2.4.1 
and  theorem  2. 2. 2. 8  imply  that  the  irreducible  factors  of  a  are 
selfparaconjugate  Schur,  and  are  thus  reactance  Schur.  The  result  then  follows 
by  appealing  to  part  (a)  of  the  present  theorem. 

Theorem  2.4.3:  (a)  The  numerator  and  denominator  of  a  discrete  reactance 
function,  written  in  irreducible  form,  are  always  reactance  Schur  polynomials, 
(b)  Conversely,  each  discrete  reactance  Schur  polynomial  b  is  the  denominator 
(and  consequently  the  numerator)  of  a  discrete  reactance  function  in 
irreducible  form. 

Proof:  (a)  If  C-b/a  is  discrete  reactance  function,  and  in  particular  a 
discrete  positive  function  in  irreducible  form,  then  due  to  theorem  2.3.5,  a  is 
a  widest  sense  Schur  polynomial  devoid  of  multiple  discrete  selfparaconjugate 
factors.  Furthermore,  since  a  is  a  widest  sense  Schur  it  cannot  contain  z.  as 

A  X 

a  factor  for  any  i.  Thus,  deg. a-deg.a  for  i  ■  1  to  k.  Since  from  assertion 
2.3.6  we  have  b/a  -  -  b/a,  in  irreducible  form,  we  must  have  a  *  ya,  where  y  is 
a  constant  i.e.,  a  is  selfparaconjugate  Schur.  Since  irreducible  factors  of  a, 
due  to  theorem  2. 2. 2. 8,  are  also  selfparaconjugate  Schur,  they  must  occur  in  a 
with  multiplicity  equal  to  one.  Therefore,  a  is  reactance  Schur.  Similar 
arguments  apply  to  b. 

(b)  Lot  b^  i  ■  l,2,...n  be  the  nontrivial  irreducible  factors  of  b,  which  are 
also  reactance  Schur  due  to  theorem  2.4.2.  Let  b^  involve  the  variable  z^j 


Since 


n 


and  consider 
from  theorem 


the  rational  function  4  -  £  4-»  where  4. 


- -  -  'i'  -  (bi>x(i)/bi  ’ 

2.3.12,  each  4^  is  a  discrete  reactance  function,  their  sum  *  is 


also  a  discrete  reactance  function.  Let  4  ■  c/b,  then  c  ■  £  ((b. ) 


i-1 


i'X(i) 


v  » .  1 

V*1 


b^).  It  then  follows  by  appealing  to  property  2.C2  in  the  appendix  that  c  and 


b  are  relatively  prime  polynomials. 


Theorem  2.4.4:  If  the  polynomial  a  is  nontrivial  and  is  scattering  Schur,  then 
the  discrete  paraeven  and  the  discrete  paraodd  parts  of  a  respectively  denoted 
by  a@-(a+S)/2  and  ao*(a-&)/2  are  relatively  prime  reactance  Schur  polynomials. 

Proof:  Since  a  is  scattering  Schur  the  rational  function  p»a/a,  due  to 
definition  2. 2. 3.1,  is  in  irreducible  form  and  satisfies  | p | <1  for  | z | <1 - 

A  A  “ 

Therefore,  if  C  -  (l-p)/(l+p)  ■  (a-a)/(a+a),  then  C  is  a  discrete  positive 

function  (i.e.  Ret  >_  0  for  |z|<^)  in  irreducible  form.  Consequently,  due  to 

theorem  2.3.5,  both  a„  -  a+a  and  a-  -  a  -  a  are  widest  sense  Schur.  The  facts 

e  u 

that  ae  and  aQ  are  respectively  paraeven  and  paraodd  follows  from  property 
2. BIO.  Therefore,  a@  and  aQ  are  selfparaconjugate  Schur.  Furthermore,  by 
appealing  to  the  second  part  of  theorem  2.3.5  it  follows  that  irreducible 
factors  of  ag  and  a^,  which  by  virtue  of  theorem  2. 2. 2. 8,  must  be 
selfparaconjugate  Schur,  must  not  occur  with  multiplicity  higher  than  one. 
Therefore,  aQ  and  a@  are  reactance  Schur  polynomials. 

It  is  now  possible  to  derive  the  converse  of  theorem  2.4.4. 

Theorem  2.4.5:  If  b  is  a  discrete  paraeven  (or  discrete  paraodd)  reactance 
Schur  polynomial  then  there  exists  a  discrete  paraodd  (or  discrete  paraeven) 
reactance  Schur  polynomial  c,  relatively  prime  with  b,  such  that  (b+-c)  is 
either  a  nonzero  constant  or  a  nontrivial  scattering  Schur  polynomial. 

Proof:  Consider  the  irreducible  discrete  reactance  function  C  -  c/b,  the 

existence  of  which  has  been  demonstrated  in  theorem  2.4.3(b).  Since  due  to 
assertion  2.3.6  we  have  (c/b)  +  ( c/b)  -0,  b  •  +  b  implies  c  -  ~  c.  Thus,  if  b 
is  paraeven  (or  paraodd)  then  c  is  paraodd  (or  paraeven).  Furthermore,  the 
rational  function  p  ■  (1-C)/(1+C)  -  (b-c)/(b4-c)  is  irreducible  form,  and 

satisfies  | p | < <1  for  |z|<l.  Consequently,  due  to  theorem  2.2.3.12(b),  (b+c)  is 
either  a  nonzero  constant  or  scattering  Schur  polynomial. 
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In  the  above  theorem  it  is  indeed  possible  to  have  b+c  -  constant.  Consider  b 
-  1  +  z^,  c  - 

The  following  is  an  yet  alternate  characterization  of  reactance  Schur 
polynomials. 


Theorem  2.4.7:  Let  a  be  a  se  1  f pa racon jugate  Schur  polynomial. 

reactance  Schur  if  and  only  if  the  set  of  (k+1)  polynomials  a,a_  ,a_ 

z 

not  have  a  proper  common  factor. 


2 


Then  a  is 
,  . . . a  do 

K 


Proof:  In  view  of  property  2.C2  in  the  appendix,  it  follows  that  the 
irreducidble  factors  of  a  are  simple  if  and  only  if  the  set  of  (k+1) 
polynomials  a,az  ,...az  do  not  have  a  proper  canmon  factor.  The  result  then 
follows  by  appealing  to  £he  definition  2.4.1  of  reactance  Schur  polynomials. 

Definition  2.4.8:  A  polynomial  a  is  said  to  be  immittance  Schur  if  it  is  the 
product  of  a  scattering  Schur  polynomial  and  a  reactance  Schur  polynomial. 

The  following  result  clearly  follows. 

Theorem  2.4.9:  (a)  Factors  of  immittance  Schur  polynomials  are  immittance 
Schur.  (b)  Products  of  immittance  schur  polynomials  that  do  not  have  a  discrete 
selfparaconjugate  factor  in  common  are  also  immittance  Schur.  (c)  The  least 
common  multiple  of  a  set  of  immittance  Schur  polynomials  is  also  immittance 
Schur. 


The  following  result  is  an  obvious  consequence  of  theorems  2.3.5. 

Theorem  2.4.10:  The  numerator  and  denomiantor  polynomials  a  and  b  of  a  discrete 
positive  function  C*a/b  in  irreducible  form  are  necessarily  immittance  Schur. 

Theorem  2.4.11:  If  a  is  a  polynomial  involving  z^  and  is  immittance  Schur  then 
the  polynomial  az  is  also  immittance  Schur. 

Proof:  Let  a_  /a  -  b/c,  with  a_  -  b.d  and  a  -  c.d,  where  d  is  the  gcd  between 

1  1  z »  z . 

az  .  Due  to  theorem  2.4.9,  d  isliramittance  Schur.  Since  the  polynomials  b  and 
c  lare  relatively  prime,  the  rational  function  b/c,  due  to  theorem  2.3.5  is  a 
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discrete  positive  function  in  irreducible  form.  Consequently,  by  virtue  of 
theorem  2.4.10,  b  is  an  immittance  Schur  polynomial.  Let  d  contain  a 
self pa racon jugate  Schur  factor  e.  Since  e  is  a  factor  of  a,  which  is 
immittance  Schur,  e  occurs  in  a  with  multiplicity  one.  Now,  if  e  involves  z ^ 
then  az  ,  by  property  2.C2,  does  not  contain  the  factor  e,  which  contradicts 
az  -  i.d.  Therefore,  d  cannot  contain  e,  and  in  particular,  b  and  d  do  not 
ha^e  discrete  selfparaconjugate  factors  involving  z^^  in  common.  On  the  other 
hand,  if  e  does  not  involve  z^  then  the  multiplicity  of  e  in  a  and  az  must  be 
the  same.  Also,  since  a»c.d  is  immittance  Schur,  c  cannot  contain  a  factor  e. 
Thus,  b  cannot  contain  the  factor  e,  beacuse  otherwise  the  multiplicity  of  e  in 
az  would  be  larger  than  in  a.  Therefore,  b  and  d  cannot  have  a  discrete 
seifparaconjuyate  factor  not  involving  z^  in  common.  Therefore,  b  and  d  cannot 
have  a  common  selfparaconjugate  Schur  factor,  thus  confirming  the  immittance 
Schur  property  of  a  -  b.d  in  view  of  theorem  2.4.9. 

Z  • 

1 

Theorem  2.4.12;  The  denominator  of  a  rational  discrete  positive  function  in 
irreducible  form  is  an  immittance  Schur  polynomial.  Conversely,  every 
immittance  Schur  polynomial  is  the  denominator  of  a  discrete  positive  function 
in  irreducible  form. 

Proof;  The  first  half  of  the  present  theorem  has  already  been  proved  in 
theorem  2.4.10.  To  prove  the  second  half  let  a  -  e.f,  where  e  and  f  are 
respectively  the  scattering  Schur  and  reactance  Schur  factors  of  a  (cf.  theroem 
2.2.6).  Notice  that  e  and  f  are  relatively  prime  due  to  theorem  2. 2. 2. 8  arid 
2. 2. 3. 2.  Since  e  is  scattering  Schur,  theorem  2.2.3.12(a)  yields  |e/e|<l  for 

A  A 

1 2! | <1  implying  Re[l+(e/e) ]>0  for  |zj<l  i.e.,  [l+(e/e)J  is  a  discrete  positive 
function.  Furthermore,  since  f  is  reactance  Schur,  by  virtue  of  theorem  2.4.3, 
there  exists  a  polynomial  g,  such  that  g/f  is  a  discrete  reactance  function 
(and  thus  a  positive  function)  in  irreducible  form.  Therefore,  l+(e/e)+(g/f )  - 

A  A 

p/a,  where  p  -  e(f+g)+ef  -  (e+e)f+eg  is  a  discrete  positive  function.  It 
remains  to  show  that  p  and  a«ef  are  relatively  prime  polynomials.  Clearly,  if 
e  and  p  had  a  factor  in  common,  then  e  and  (ef)  would  have  a  factor  in  common. 
However,  e,  being  scattering  Schur,  is  relatively  prime  with  e  ,  and  cannot 
have  a  factor  in  common  with  f  as  stated  earlier.  On  the  other  hand,  if  f  had 
a  factor  in  common  with  p  then  it  would  have  a  factor  in  common  with  eg.  Since 
g/f  is  an  irreducible  rational  function,  this  later  situation  implies  that  e 
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and  £  would  have  a  factor  in  common,  which  is  again  impossible.  Thus,  p/a  is 
an  irreducible  rational  function. 
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2.5.  EXTENSIONS  TO  PASSIVE  MULTIPORT  SCATTERING  AM)  imiTTANCE  FUNCTIONS: 


The  following  results  generalize  the  above  discussion  on  the  scattering  or 
immittance  description  of  k-D  discrete  passive  multidimensional  systems  to 
multiport  systems. 


Definition  2.5.1:  An  (NxN)  rational  matrix  H^H(z)  is  said  to  be  discrete 
bounded  if  H  is  holomorphic  in  |z|<l  and  (I^-HH  )  is  nonnegative  definite  at 
all  regular  points  of  H  in  the  domain  |z|<l. 

We  shall  need  the  following  extensions  of  lemma  2. 2. 3.^1  for  further 
developments  to  follows.  Let  us  consider  the  function  p«p(z)«  l  |p^(z)|,  where 
each  pi-bi(z)/ai(z)  with  a^-a^z),  b^-b^z)  polynomials  in  -Variables  z.  We 
note  that  since  p  is  not  an  analytic  function  of  z  neither  the  maximum  modulus 
principle  for  functions  of  several  complex  variables  [9]  nor  lemma  2.2.3.11  of 
the  present  paper  applies  to  p.  However,  the  following  result  is  true. 

Lemma  2.5.2:  Let  p-p(z)  be  as  defined  above.  If  ai»ai(z)  for  each  i-l,2...N  are 
widest  sense  Schur  polynomials  and  p<l  for  |z|-l,  whenever  a^O  for  each 
i«l,2...N,  then  |p|<1  for  |z|<l.  Furthermore,  if  p-1  for  some  z  in  |z  |-1  then 
p<1  for  |z|<l  unless  each  p^  is  constant. 

Proof:  Let  p(Zq)>1  with  |Zq|<1.  Consider  the  rational  functions  ( z)-C^p^(z) , 

i«l,2...N  where  each  C^  is  a  constant  such  that  |C^|*1,  and  r^(Zg)  is  a 

positive  real  number  i.e^,  I  Pj^ZQ^Mr^Zg)  l-r^z^) .  Also,  consider  the 

rational  function  R(z)«  E  r^(z).  Obviously,  then  |R(z)|<  E  |r^(z)|« 

E| p. <z) |«p(z)  for  all  z.  lA”^articular,  |R(z) |<p(z)<l  for  |z|-l,  wh4n£ver  a./O 
)  1  —  1 

tor  each  i-l,2...N.  Thus,  by  invoking  lemma  2.2.1.11  with  R(z)  it  follows  that 
either  |R(z)|<l  for  |z|<l  or  R(z)  is  a  constant.  The  first  of  the  last  two 
situations  is,N  however,  in  contradiction  with  the  fact  that 
1  Ipj^IZq)  |-p(Zo)>1*  On  the  other  hand,  if  R(z)  is  a  contant 
follow4“*that  1<|p(Zq) l-jR(z^) |-|R(z) |<p(z)  for  all  z.  However,  since 
|p(z)  ;>1  for  |z|«l  whenever  a^O  for  each  i,  the  last  stated  chain  of  weak 
inequalities,  in  fact  hold  with  inequalities  replaced  by  equalities.  Thus,  we 
have  that  p(z)«l  for  all  z.  Furthermore,  a  conqparison  of  the  equality 


R(Zq)-  E 
then  If1 
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N  N  N 

£  |r^(z)|-  l  p^(z)  |-|  p(z)  |-|R(z)  |  with  Z  r^(z)-R(z)  clearly  shows  that  r^z) 
ioi  each  i-i,2...N  is  a  real  positivi“Aumber  for  all  z.  It  then  follows  by 
invoking  a  standard  result  from  the  theory  of  function  of  several  complex 
variables  that  r^(z),  for  each  i,  is  a  constant  i.e.,  independent  of  z. 

The  following  result  can  be  viewed  as  a  multiport  extension  of  the  maximum 
modulus  theorem  for  function  of  several  coapex  variables  useful  in  the  present 
context  of  passive  or  lossless  networks. 

Theorem,.  2.5.3:  Let  H  be  an  (NxN)  rational  matrix  holomorphic  in  |z|<l  and 
-  £ 

(In~HH  )  is  nonnegative  definite  for  |z|«l,  wherever  H  is  well  defined.  Then  H 
is  a  discrete  bounded  matrix. 


Proof :  Let  x- 


(x1(x2, 


.x^)  be  any  (lxN)  constant  vector  and  • • -yk)  be 

a  (lxN)  vector  defined  ^ia  y-Hx.NThen  it  follows  frgm  the  nonnegative 

-  -  -  *  -  .2  -  - 


definiteness  of  )  that  £  |y.  |  <  £  |x.  |  ,  whenever  H  and  thus,  each 

i-1  1  i-1 

yi*Yi(z)  is  well  defined  on  |z|-l.  Clearly,  the  property  of  holomorphy  in  |z|<l_ 


is  inherited  by  each  y.  from  H.  Then  by  invoking  lemma  2.5.2,  it  follows  that 
N  ,  N  $ 

£  |y.  |z  <  £  | x .  }  for  | z | <1 .  Since  x  is  arbitrary  the  last  conclusion 

i-1  1  i-1  ~  ~ 

yields  that  (IN-HH  )  is  nonnegative  definite  for  |z|<l,  which  when  combined 
with  the  fact  that  ( IN_HH  )  is  nonnegative  definite  whenever  H  is  holomorphic 


in  I z|*l.,  yields,  in  view  of  definition  2.5.1,  that  H  is  a  discrete  bounded 
matrix. 


Definition  2.5.4:  A  rational  matrix  of  size  (NxN)  is  said  to  be  discrete  k-D 
lossless  bounded  if 

(a)  the  entries  -  H^(z)  of  H  -  H(z)  are  holomorphic  and 
unimodularly  bounded  i.e.,  |H^|  <  1  for  |zj  <  1. 

(b)  H  H  -  IN. 

we  note  that,  as  expected,  the  class  of  multidimensional  discrete  lossless 
bounded  matrices  form  a  subclass  of  the  class  of  multidimensional  discrete 
bounded  matrices. 
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Property  2.5.5:  The  denominators  of  each  entry  of  discrete  k-D  lossless  bounded 
matrix  H,  when  written  in  irreducible  form,  are  scattering  Schur  polynomials. 

Proof:  Follows  from  theorem  2.2.3.12(b). 

Property  2.5.6:  Determinant  of  H  can  be  expressed  in  irreducible  form  as  det  H 
■  d(a/4),  where  a  is  a  scattering  Schur  polynomial,  and  d  is  a  unimodular 
constant  i.e.,  |d|  -  1. 


Proof:  Let  det  H  -  a/h  ,  where  a  and  h  are  relatively  prime  polynomials. 
Since  due  to  theorem  2. 2. 3. 9,  products  as  well  as  factors  of  scattering  Schur 
polynomials  are  also  scattering  Schur,  from  property  2.5.5  it  follows  that  h  is 
a  scattering  Schur  polynomial.  Now,  from  definition  2.5.4(b)  of  lossless 
bounded  matrices  it  obviously  follows  that  det  H  det  H  -  1  i.e.,  (a/h)(3/Ti)* 
1,  thus  implying  (2.5.1a)  below.  Since  in  (2.5.1a),  the  polynomial  h  is 
relatively  prime  with  a,  and  h,  being  scattering  Schur,  cannot  have  a  factor  z ^ 
for  any  i,  the  polynomial  4  must  contain  h  as  a  factor  i.e.,  (2.5.1b)  follows, 
where  d  is  a  polynomial  in  z. 


n.  -  na 

a  4  z  -  hh  z  ,  4  -  h.d 


(2.5.1a,b) 


Substituting  (2.5.1b)  in  (2.5.1a)  and  considering  the  discrete  paraconjugate 
of  the  resulting  equation  it  follows  from  Property  2.A5  that  ad-h.  The  latter 
equality,  in  view  of  scattering  Schur  property  of  h  and  property  2.A6,  imply 
(2.5.2). 


a  .  d  -  h  (2.5.2) 

A 

Substituting  for  h  from  (2.5.2)  into  (2.5.1b),  we  obtain  dd*l,  which  in  turn 

*  *  if 

imply  that  d  is  a  constant.  Consequently,  d«d  ,  and  dd  »|d|»l.  It  then  follows 
from  (2.5.2)  that  det  h-  a/h  -  a/(ad)  -  d(a/4). 

Remark  2.5.7:  Although  a  is  scattering  Schur  a  is  not  necessarily  so. 

Consider  H  -  z. 


Property  2.5.8:  A  discrete  lossless  k-D  rational  matrix  H  can  be  written  as  H 
-  P/4,  where  P  ■  P(z)  is  a  polynomial  matrix  in  z,  and  a  is  as  in  property 
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2.5.6  above. 


Proof ;  Fran  condition  (b)  of  the  definition  2.5.4  and  property  2.5.6  it  follows 
that  H  -  H  ^  »  d  (Adj  H)(S/a),  which  in  turn  imply  that  H  -  H  -  d(S/a) (AdjH) . 
Furthermore,  it  is  straightforward  to  verify  that  3/5  -  a/S  for  any  polynomial 
a.  Hie  last  two  equlities  combined  togather  yield  (2.5.3)  in  the  following. 

S.H  ■  d.a. (Adj  H)  (2.5.3) 

Due  to  property  2.5.5  and  theorem  2.2.3.9(b)  the  entries  of  Adj  H,  when 
written  in  irreducible  form,  have  scattering  Schur  denominators,  and  therefore, 
are  holomorphic  in  |z|<l.  Hence  in  (2.5.3)  the  right  hand  side  is  holomorphic 
in  |2|>1.  The  left  hand  side  is  either  a  polynomial  matrix  P*P(z)  or  the 
entries,  due  to  property  2.5.5  and  theorem  2.2.3.9(a),  have  scattering  Schur 
denominators  when  written  in  irreducible  form.  In  the  latter  case,  however, 
SH,  by  virtue  of  theorem  2. 2. 3. 6,  is  not  holomorphic  in  |£|>1.  Thus,  aH-P, 
where  P-P(z)  is  a  polynomial  matrix. 


For  N*2  by  considering  the  fact  that  H  -  H-1,  (2.5.4)  and  (2.5.5)  below 
follows  from  property  2.5.6. 


Hn  -  d*(S/a)H22  ;  H21  -  -  d*(S/a)H12  (2.5.4a,b) 

H12  -  -  d*(S/a)Hn  ;  H22  -  d*(S/a)H11  (2.5.5a,b) 


We  note  that  (2.5.4a)  and  (2.5.4b)  are  respectively  identical  with  (2.5.5b)  and 
(2.5.5a)  via  the  operation  Also,  H^2  in  (2.5.5a)  is  holomorphic  in  |z|>l_, 
whereas  H21,  due  to  property  2.5.5,  must  have  scattering  Schur  denominator,  and 
therefore,  due  to  theorem  2. 2. 3. 6  has  singularities  in  (z|>l.  Hence  the 
denominator  of  H^  divides  3,  i.e.,  H2^  -  c/3  (not  necessarily  in  irreducible 
form).  Therefore^  from  (2.5.5a)  it  follows  that  H^2  -  -  d*(c/a)  i.e.,  H12  *  - 
d(2/5)  ■  -  It  follows  from  (2.5.5b)  via  a  similar  argument  that  H^  « 
b/3  (not  necessarily  in  irreducible  form),  and  H22  -  d(b/3)z_a.  Thus,  H  can  be 
written  as  in  (2.5.6)  below. 
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~  Ha 
-dc  2 


+db  z-* 


(2.5.6) 


Also,  by  considering  the  determinant  of  (2.5.6),  it  follows  via  the  use  of 
property  2.5.6  that 


bb  +  cc  -  aa 


(2.5.7) 


-  n 

Furthermore,  since  H^2  ”  -  dcz 
it  is  necessary  that  deg, c  <  n  . 

1  al 


-  -  n  - 

/a  and  H22  *  dbz  /a  are  holomorphic  in  | z | <1, 
and  deg . b  <  n  .  for  i«  1  to  k. 

1  al 


The  following  result  can  be  considered  to  be  the  discrete  multidimensional 
counterpart  of  the  Belevitch  canonical  form  for  the  representation  of  lossless 
two-port  networks  well  known  in  classical  network  theory.  The  multidimensional 
continuous  version  of  this  result  has  recently  been  discussed  in  [2]. 


Theorem  2.5.9:  A  (2x2)  rational  matrix  H  -  H(z)  is  a  lossless  bounded  matrix 
if  and  only  if  there  exists  polynomials  a,b,c  and  a  unimodular  complex  constant 
d  such  that: 


(i)  H  can  be  written  as  (2.5.6) 

(ii)  5  is  a  scattering  Schur  polynomial 

(iii)  a,b,c  are  related  as  in  (2.5.7) 

(iv)  deg^^deg^b,  and  degia>degic  for  all  i-1  to  k. 

Proof:  Only  the  sufficiency  part  of  the  theorem  needs  to  be  proved,  the 
neceesity  being  already  established.  We  need  to  show  that  under  the  conditions 
stated  in  the  present  theorem,  condition  (a)  and  (b)  of  definition  2.5.4  are 
satisfied.  Straightforward  computation  along  with  (2.5.7)  and  the  identitity  55 
■  3a  yields  HH  -  I2,  which  verifies  condition  (a)  of  definition  2.5.4. 

Furthermore,  a  second  use  of  the  identity  55  -  aa  along  with  (2.5.7)  yields 
(2.5.8)  below. 

( b/S ) ( B/S )  +  ( c/S ) ( 5/5 )  -  1,  ( b/S ) ( B/S )  +  (c/S)(c/S)  -  1  (2.5.8a,b) 
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It  clearly  follows  from  (2.5.8b)  that  |b/8|^  +  |c/8|^  -  1,  and  thus  |b/S|<l, 

| c/S | <^j,  for  |zi*l^  whenever  &  /  0.  However,  since  &  is  scattering  Schur,  both 
(b/ajz-3  and  (a/ajz-3  are  holomorphic  in 

|z|  <1,  and  t  as  an  application  of  lerana  2.2.3.11  yields  that  |H12I<1  and 

I H22 1  — 1  ^or  I— I  <  s^^ar  arguments  with  (2.5.8a)  yield  |H^|<1  and 
thus  completing  the  proof  of  the  present  theorem. 


Definition  2.5.10:  An  (NxN)  rational  matrix  Z  -  Z(z)  is  said  to  be  discrete 

positive  if  each  entry  Z.  .(z)  is  holomorphic  in  |z|<l  ,  and 
*  t  ^  J 

Z(z)  +  Z  (z)  is  nonnegative  definite  for  |z|<l. 

Lemna  2.5.11:  Let  Z-Z(z)  be  a  discrete  positivefcmatrix  in  k-variables.  If  Z(z) 
is  regular  at  z^  with  |Zg|<l  then  Z( z^)+(Z( z^) )  is  nonnegative  definite. 

Proof:  If  |Zq|<1  then  the  result  is  obvious  in  view  of  definition  2.5.10.  If 

I but  | Zq Ifl  then  consider  the  neighbourhood  Nc-{z:  |  (z-z^ j  | <c}1of  Zg  such 

that  Z(z)  is  holomorphic  in  N  .  The  existence  of  such  a  N  follows  from 

property  2.A5  in  the  appendix.  Thus,  in  particular,  each  element  of  the 

matrices  Z(z)  and  Z  (z)  are  continuous  functions  of  z  in  the  connected  set 

t 

N-N  n{z: |z|cl} .  Consequently,  for  any  constant  vector  x,  x(Z+Z  )x  is  a 
*  ¥ 
continuous  faction  of  z  in  N.  Furthermore,  x(z+z  )x  >0  for  zeN  due  to 

-  t 

definition  2.5.10.  It  then  follows  by  invoking  the  continuity  of  x(Z+Z  )x  in 
t, 

N  that  x(Z+Z  )x  >0  for  z~Zq  and  for  arbitrary  x. 

The  following  result  is  a  generalization  of  theorem  2.3.5  to  multiports. 
Property  2.5.12:  Let  Z(z)-[Z^ j(z) ]  be  a  discrete  positive  matrix.  Then  the 

denominators  of  each  entry  Z^(z),  when  expressed  in  irreducible  form,  are 
immittance  Schur  polynomials. 

Proof:  The  fact  that  the  denominators  of  each  entry  of  Z.,  j(z)  is  widest  sense 

Schur  obviously  follows  from  definition  2.5.10.  Let  Z'fz^)  -  Zfz^z^),  where  Zg 
is  such  that  |z'g(-l  and  let  a»a(z)  be  the  l.c.m  of  the  denominator  polynomials 

1  ^  2 
||.||  denotes  the  Eucledean  norm  i.e.,  ||z||-  £  |z.  |  . 

i-1  1 
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of  the  elements  of  the  matrix  Z(z).  Let  v  be  the  multiplicity  of  a 
selfparaconjugate  factor  c  in  the  denominator  of  Z^(z)  in  irreducible  form. 
Then  due  to  property  2.A4,  c'-c'U^-cU^z^)  will  be  a  factor  of  multiplicity 
v  in  the  denominator  of  Zf j(z^)-z! j(zlfz^)  in  irreducible  form  for  any  choice 
of  z'0  on  |z'|-l  from  a  sequentially  almost  complete  set  a  of  (k-1)  tuples  of 
order  (k-1).  We  note  the  polynomial  a' -a' (Zj)-a(z^,z'g)  cannot  be  identically 
zero  for  all  z^  in  Q,  beacuse  otherwise  a  would  be  identically  zero  due  to 
property  2.A2  in  the  appendix.  Consequently,  there  exists  z^  c  Q  with  |Zq|«1 
such  that  Zj(z^)  isfc  well  defined.  By  invoking  lemma  2.5.11  it  then  follows 
that  Z'(z^)+(Z'(z^))  is  a  discrete  positive  matrix  function  in  the  variable 
z ^  only.  Since  c'  is  a  selfparaconjugate  factor  in  the  denominator  of  Z?  ^  it 
immediately  follows  from  known  one-variable  results  that  v  -  1. 


Corollary  2.5.13: 
common  denominator 
polynomial  i.e., 
polynomial  matrix. 


If  Z-Z(z)  is  a  discrete  positive  matrix  then  the  ast 
of  all  entries  Z^(z)  of  Z(^)  is  an  immittance  Schur 
there  exists  an  inmittance  Schur  polynomial  a  such  that  aZ  « 


Proof;  Follows  from  the  fact  that  l.c.m.  of  a  set  of  immittance  Schur 
polynomials  is  also  an  immittance  Schur  polynomial  (theorem  2.4.9c). 
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2.6.  RELATIONSHIP  BETWEEN  MULTIDIMENSIONAL  SCHUR  POLYNOMIALS  AND 
MULTIDIMENSIONAL  HUIWIT2  POLYNOMIALS: 


Let  g(p)  and  a(z)  be  polynomials  in  the  variables  p  and  z  respectively.  Also 
let  the  polynomial  a(z)  be  obtained  from  g(p)  via  the  bilinear  transformations 
p^  -  (l-zi)/(l+zi)  for  i-1  to  k  as: 


1— z.  l-z7 


k 

n  d+2. ) 
i-l  1 


m. 

l 


(2.6.1) 


where,  nu  -  deg^  for  i  «  1  to  k. 

Then  the  polynomial  a  -  a(z)  has  been  called  [7]  the  associate  of  the 
polynomial  g  -  9<£)-  On  the  otherhand  if  a  and  g  are  related  as  in  (2.6.2) 
below  then  g  is  called  [7]  the  associate  of  a. 


9(E) 


i-Pi  i-p2 

a,KpI'I^' 


Ifk 

I+pk 


k  n. 

n  (1+p. ) 

i-1  1 


(2.6.2) 


where  deg^  -  n^  for  i  -  1  to  k. 

General  properties  of  associates  of  polynomials  in  z  and  in  p  variables  and 

their  interrelations  are  discussed  extensively  in  (7).  We  only  note  here  that 

two  polynomials  may  not  be  associates  of  each 

other.  Consider,  for  exmple,  a  -  FI  (1+z.)  and  g  -  2k. 

1-1  1 

In  this  context  it  has  been  shown  in  [7]  that  any  scattering  Schur  polynomial 

is  the  associate  of  a  scattering  Hurwitz  polynomial.  Conversely,  any 

scattering  Hurwitz  polynomial  devoid  of  factors  (1+p^)  for  any  i  -  1  to  k,  is 
also  associate  of  some  scattering  Schur  polynomial.  Stated  alternatively,  the 
bileaner  transformation  along  with  its  inverse  transformation  sets  up  a 

one-to-one  correspondence  between  the  scattering  Schur  polynomials  and 

scattering  Hurwitz  polynomials  devoid  of  factors  of  the  type  (1+p.).  The 

scattering  Hurwitz  polynomials  containing  factors  of  the  type  (1+p^),  when 

bilinearly  transformed  via  pi  -  (l-z^/d+z^  also  yield  scattering  Schur 

polynomials.  This  is  in  contrast  to  the  situation  discussed  in  (4),  where  it 
is  shown  that  the  associates  of  strict  sense  Hurwitz  polynomials  need  not  be 
strict  sense  Schur.  In  what  follows  answers  to  questions  such  these  are  sought 
in  the  context  of  other  classes  of  multidimensional  polynomials. 
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We  first  note  that  the  associate  of  a  discrete  selfparaconjugate  polynomial  is 
necessarily  selfparaconjugate  i.e.  if  a  -  ya  then  g*  -  yg,  but  the  converse  is 
not  true  i.e.,  the  associate  of  selfparaconjugate  polynomial  is  not  necessarily 
discrete  sglfpapcon jugate.  However,  it  is  true  that  if  g*  -  yg  then 
a  -  y.a  H  zi  1 ,  where  r^  -  deg^  -  deg  .a  -  multiplicity  of  the  factor  (1+p^ 
in  g.  H&wiver,  if  g  is  selfparaconjugate  Schur  then  its  irreducible  factors 
must  also  be  so,  and  therefore,  g  cannot  have  a  factor  (1+p.).  Thus,  if  g  is 
selfparaconjugate  then  its  associate  a  is  also  discrete  selfparaconjugate. 
Furthermore,  the  associate  of  a  widest  sense  Hurwitz  polynomial  is  widest  sense 
Schur,  and  conversely  the  associate  of  a  widest  sense  Schur  polynomial  is 
widest  sense  Hurwitz  (7).  It,  therefore,  follows  that  there  is  a  one-to-one 
correspondence  under  the  bilinear  transformation  and  its  inverse  between  the 
selfparaconjugate  Schur  polynomials  devoid  of  factors  of  the  type  (1+z^and  the 
selfparaconjugate  Hurwitz  polynomials.  Furthermore,  selfparaconjugate  Schur 
polynomials  containing  factor  of  the  type  (l+zi>,  under  the  action  of  the 
transformation  pi  -  (1-z^ )/(l+z^ )  i-1  to  k  necessarily  yield  selfparaconjugate 
Hurwitz  polynomials.  Since  it  is  known  [7,  lemma  9]  that  if  two  polynomials 
are  relatively  prime  their  associates  are  also  relatively  prime  the  above 
comments  also  apply,  in  particular,  to  reactance  Hurwitz  polynomials. 

Furthermore,  since  the  immittance  Hurwitz  (Schur)  polynomials  are  products  of 
reactance  Hurwitz  (Schur)  polynomials  and  scattering  Hurwitz  (Schur) 
polynomials  the  following  conclusions  follow.  Hie  whole  class  of  immittance 
Schur  polynomials  devoid  of  factors  of  the  type  (1+z^)  can  be  identified  as  the 
class  of  associates  of  the  whole  class  of  immittance  Hurwitz  polynomials.  On 
the  other  hand,  the  whole  class  of  immittance  Hurwitz  polynomials  devoid  of 

factors  (1+p^)  for  i  ■  1  to  k  can  be  identified  as  the  class  of  associates  of 

the  class  of  immittance  Schur  polynomials.  Furthermore,  if  a  is  any  immittance 
Schur  polynomial  not  containing  (1+z^  for  any  i  -  1  to  k,  as  a  factor  then 
there  exists  an  immittance  Hurwitz  polynomial  g  without  the  factor  (1+p.)  such 
that  a  and  g  are  associates  of  each  other.  Conversely,  for  any  immittance 

Hurwitz  polynomial  g  devoid  of  the  factor  (1+p^)  for  all  i  ■  1  to  k  there 

exists  an  immittance  Schur  polynomial  a  such  that  g  and  a  are  associates  of 
each  other.  Thus,  there  is  a  one-to-one  correspondence  between  the  members  of 
the  class  of  immittance  Schur  polynomials  devoid  of  factors  of  the  type  (1+z.) 
for  i  -  1  to  k  and  the  members  of  the  class  of  immittance  Hurwitz  polynomials 
devoid  of  factors  of  the  type  (1+p^)  for  i  -  1  to  k. 
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The  above  discussion  can  be  summarized  in  the  form  of  theorem  2.6.1  for  the 
purpose  of  which  the  following  concise  notation  is  adopted.  The  set  of 
scattering  Hurwitz  (Schur ) ,  selfparacon jugate  Hurwitz  (Schur),  reactance 
Hurwitz  (Schur)  and  iramittance  Hurwitz  (Schur)  polynomials  are  denoted 
respectively  by  SH  (SS),  PH  (PS),  RH  (RS)  and  1H  (IS).  Also,  the  subclass  of  a 
certain  class  of  polynomials  devoid  of  factors  of  the  type  (1+p^  or  (1+z^  for 
all  i,  is  to  be  denoted  by  adding  the  subscript  l^p  or  1+z,  as  the  case  may  be, 
to  the  symbol  designating  the  corresponding  class  of  polynomials.  For  example, 
the  class  of  scattering  Hurwitz  polynomials  devoid  of  factors  of  the  type  l+pif 
for  all  i  will  be  denoted  by  SH1+p,  whereas  the  class  of  selfparacon jugate 
Schur  polynomials  deviod  of  factors  of  the  type  is  to  be  denoted  by  K>1+z* 
Furthermore,  the  notation  A(C)  will  denote  the  the  set  of  polynomials  obtained 
by  considering  the  associates  of  all  polynomials  belonging  to  a  class  C. 


Theorem  2.6.1:  The  following  set  inclusion  relations  hold  true. 


Theorem  2.6.1:  ' 

(i) 

A(SHi+p) 

and  those 

(ii) 

A( PH)  -  ; 

(iii) 

and  PS1+z 
A(RH)  -  1 

(iv) 

and  RS1+Z 
A(IHl+p>  ' 

A(  SH)  -  SS 


1+z 


SS;  A(SS)  -  SHi+p  c  SH*  T*1e  elements  of  SH^+p 


PSi+z  c  PS;  A(PS1+z) 


RSi+z  c  RS;  A(RS1+Z) 


-  A (PS)  -  PH  -  PH 


A(RS)  -  RH 


^l 


1+p.  The  elements  of  Pft^+p 
The  elements  of  RH, 


+P 


1+P 


The  elements  of  IH^+p  and  those  of  IS^+z  are  in  1-1  correspondence. 
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2.7.  SWWAKY: 


Properties  of  widest  sense  Schur,  selfparacon jugate  Schur  and  scattering  Schur 
polynomials  previously  not  published  in  the  literature  are  discussed. 
Elementary  properties  of  discrete  multidimensional  positive  functions  are  then 
studied.  Hie  reactance  Schur  and  the  immittance  Schur  polynomials  occurring 
respectively  as  the  numerators  and  denominators  of  discrete  reactance  functions 
and  discrete  positive  functions  are  introduced  for  the  first  time,  various 
alternate  characterizations  of  these  polynomials  are  then  suggested.  The 
heirarchial  ralationship  between  these  different  classes  of  Schur  polynomials 
can  be  summarized  diagraroatically  as  is  done,  for  example,  in  the  case  of  the 
corresponding  classes  of  Hurvitz  polynomials  in  (2).  The  role  of  these 
polynomials  in  scattering  or  immittance  descriptions  of  discrete  (k-D)  passive 
multiports  are  investigated.  In  particular,  a  discrete  (k-D)  counterpart  of 
Belevitch  canonical  form  well  known  in  classical  network  theory  having 
potential  applications  in  passive  multidimensional  digital  filter  design  is 
derived.  Relationship  between  the  classes  of  multidimensional  Schur  polynomilas 
and  the  corresponding  classes  of  multidimensional  Hurwitz  polynomials  are  also 
discussed  within  the  present  context,  and  are  shown  diagrammatically  in  fig 
2.1. 
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APPENDICES 


2. A.  The  proofs  of  the  results  stated  below  are  ommitted  for  brevity  (2). 


Definition  2.A1:  Let  P  be  a  set  of  k-tuples,  z  -  (z^,z2, . . .z^) ,  where  all  z ^ 
belong  to  a  same  number  field  8  (e.g., either  the  field  of  real  numbers  or  the 
field  of  complex  numbers).  We  designate  by  8^,  i  ■  1  to  k,  subsets  of  8  and  we 
say  that  8^  is  almost  complete  if  it  comprises  almost  all  elements  of  8  (  i.e., 
all  elements  except  at  most  finitely  many),  we  then  say  that  P  is  a 
sequentially  almost  complete  set  of  order  m  >  1,  with  m  <  k,  if  there  exists  a 
permutation  i.,^,...^  of  the  integers  1,2,...  k  such  that  all  z  c  P  can  be 
generated  in  the  following  way.  There  exists  an  almost  complete  8^  such  that 
any  z^  e  8^  may  be  chosen.  For  any  choice  thus  made,  assuming  m  >  2,  there 
there  eiists  an1almost  conplete  8^^  (possibly  depending  on  the  particular  z^^  e 
8^  selected)  such  that  any  z^  c28^  may  be  chosen.  Again  for  any  choice  tAus 
maAe,  assuming  m  >  3,  there  exists  a A  almost  complete  8.  (  possibly  depending 
on  the  particular  z i  and  z ^  selected)  such  that  any  z.^  e  8.^  may  be  chosen 
etc.  If  m  «  k  this  process  Is  continued  until  we  have  rlached3p.  .  If  m  <  k, 
once  we  have  reached  i  there  exists  at  least  one  (k-m)  tuple  (z^  , ...z^  ) 
(possibly  depending  on  the  particular  zi  to  z^  selected)  that  may1^  choseK. 
Finally,  we  may  extend  the  above  definition  t8  the  situation  m  -  0  by  saying 
that  in  this  case  P  is  not  empty. 


The  set  P  is  said  to  be  sequentially  infinite  if  in  the  above  definition  the 
term  almost  complete  is  replaced  by  infinite. 

Theorem  2  .A 2;  If  a  is  a  polynomial  in  k-variables  such  that  the  set  of  zeros  of 
a  comprises  a  sequentially  infinite  set  of  order  k  then  a  is  identically  equal 
to  zero. 


Theorem  2. A3:  If  and  b  are  polynomials  in  k-variables  then  a  and  b  have  a 
proper  common  factor  if  and  only  if  the  set  of  common  zeros  that  are  common  to 
a  and  b  is  a  sequentially  infinite  set  of  order  (k-1). 

Theorem  2 . A4 :  Let  a  nd  b  be  two  relatively  prime  polynomials.  For  any  m  such 
that  l<nKk  let  us  freeze  m  of  the  variables  z^,  say  for  i^  to  iffl  at 


corresponding  values  ziQ.  Let 
remaining  variables.  Then  there  exists  a  sequentially  almost  complete  set  2ffl  of 


^  and  b^  be  the  resulting  poynomials  in  the 

st  complete  set  Q  of 
m 

the  polynomial  a^  and 


m-tuples  of  order  m  such  thatfor  (zi  ,z ^  , ...z^  )c2m*  the  polynomial  a1  an 

b ^  are  still  relatively  prime.  Furthe^&ore^any  orS2ring  may  be  chosen  for  2^. 


Theorem  2.A5:  If  a(z)  is  a  polynomial  and  a(z^)^0  for  some  Zq  then  there  exists 
c>0  such  that  a(z)/0  for  all  z  in  the  neighbourhood  ||z-z0||<e,  where  ||.|| 
denotes  the  Euclidean  norm. 


2.B.  Various  properties  of  the  operations  and  as  defined  in  the 

introductory  section  of  the  present  paper,  are  derived  in  this  section. 

Property  2.B1:  If  a  -  b.c  then  a-bc  and  a  -  b.c,  where  b  and  c  are 

polynomials. 

Property  2.B2:  If  a  is  any  polynomial  then  a  cannot  contain  the  factor  zi  for 
any  i. 

Property  2.B3:  If  the  polynomial  contains  the  factor  zi  with  exact  multiplicity 
Pj^  then  (degia-degia)  -  p^ 

Property  2.B4:  If  the  polynomial  a  does  not  have  z,.  for  any  i  as  a  factor,  and 
b  -  a  then  b  -  a. 

Property  2.B5;  If  a  is  a  widest  sense  Schur  polynomial,  and  b  -  a  then  b  -  a. 

Property  2.B6:  If  the  polynomial  d  is  the  greatest  common  factor  between  the 
polynomial  a  and  its  discrete  paraconjugate  a  then  d  is  a  discrete 
selfparacon jugate  polynomial. 

Property  2.B7:  Let  a'«a'(£')  be  the  polynomial  obtained  from  by  freezing  in  the 

polynomial  a-a(z)  of  the  variables,  say  z^,  at  z1*z1q  on  |z1|-l.  Then 

(aj_  _  «  r.i  a'  n  z.  x,  where  p.«deg.a  -  deg.a'X)  for  all  i^O. 

Z1  z10  iU  i-2  1  11  1 

Property  2.B8;  Let  a  be  a  discrete  selfparaconjugate  polynomial  with  a-ya 
involving  the  variable  z If  a  is  expressed  as  in  (2.1.1)  then  for 
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\*-0,l, . .  .n^  (2.A1)  holds  true. 

n 

ra^  -  (an._v)nzr  *  where  the  product  extends  over  all  t/i  (2.A1) 

Property  2.B9:  If  a  is  a  polynomial  and  deg^(a+a)  <  deg  .a  for  some  i,  then 
(a+a)  must  contain  the  factor  z i. 

Property  2. BIO:  If  the  polynomials  a  and  b  are  such  that  b  •  a+a  is  widest 
sense  Schur  then  b  -  +b,  i.e.,  b  is  either  paraeven  or  paraodd  respectively. 

2.C.  The  polynomial  az  associated  with  a  is  obtained  from  a  via  a  formal 
algebraic  operation  as  defined  in  section  2.1.  It  turns  cut  that  this  operation 
plays  the  corresponding  role  of  considering  the  derivative  of  a  polynomial  in 
the  context  of  continuous  systems.  Related  properties  of  a  are  studied  in  the 
following. 

Property  2. Cl:  If  f  -  gh,  where  g  and  h  are  polynomials  in  z  with  their 

coefficients  as  polynomials  over  the  field  of  complex  numbers  then 
fz  -  hgz  +  ghz. 

Property  2.C2:  If  f  contains  a  factor  g  of  multiplicity  n  (n>l)  then  f_ 

Z 

contains  the  factor  g  with  exact  multiplicity  (n-1). 

Property  2.C3:  If  the  coefficients  of  the  polynomial  f  belong  to  the  field  of 
complex  numbers  and  f#0  for  |z|<l,  then  Re(fz/f)>0  for  |z|<l. 
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CHAPTER  3 

ON  A  GENERALIZED  FACTORIZATION  PROBLEM  FOR  THE  SYNTHESIS 
OF  QUARTER  PLANE  TYPE  MULTIDIMOBICNAL  DIGITAL  FILTERS 


3.1.  INTRODUCTION: 

Various  synthesis  schemes  such  as  the  Darlington  synthesis  scheme  for 
synthesizing  lossless  transfer  functions  as  a  cascade  interconnection  of 
elementary  lossless  building  blocks  such  as  inductors,  capacitors,  gyrators 
etc.  in  the  continuous  domain  are  well  known  in  classical  network  theory.  The 
corresponding  problem  in  the  discrete  domain,  namely  that  of  synthesizing  a 
discrete  lossless  bounded  (or  positive)  transfer  function  as  a  structurally 
passive  interconnection  of  elementary  lossless  building  blocks  was  first 
resolved  via  transformation  from  prototype  problems  in  the  continuous  domain, 
and  the  resulting  class  of  filter  structures  are  now  known  as  the  wave  digital 
filters  (1),  (17].  Recently,  however,  successful  attempts  to  derive  these  and 
similar  other  discrete  domain  results  without  making  explicit  use  of  tools  of 
classical  network  theory  have  been  made.  Notable  among  these  are  the  orthogonal 
filters  [2], [4], [14], [15]  and  the  class  of  filters  referred  to  as  the  lossless 
bounded  real  (LBR)  filters  described  in  [3],  [13]  and  in  related  other 
publications. 

In  view  of  interest  in  the  synthesis  of  multidimensional  (k-D)  wave  digital 
filters,  the  problem  of  synthesis  of  k-D  lossless  two-port  scattering  transfer 
matrix  via  the  bisection  of  a  prescribed  two-port  into  a  cascade  connection  of 
two  lossless  two-port  sections  of  smaller  "degree"  has  been  addressed  in  the 
continuous  domain  in  [5].  Factorability  of  continuous  domain  two-port 
scattering  matrices  has  also  been  studied  recently  [11]  in  the  multidimensional 
context.  An  attempt  to  develop  a  complete  and  self  consistent  theory  for  the 
synthesis  of  k-D  structurally  passive  quarter-plane  causal  type  digital  filters 
independent  of  the  continuous  domain  methods  have  already  been  initiated  in 
[8],  [9]  and  [10]  by  discussing  the  discrete  domain  stability  properties  of  a 
class  of  multidimensional  polynomials.  In  the  present  paper  the  problem  of 
synthesizing  a  k-D  discrete  quarter-plane  causal  type  lossless  two-port  as  a 
structurally  passive  interconnection  of  more  elementary  digital  building  blocks 
directly  in  the  discrete  domain  is  approached  by  the  methods  of  factoring  the 
chain  matrix,  the  hybrid  matrix  and  the  transfer  function  matrix  associated 
with  a  lossless  two-port.  By  following  recent  results  in  [9]  it  can  be  shown 
that  each  of  these  matrices  can  be  uniquely  expressed  by  means  of  a  set  of 
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three  polynomials  in  a  fora  analogous  to  the  Belevitch  canonical  fora  [16]  of 
classical  circuit  theory.  For  the  purpose  of  unified  presentation  of  results,  a 
matrix  referred  to  as  the  (multidimensional)  generalized  lossless  two-port 
matrix,  which  can  be  viewed  as  a  generalization  of  the  multidimensional  chain 
matrix,  the  hybrid  matrix  and  the  transfer  function  matrix  has  been  introduced. 
Interestingly,  in  1-D  this  matrix  can  be  categorized  under  the  class  of  sigma 
lossless  rational  matrices  considered,  for  example,  in  [12].  The  problem  of 
factorizing  this  generalized  lossless  two-port  matrix  into  the  product  of  two 
matrices  of  identical  type  can  then  be  viewed  as  a  problem  of  structurally 
passive  synthesis  of  multidimensional  two-ports.  It  must  be  noted  that 
factorization  of  chain  matrix,  when  feasible,  yield  networks  having  cascade 
structures  as  shown  in  Figure  3.1,  whereas  the  factorization  of  hybrid  matrix 
and  transfer  function  matrix,  when  feasible,  yield  networks  having  the 
topological  structure  as  shown  in  Figures  3.2  and  3.3. 

Necessary  and  sufficient  conditions  for  this  generalized  factorization  problem 
so  introduced  to  be  solvable  are  obtained  in  the  present  paper  via  constructive 
techniques.  As  expected  from  our  previous  study  of  analogous  problems  in  the 
continuous  case  [5], [11]  it  turns  out  that  the  factorization  may  not  be 
feasible  in  a  generic  multidimensional  (k>2)  situation.  However,  the 
impossibility  of  factorization  of  the  chain  matrix  does  not  by  any  means  rule 
out  the  feasibilty  of  factorization  of  the  transfer  function  matrix.  Exactly 
same  comment  also  applies  if  the  role  of  three  types  of  matrices  (i.e.,  chain, 
hybrid  and  transfer  function)  are  permuted  in  any  posssible  manner  (cf.  Section 
7).  Furthermore,  in  order  for  the  factorizations  under  consideration  to  yield 
computable  digital  filters  the  structures  resulting  from  the  factorization  may 
not  have  any  delay  free  loop.  Apparently,  this  inposes  a  further  constraint  on 
the  factorization  not  present  in  the  corresponding  continuous  domain  problems 
discussed  in  [5]  and  [11].  However,  as  shown  in  Section  3.5,  this  constrained 
problem  cam  always  be  solved  if  and  only  if  a  solution  to  the  unconstrained 
problem  exists. 

In  the  special  case  of  1-D,  the  criterion  for  factorability  is  always  seen  to 
be  satisfied,  thus  guaranteeing  the  feasibility  of  factorization.  Additionally, 
the  factorization  is  seen  to  be  nonunique.  Our  algorithm  for  computing  these 
factors,  however,  enjoys  two  remarkable  properties.  First,  it  encompasses  the 
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entire  family  of  possible  solutions.  From  this  point  of  view  it  may  be  remarked 
that  although  synthesis  in  cascade  type  structures  has  previously  been 
considered,  for  example,  in  (2),(14],{151  and  (3),  [13J  our  method,  namely  that 
of  factoring  the  corresponding  chain  matrix,  is  somewhat  more  general.  A  second 
important  property  is  that  the  factors  can  be  computed  essentially  by  solving  a 
highly  structured  set  of  linear  simultaneous  equations,  and  thus  can  be 
potentially  computed  in  a  fast  manner.  On  the  other  hand,  although  similar 
topological  structures  have  been  mentioned  in  (1)  in  the  context  of  wave 
digital  filters  obtainable  from  analog  prototypes,  we  are  unaware  of  any 
previous  work  on  discrete  domain  schemes  for  internally  passive  synthesis, 
which  yield  structures  as  a  result  of  repeatative  application  of  the 
decomposition  indicated  in  Figures  3.2  and  3.3.  Our  discussion,  even  in  the  1-D 
context,  thus,  yields  a  set  of  new  algorithms  for  structurally  passive 
synthesis  of  1-D  lossless  digital  filters  previously  not  discussed  in  the 
literature. 

Most  multidimensional  filtering  tasks  require  the  filter  to  have  certain 
symmetries  in  their  frequency  response  characteristics  [7].  However,  this 
demands  that  the  two-port  be  either  a  symmetric  or  a  (quasi)  antimetric 
two-pcrt.  Motivated  by  potential  applications  of  the  results  developed  in  the 
present  paper  to  the  design  of  multidimensional  structurally  passive  digital 
filters,  special  attention  to  the  synthesis  of  symmetric  and  (quasi) 
antimetric  lossless  two-ports  have  also  been  paid. 

In  section  3.2  a  precise  formulation  of  the  problem  along  with  some  notation 
and  terminology  are  introduced.  In  section  3.3  it  is  shown  that  the 
factorization  problem  so  introduced  is  essentially  algebraic  in  nature.  An 
elementary  step  towards  the  general  factorization  problem  is  also  taken  here. 
In  section  3.4  Some  properties  of  the  fundamental  equation  which,  in  fact,  is  a 
linear  version  of  the  algebraic  problem  and  is  central  to  our  study,  are 
examined.  Necessary  and  sufficient  conditions  for  factorability  and  an 
algorithm  for  obtaining  the  factors,  when  they  exist,  are  obatained  in  section 
3.5.  Section  3.6  discusses  how  the  results  so  obtained  yield  new  as  well  as 
known  internally  passive  1-D  digital  filter  structures.  In  section  3.7  remarks 
are  made  on  computational  considerations  of  the  algorithm  for  synthesis, 
examples  are  worked  out  to  demonstrate  the  need  for  factorability  of  three 
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different  kinds  of  matrices  associated  with  discrete  lossless  two-ports,  and 
the  special  cases  of  symmetric  and  (quasi)  antimetric  discrete  two-ports  are 
dealt  with.  Finally,  conclusions  are  drawn  in  section  3.8. 


3.2.  NOTATIONS,  TERMINOLOGY  AND  PROBLEM  FORMULATION: 


English  capital  letters  are  used  to  denote  polynomials  and  rational  functions 
«  (z^,  Z2»...,  z^).  The  notation  z-  is  used  to  denote  the 
.  ...z,k)  n  being  the  k-tuple  of  nonnegative  intergers  (n., 
^z...,  n^).  Also,  A  «  A  (z^  ,  Z2  ,  . . . ,  z^  )  where  *  denotes  complex 
conjugation.  The  notation  deg^A  will  be  taken  to  mean  the  partial  degree  of  the 
polynomial  A  in  the  variable  z^.  Occassionally  we  shall  also  use  the  notation  A 
to  denote  Az—  with  n-tn^,^, . .  .n^) ,  where  n^-deg^A.  Similarly  for  B  and  C  etc. 
The  notation  \z\  <  1  denotes  |z^|  <1  for  i«l  to  k.  Similar  notations  with  < 

replaced  by  <,  >,  >,  -  etc.  are  also  used. 

The  transfer  function  matrix  I  associated  with  a  k-D  lossless  two-port  can  be 
represented  [9]  as  in  (3.2.1).  Note  that  (3.2.1)  is  slightly  different  although 
an  equivalent  version  of  the  corresponding  representation  in  [9).  Consequently, 
the  chain  matrix  0  and  the  hybrid  matrix  T  can  also  be  represented  as  in 
(3.2.2)  and  (3.2.3)  respectively.  Also,  note  that  a  given  E  can  be  uniquely 
represented  as  in  (3.2.1)  by  requiring  that  A(0)-1.  The  representations 
(3.2.1),  (3.2.2)  and  (3.2.3)  can  be  regarded  as  canonic  in  this  sense. 


B  yCz- 

[a  rBZ- 

A  — yCz-^ 

E  -(1/A) 

”  1 

e-d/cr 

r-d/B) 

c  -yBz- 

B  yAZ- 

C  -yAz-  . 

~ 

— 

_  J 

(3.2.1,  2.2,  2.3) 


where  (i)  A  is  a  scattering  Schur  polynomial  [9)  (3.2.4a) 

(ii)  r  is  a  unimodular  constant,  i.e.,  |y|  -  1  (3.2.4b) 

( iii )  AA  -  BB  +  CC  (3.2.4c) 

(iv)  deg^^  B  <  n^ ,  degi  C  <  n^  for  all  i»l  to  k  (3.2.4d) 


Note  that  as  a  consequence  of  (3.2.4c)  and  (3.2.4d)  we  also  have: 

(v)  degi  A  <  n^  for  all  i»l  to  k.  (3.2.4e) 

For  the  purpose  of  a  unified  presentation  of  the  discussions  that  will  follow, 
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matrices  0,  r  and  E  associated  with  a  lossless  two-port  will  be  viewed  as  a 
matrix  ♦  associated  with  the  lossless  two-port  as  expressed  in  (3.2.5)  while 
properties  (Pi)  through  (P4)  hold  true. 


Fx  <tY2- 

*  -  (1/W)j 

j Y  pcXz- 


(3.2.5) 


Property  1  (Pi) :  X,  Y  and  W  are  polynomials;  a  ■  constant,  |c|~l,  |p|=+l. 
Property  2  (P2):  XX  -  pYY  *  WW. 

Property  _3  (P3) :  degi  W  <  n^,  deg^^  X  <  n^  and  deg^^  Y  <  n^  for  all  i-1  to  k. 
Property  4  (P4) :  X  is  scattering  Schur  when  p  *  1, 

whereas  W  is  scattering  Schur  when  p  ■  -1. 


A  (2  X  2)  matrix  such  as  the  one  in  (3.2.5)  will  be  said  to  be  a  generalized 
lossless  two-port  matrix  if  the  above  properties  (Pi)  through  (P4)  hold  true. 
Since  it  can  be  shown  that  *  satisfy  diag(l,  -p)  «  *.[diag(l,  -p)].*  and 
diag(l,  -p)  -  *.[diag(l,  -p)].*  >0  (i.e. , non-negative  definite)  in  |z|<l, 
where  *  denotes  the  Hermitian  transpose,  the  one-dimensional  counterpart  of  * 
thus,  falls  into  the  class  of  sigma-lossless  transfer  functions  studied  in 
[12].  Since  factorability  of  the  ^-matrices  form  the  major  topic  discussed  in 
the  present  paper,  out  results  can  also  be  viewed  as  a  multidimensional 
generalization  of  the  factorability  of  1-D  sigma-lossless  transfer  function 
elucidated  in  [12]. 


We  note  the  following  identification  of  the  parameters  of  the  matrix  *  in  terms 
of  the  parameters  of  the  chain  matrix  0,  the  hybrid  matrix  r  or  the  transfer 
function  matrix  Z. 


(i)  If  *  -  ©  then  W  «  C,  X  -  A,  Y  «  B,  a  =  y,  p  -  1 

(ii)  If  *  -  r  then  W  -  B,  X  -  A,  Y  -  C,  cr  =  -y,  p  -  1 

(iii)  If  *  -  E  then  W  -  A,  X  -  B,  Y  -  C,  a  -  y,  p  «  -1 

An  Advantage  of  the  above  formulation  is  that  the  problem  of  factoring  the 

chain  matrix,  or  the  hybrid  matrix  or  the  transfer  function  matrix  into  a 
product  of  two  non-trivial  matrices  of  identical  kind  can  be  conveniently 
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formulated  in  a  unified  fashion  as  a  single  problem  of  factoring  the 
generalized  lossless  two-port  matrix  *  as  *  -  where  *'  and  *"  are  valid 
generalized  lossless  two-port  matrices  represented  as  in  (3.2.6a)  and  (3.2.6b) 
with  conditions  analogous  to  (Pi)  through  (P4)  satisfied  for  #'  as  well  as  +". 


(1/W') 


X' 

a'Y'z^ 

X" 

~  n«“ 
o"Y"z- 

. 

t 

*”  -  (1/W"), 

Y' 

•ft 

pa'X'z-' 

Lr 

pff"X"z-!' 

(3.2.6a,  b) 


Consider  first  two  generalized  lossless  two-port  matrices  *' ,  ♦"  expressible 
respectively  in  terms  of  X',  Y',  W' ,  a'  and  X",  Y",  W"  ,  a"  which  satisfy 
properties  analogous  to  those  satisfied  by  X,  Y,  W,  o  in  *  as  in  (Pi)  to  (P4). 
Specifically,  one  obtains: 


(i)  X' ,  Y' ,  W'  and  X",  Y",  W"  are  polynomials; 

a'  and  a"  are  unimodular  constants.  (3.2.7) 

(ii)  X'X'  -  pY'Y'  -W'W',  X"X"  -  pY"Y"  -  W"W"  (3.2.8a,  b) 

(iii)  deg^^  W'  <  n?,  deg^  W"  <  nV  for  all  i«l  to  k  and  (3.2.9a) 

deg^  Y'  <  nJ,  deg^  X'  <  nf,  deg^  Y"  <  nV,  degi  X"  <  nV 

for  all  i-1  to  k  (3.2.9b) 

(iv)  X'  and  X"  are  scattering  Schur  when  p-1, 

whereas  W'  and  W"  are  scattering  Schur  when  p«-l.  (3.2.10) 


Then  the  following  fact  holds  true. 

Fact  3.2.1:  If  *'  and  #"  are  generalized  lossless  two-port  matrices  then  *  - 
♦'*"  is  also  a  generalized  lossless  two-port  matrix. 

Proof:  Given  *'  and  *"  define  X,  Y,  W  and  a  and  n-fn^,^, . .  .n^)  as  follows. 

W  -  W'W"  (3.2.11) 

X  -  X'X"  +  <t'Y'Y"z-\  Y  -  Y'X"  +  pv'X'Y" z-  (3.2.12a,  b) 

a  -  pv'ff"  (3.2.13) 

ni  -  nf  +  n?  for  each  i-1  to  k  (3.2.14) 
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It  then  follows  from  in  a  straightforward  manner  that  ♦  can  be  expressed 
in  terms  of  X,  Y,  W  and  a  as  in  (3.2.5).  Clearly,  (3.2,12),  (3.2.7)  and 
(3.2.9b)  show  that  *  satisfy  (Pi)  and  straightforward  algebraic  manipulations 
involving  (3.2.12),  (3.2.8),  (3.2.11)  and  p  -  +1  show  that  ♦  satisfy  (P2). 
Considering  the  degree  restriction  imposed  by  (3.2.9)  and  (3.2.14)  on  (3.2.11) 
and  (3.2.12)  likewise  shows  that  (P3)  is  satisfied.  Note  that  when  p  *  -1  i.e., 
W'  and  W"  are  scattering  Schur  polynomials  then  W  -  W'W"  is  clearly  scattering 
Schur  [9],  and  (P4)  is  thus  obviously  satisfied.  This  completes  the  proof  that 
when  p»- 1,  *-*'*"  as  expressed  in  (3.2.5)  via  (3.2.11)  through  (3.2.14)  is 
indeed  a  generalized  lossless  two-port  matrix.  On  the  other  hand,  when  p  -  1 
i.e.,  X'  and  X”  are  scattering  Schur,  X  as  in  (3.2.12a)  need  not  necessarily  be 
scattering  Schur,  but  can  be  shown  to  be  immittance  Schur  [9). 

However,  to  prove  that  I  is  a  generalized  lossless  two-port  matrix  we  consider 
the  rational  function  F  -  X/(X'X").  In  view  of  (3.2.12a)  we  also  have: 

F  -  X/(X'XM)  -  1  +  <r'(Y'/X'  )(Y"/X")z2'  (3.2.15) 

It  follows  from  (3.2.8)  with  p-1  that  on  |z|  ■  1  we  have  |Y'/X'|  <  1,  whenever 
X'  *  0  and  |Y"/X" |  <  1,  whenever  X"  *  0.  Thus,  Re  F  >  0  for  |z|  «  1,  whenever 
Vi'*  0  and  X"  *  0.  Next,  since  X'  and  X"  are  scattering  Schur,  by  invoking  Lemma 
3b  in  [9],  it  follows  that  Re  F  >  0  for  |z|  <  1  i.e.,  F  is  a  discrete  positive 
function.  Consequently,  the  numerator  polynomial  of  F,  in  irreducible  rational 
form,  is  a  immittance  Schur  polynomial  19).  Note  that  any  possible  factor 
common  to  X'X"  and  X  must  be  scattering  Schur,  because  X'  and  X"  are  so.  Thus, 
X  is  immittance  Schur  and  can,  therefore,  be  expressed  as  the  product  of  a 
scattering  Schur  factor  X^  and  a  reactance  schur  factor  D  19).  Let  be  any 
irreducible  (thus,  reactance  schur  19])  factor  of  D  and  note  that  there  exists 
a  sequentially  almost  complete  set  2  of  order  (k-1)  of  unimodular  complex 
numbers  (9)  such  that  «  X  •  0  for  any  z^  e  2.  Consequently,  in  view  of  (P2) 
we  conclude  that  Y  »  W  «  0  for  all  z^  e  2.  Since  ,  X,  Y,  W  have  a 
sequentially  almost  complete  (and  thus  sequentially  infinite  (9))  set  of  common 
zeros  of  order  (k-1)  and  D1  is  assumed  irreducible,  must  be  a  factor  of  X,  Y 
and  W  (6).  Since  is  any  irreducible  factor  of  D,  we  then  have  that  X  ■  X^D, 
Y  -  Y^D,  W  -  WjD,  where  X^,  Y^,  are  polynomials.  Since  D  is  reactance  Schur 
D  -  <jdD  for  some  unimodular  constant  ®D.  Clearly,  then  Xz-  -  (DppMX^-)  and 
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Yz—  -  (Dffp)  ( Y^z— ) ,  where  nHm^mj, . .  .m^)  with  n^-n^-  degiD  for  all  i-1  to  k. 
Thus,  after  cancelling  the  common  factor  D  from  the  numerator  and  denominator 
of  each  entry  of  (3.2.5),  ♦  can  be  written  as  in  (3.2.16),  where  -  aa^. 


* 


(1/V^) 


fflYl^2 


(3.2.16) 


Since  is  scattering  Schur,  Property  (P4)  is  satisfied  by  the  representation 
(3.2.16)  for  ♦.  It  can  be  further  shown  via  trivial  algebraic  manipulations 
that  Xj^,  Y^ ,  in  (3.2.16)  satisfy  properties  corresponding  to  (Pi)  to 
(P3)  because  X,  Y,  W,  a  has  been  shown  to  satisfy  the  same  properties. 


Note  that  in  the  case  p-1  if  D  is  a  nonconstant  polynomial  involving,  say,  , 
then  the  two-port  associated  with  *  is  degenerate  in  the  sense  that  in  (3.2.16) 
nn  <  n^-nf  +  n?.  The  main  problem  addressed  in  the  present  paper,  however,  is 
the  converse  problem  of  finding  a  non-degenerate  factorization  *-*'*"  of  a 
prescribed  generalized  lossless  two-port  matrix  *  into  two  factors  of  its  own 
kind.  More  specifically,  we  have  the  following  problem. 

MAIN  PROBLEM:  Given  a  generalized  lossless  two-port  matrix  ♦  as  in  (3.2.5),  two 
constants  a',  a"  such  that  |o'|  -  |c"|  -  1  and  o  -  o'  a",  and  the  polynomial 
factorization  W  -  W'W"  along  with  two  k-tuples  of  nonnegative  intergers  n'-(n£, 

n2'”‘nk^  an<*  n2*  sucl1  t^iat  ^e9i  w'  -  ni'  ^e^i  w"  ^  ni  arK*  ni  “ 

nj  +  n?  for  all  i-1  to  k,  we  seek  a  factorization  or  equivalently, 

find  polynomials  X',  Y',  X"  and  Y"  such  that  (3.2.12)  along  with  (3.2.8)  and 

(3.2.9)  hold.  Furthermore,  if  p-1  (or  p— 1)  then  we  require  X'  and  X"  (or  W' 

and  W")  to  be  scattering  Schur. 

It  proves  to  be  convenient  to  introduce  the  following  two  definitions: 

Definition  3.2.1:  The  pair  of  polynomial  two-tuples  {X',  Y'}  and  {X",  Y"}  is 
said  to  be  a  solution  to  the  algebraic  equation  if  (3.2.12)  along  with  (3.2.8) 
and  (3.2.9b)  are  satisfied. 
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Note  that  the  restrictions  that  the  polynomials  X'  and  X"  or  W'  and  W"  be 
scattering  Schur  polynomials  are  not  imposed  at  all  in  the  above  definition. 

Definition  3.2.2:  A  polynomial  triple  {X',  Y' ,  Y"}  is  said  to  satisfy  the 
fundamental  equation  if  (3.2.17)  along  with  (3.2.18)  holds  true. 

YX'  -  XY'  -  pc' Y"W'W'£— *  (3.2.17) 

deg.  X'  <  nf,  degj^  Y'  <  nr,  deg.  YH  <  n?  for  all  i-1  to  k  [ 3.2.18a,b,c) 

Note  that  (3.2.17)  is  obtained  by  adding  the  product  of  Y'  and  (3.2.12a)  to  the 
product  of  (3.2.12b)  and  (-X')  and  subsequently  by  using  (3.2.8a).  Obviously, 
then  any  solution  of  the  algebraic  equation  also  satisfies  the  fundamental 
equation.  However,  the  converse  statement  is  false,  consider  e.g.,  X'  «  0,  Y'  * 
0,  Y"  ■  0.  Note  further  that  the  algebraic  equations  (3.2.12)  along  with 
(3.2.8)  and  (3.2.9b)  constitute  a  highly  constrained  nonlinear  problem.  It  is 
shown  in  Section  3.4  that  due  to  the  inherent  structures  underlying  the  problem 
under  consideration,  solutions  to  this  nonlinear  equations  can  be  obtained  from 
a  certain  subclass  of  solutions  to  the  fundamental  equation,  which,  in 
contrast,  is  clearly  linear.  Thus,  solutions  to  the  algebraic  equation  can  be 
conveniently  characterized  in  terms  of  the  solutions  of  the  fundamental 
equation. 
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3.3.  SOLUTION  TO  THE  ALGEBRAIC  EQUATION: 

Clearly,  any  solution  to  the  problem  of  factorization  of  *  -  *'♦”  is  also  a 
solution  to  the  algebraic  equation.  The  converse  statement  is  obvious  if  *  is 
such  that  p—1  i.e.,  W  is  scattering  Schur.  The  validity  of  the  converse 
statement  when  p-1  i.e.,  when  X  associated  with  ♦  is  scattering  Schur  is 

nontrivial  but  can  be  proved  as  follows  in  Theorem  3.3.1.  Thus,  the  problem  of 
factoring  *  reduces  to  that  of  solving  a  purely  algebraic  problem  namely  that 
of  finding  a  solution  to  the  algebaic  equations. 

Theorem  3.3.1:  Let  the  pair  of  polynomial  two  tuples  {X',  Y'}  and  {X”,  Y"} 
constitute  a  solution  to  the  algebraic  equation.  If  p  -  1  (or  p  «  -1)  then  the 
polynomials  X'  and  X"  (correspondingly,  W'  and  W")  are  scattering  Schur.  Thus, 
any  solution  to  the  algebraic  equations  is  a  solution  to  the  problem  of 
factoring  *  - 

Proof:  The  case  when  p  -  -1  trivially  follows  from  scattering  Schur  property 

of  W  -  W'W”.  When  p  «  1  i.e.,  X  is  scattering  Schur,  consider  the  rational 
function  defined  as: 

F  -  (X'X")/X  (3.3.1) 

Furthermore,  by  adding  the  product  of  (3.2.12a)  and  (X')  to  the  product  of 

(3.2.12b)  and  (-pY')  and  subsequently  by  using  (3.2.8a),  one  obtains 

X"  ■  (XX'  -  pYY' )/(W'W' ) .  Substituting  the  last  expression  in  (3.3.1) 
straightforward  manipulation  yields  the  following: 

F  -  (X'X'/W'W' ) { 1  -  (pYY'/XX')]  (3.3.2) 

It  follows  respectively  that  from  (P2)  and  (3.2.8a)  that  on  |z|-l  we  have  |Y/X| 
<  1,  whenever  X  *  0  and  |Y'/X'|  <  1,  whenever  X'  *  0.  An  examination  of  (3.3.2) 

yields  that  ReF  >  0  for  |z|  ■  1,  wherever  F  is  well  defined.  Thus,  from  Lemma 

3b  in  (9]  it  follows  that  F  is  a  discrete  positive  function.  Consequently,  the 
numerator  polynomial  of  F,  in  irreducible  form,  is  an  immittance  Schur 
polynomial.  Note  that  any  possible  factor  conmon  to  X'X"  and  X  must  be 
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scattering  Schur,  because  X  is  so  [9],  and  thus,  X'X"  is  widest  sense  Schur 
(more  specifically,  immittance  Schur).  Next,  if  for  some  Zq  on  the 
distinguished  boundary  \z\  -  1,  we  have  x"»  0,  then  from  (3.2.8b)  it  follows 
that  Y"«  0,  which  in  turn  due  to  (3.2.12a)  imply  that  X  -  0.  Consequently,  if 
X"-  0  for  all  Zq  c  2,  where  2  is  sequentially  infinite  set  (6]  of  order  (k-1) 
of  unimodular  complex  numbers  then  X  -  0  for  all  z^  e  2.  However,  this  is 
impossible  if  X  is  scattering  Schur  (9).  Therefore,  X"  cannot  have  a 
sequentially  infinite  set  of  zeros  of  order  (k-1)  on  the  distinguished 
boundary.  The  scattering  Schur  property  of  X"  is  thus  established  in  view  of 
Theorem  9  in  [9].  Similar  arguments  hold  for  X'. 

Proposition  3.3.2:  Any  generalized  lossless  two-port  matrix  #  can  be  decomposed 
as  ♦  «  *l*o*2'  where  *i»  *o'  *2  are  generalized  lossless  two-port 
matrices  such  that  ,  *2  are  diagonal  and  the  polynomials  XXz—  and  WWz- 
associated  with  are  coprime. 

Proof:  Let  H  ■  gcd(W,  X),  where  W-HWj,  X^HX^.  Since  W  or  X  is  scattering  Schur 
H  is  also  so.  It  then  follows  from  property  (P2)  of  (3.2.5)  that  H  is  a  factor 
of  YY.  Let  H-H'H",  where  H',  H"  are  factors  of  Y,  Y  respectively.  Since  H" 
divides  Y,  H"  divides  Y.  Due  to  the  scattering  Schur  property  of  H"  inherited 
from  H,  H"  and  H"  are  coprime  and  thus  Y«(H"H")Y^  for  some  polynomial  Y^ .  A 
direct  substitution  of  the  last  equation  along  with  W-HW^ ,  X-HX^  in  XX  -pYY-ww 
yields  (3.3.3a)  in  the  following. 

*1*1  ~  pYlYl  “  WA'  ¥c  -  pYcYc  "  wc«c  ( 3.3.3a,b) 

Next,  let  F  -  gcd(X^,  W^),  where  X^  «FXc,  -  FWC  and  the  monomial  factors  of 
of  maximal  degree  in  XQ  and  X^  are  identical  (note  that  this  uniquely  defines 
Xc  upto  a  constant  multiplier).  Clearly,  F,  being  a  factor  of  X^^  cannot  have  a 
monomial  factor  and  thus,  from  (3.3.3a),  F  must  divide  Y.Y,.  Let  F-F'F",  where 
F'  divides  Y^  and  F"  divides  Y^.  The  last  requirement  implies  F"  divides  Y^.  If 
p—  1  then  F  is  scattering  Schur  since  W,  is  also  so.  On  the  other  hand,  if  p»l 
since  X^-FXQ,  then  F  is  scattering  Schur,  because  X-HX^  is  also  scattering 
Schur.  Thus,  either  both  F'  and  F"  or  both  F'  and  F"  are  scattering  Schur. 

A  A 

Consequently,  F'  and  F"  are  coprime,  and  Y.  -  (F'F")Y,.  for  some  polynomial  Y  . 
Next,  if  we  define  matrices  ♦-  -  Diag(H"F',  H"F'),  ♦  -  Diag(H'F",  H'F")  and  ♦ 


asj  l*clll"x</l,c# 

« •  •  •  ) »  ® 

*“*r*c*f  anc*  eac^  t*'e  “^trices  so  defined  is  a  generalized  lossless  two-port 
matrix  with  same  p.  In  particular,  (3.3.3b)  holds  true,  where  is  at  least  as 
large  as  deg^,  deg^,  deg^Wc  for  each  i-1  to  k. 

Clearly,  Wc  is  coprime  with  Xc.  On  the  otherhand,  since  it  follows  from  the 
definition  of  Xc  that  X^-FXc,  Wc  and  Xc  cannot  have  a  common  factor,  because 
otherwise,  due  to  W^-FW,,  V?1  and  X^  would  not  be  coprime.  Consequently,  XcXc 
and  W  W  are  coprime.  If  m.«deg.W„  for  all  i-1  to  k  then  the  last  conclusion 
also  implies  that  the  polynomials  WcWcz~  and  X^^  are  coprime,  the  proof  of 
the  present  theorem  is  complete.  Otherwise,  diagonal  lossless  two-port  matrices 
from  left  and/or  right  of  ♦  needs  to  be  further  extracted  to  satisfy  the 
coprimeness  requirement. 

For  this  purpose,  note  that,  due  to  (3.3.3b),  any  monomial  factor  of  X  present 
~  - n  **'  _  c 

in  WcWcz—  must  also  be  a  factor  of  Y^z-.  Let  S-SjS2  be  such  a  factor  of 

maximal  (total)  degree,  where  the  monomial  Sj  divides  Yc,  and  the  monomial  S2 

divides  Ycz~.  Consider  polynomials  X^-X^S,  Y^-Y^ys^,  Wc,-Wc  and  the  integer 

k-tuple  m'»(m£,m£, , .  .m£)  with  mi  -  nu-deg^S.  Then  clearly  Xc,Xc,  -  pYc,Yc,  - 

Wc,Wc,  with  degiXc,,  degiYc,  and  deg^Wc,  upper  bounded  by  m!  for  all  i«l  to  k 

holds  true.  Thus,  any  monomial  factor  of  X  ,z-  present  in  W  ,W  ,z-  must  also 

*■*  c  c  c 

be  a  factor  of  Yc,Yc,z-  .  Let  T»TjT2  be  such  a  factor  of  maximal  (total) 
degree,  where  the  monomial  divides  Yc,,  and  the  monomial  T2  divides  Yc,zi— . 
Next,  consider  XQ-Xc , ,  Yq-Yc,/T  and  the  integer  k-tuple  n^-tn^,  nQ2 , . . . no)( ) , 
where  nQj“  mj-deg^T.  By  letting  WQ«Wc  it  is  then  routinely  verified  that  *c  - 
Diag(S2,  Ti> ,*Q.Diag(S^,  T2),  where  *o  is  a  generalized  lossless  twc^port 
matrix  described  „by  X  .  Y  .  W  .  n  and  o  as  in  (3.2.5)  such  that  w  w  z~°  is 
coprime  with  XqX^-0 .  The  proof  of  the  theorem  is  then  completed  by  setting: 

-  Diag(S2,  T^.*^  and  *2  -  Diag(S1,  T2).*  . 


m-  ( m^  t  ^2 


l«c>21  'Ve'  (,c>12  -(Yci")/WC  I‘c)22-|xc5.-|/Mc'  uhere 
^-n^-deg^(HF)  it  follows  in  a  straightforward  manner  that 
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3.4.  PROPERTIES  OF  THE  FUNDAMENTAL  EQUATION. 


In  this  section  certain  properties  of  the  fundamental  equation  (3.2.17)  crucial 
to  the  development  are  studied  under  the  assumption  that  the  prescribed 
generalized  lossless  two- port  matrix  *  be  such  that  the  polynomials  XXz—  and 
WWz—  are  relatively  prime.  As  shown  in  Proposition  3.3.2  no  loss  of  generality 
is  incurred  due  to  this  assumption. 

Lemma  3.4.1:  If  the  polynomial  triple  {X',  Y',  Y")  is  a  solution  to  the 
fundamental  equation  then  there  exist  a  polynomial  X"  given  by  (3.4.1)  such 
that  the  polynomial  triple  {pY'z-  ,  X'z-  ,  -po'  X")  is  also  a  solution  to  the 
fundamental  equation. 


where 


X"  -  P/( Xz— )  -  -^(W’W'z-' ) 


(3.4.1a) 


(X'W"W"z—  +  a'  *Y"z— "y)  ,  Q 


Y(pY'Z-' )  -  X(X'z-' ) 


(3.4.1b,  C) 


Proof:  One  obtains  (3.4.2)  by  adding  the  product  of  tilde  of  (3.2.17)  and  pY  to 
the  product  of  (P2)  and  X'  and  subsequently  by  using  (3.2.11),  (3.2.14)  and 
trivial  manipulations. 


(Xz-)Q  -  -(W'W'z-'  )P 


(3.4.2) 


Due  to  the  upper  bounds  on  the  degrees  of  X',  Y' ,  Y"  imposed  by  (3.2.18)  and 
deg^^  W"<  n?,  it  follows  from  (3.4.1b,  c)  that  P  and  Q  are  polynomials.  The  fact 
that  X"  in  (3.4.1a)  is  a  polynomial  then  follows  from  (3.4.2)  in  view  of 
relative  primeness  of  Xz-  with  W'W'£-  .  Since  Q  -  Y(pY'£-  )  -  X(X'£-  )  - 
po' (-po' *X"  )W'W'£— * ,  the  triple  {pY'z-* ,  X'z^' ,  -o'*X")  satisfies  the 

fundamental  equation. 


**  rtf  *’  «  ^ 

The  fact  that  deg^(pY'£-  )  and  deg^(X'z-  ),  for  each  i-1  to  k,  is  upper  bounded 
by  n;  is  obvious.  In  order  to  prove  that  deg^X"  <  n7  we  first  note  that  it 
follows  from  (3.4.1b,  c)  and  the  upper  bounds  on  the  degrees  of  X,  Y,  X',  Y' , 
X"  and  Y"  that  for  all  i-1  to  k  we  have: 
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deg^P  <  n^n?,  deg^  <  (3.4.3a,  b) 

It  is  then  necessary  to  distinguish  between  the  following  two  cases. 

(i)  If  p-1  i.e.,  X  is  scattering  Schur  then  deg^(Xz-)  -  n^  for  all  i-1  to  k. 
It  then  follows  from  the  first  equality  in  (3.4.1a)  and  (3.4.3a)  that  degiX"  < 
n?  for  all  i*l  to  k. 

(ii)  If  p— 1  i.e.,  W,  and  thus  W' ,  is  scattering  Schur  we  consider  two  sets  of 
indices  1^,  Ij  such  that  i  e  1^  if  deg^'  ■  nf,  whereas  i  e  Ij  if  deg^W'  <  nr. 
If  i  e  l1  then  due  to  scattering  Schur  property  of  W'  we  have  deg^(W'W'^—  )  - 
2nr.  The  desired  result  then  follows  from  (3.4.3b)  and  second  equality  in 
(3.4.1a).  On  the  other  hand,  if  i  e  ^  then  from  (3.2.9a),  (3.2.11)  and 
(3.2.14)  it  follows  that  Wz-  must  have  a  factor  z^,  and  thus,  X  does  not  have  a 
factor  z^  because  Wz-  and  X  are  assumed  to  be  coprime.  Consequently,  deg^Xz— 
-n^  The  result  then  follows  from  (3.4.3a)  and  first  equality  in  (3.4.1a). 


Lemma  3.4.2:  If  (X£,  Y£,  Yj}  and  {X£,  Y£,  Y^}  are  two  polynomials  triples 
satisfying  the  fundamental  equation  then  the  identity  (3.4.4)  holds  and  is 
equal  to  a  constant. 

N  -  pa'(YpU,  -  X'Y£)/X  »  (X'Y£  -  X^Y' >/(W'W' z-  )  (3.4.4) 

Proof:  One  obtains  an  equivalent  form  of  (3.4.4)  by  adding  the  product  of  the 
fundamental  equation  for  {X£,  Y£,  Y£}  and  X£  to  the  product  of  the  fundamental 
equation  for  (X^Y^Yfp  and  (-X£).  Since  X  is  assumed  coprime  with  W'W'z—  and 
(X'Y£  -  X^Y£)  is  a  polynomial,  it  follows  that  X  must  divide  (YpC^  -  X^Yp. 
Thus,  N  -  pa'  ( YjX^-XjYlp/X  in  (3.4.4)  is  a  polynomial.  To  prove  that  N  is  a 
constant,  note  that  the  following  inequalities  hold  true  for  all  i*=l  to  k. 

degi ( Y£X£  -  X£Y!|)  <  n.,  deg.(X'Y£  -  X£Y')  <  2nJ  (3.4.5a,  b) 

Consider  two  sets  of  indices  1^,  ^  such  that  i  e  1^  if  deg^X  »  n^,  whereas  i  e 
I2  if  deg^X  <  m.  If  i  e  I ^  then  from  (3.4.5a)  and  first  equality  (3.4.4)  it 
follows  that  N  does  not  involve  z. .  If  i  el-  then  Xz-  has  the  a  factor  z. . 

X  4  X 
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Consequently,  due  to  the  assumed  relative  primeness  of  XXz-  and  WWz-  it  follows 
from  (3.2.9a),  (3.2.11)  and  (3.2.14)  that  neither  W'  nor  W'z-  may  have  the 

"  n r 

factor  which  in  turn  respectively  imply  that  deg^W'z-  )»n!  and  deg^W'»n?. 
Therefore,  due  to  (3.4.5b)  and  the  second  equality  (3.4.3a)  N  may  not  involve 
Thus,  N  ■  constant. 

t -omnia  4.3:  If  {X',  Y' ,  Y")  is  a  polynomial  triple  satisfying  the  fundamental 
equation  then  the  expression  given  in  (3.4.6)  is  equal  to  a  real  constant. 

K  -  (X'X'  -  pY'Y'  )/(W'W'  )  -  (ff'Y'z-V  +  X'X'*)/X  (3.4.6) 

Proof:  Consider  in  view  of  Lemma  3.4.1  two  solutions  X£  «  X' ,  Y£  ■  Y' ,  Y£  ■  Y" 
and  X^  -  pY'z— * ,  Y^  -  X'z-  ,  Yi|  -  -po'*X"  to  the  fundamental  equation.  It  then 
immediately  follows  from  Lemma  3.4.2  that  K  in  (3.4.6)  is  a  constant.  Since  for 
|z|-l  we  have  X'X'«|X'|2,  Y'Y'-|Y'|2  and  W'W'«|W' |2,  K  is  a  real  constant. 

Lemma  3.4.4:  If  the  polynomial  triple  {X',  Y' ,  Y")  is  a  solution  to  the 

fundamental  equation  then  there  exists  an  X"  as  given  by  Lemma  3.4.1  such  that 
(oX'  +  (SpY's-  ,  aY'  +  px'z-  ,  aY"  -  0pc'*X"}  is  also  a  solution  to  the 

fundamental  equation,  where  a  and  0  are  arbitrary  complex  numbers. 

Proof:  Follows  clearly  from  Lemma  3.4.1  and  the  fact  that  the  fundamental 
equation  is  linear. 
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3.5.  FACTORIZATION  OF  *: 


A  solution  {X',  Y' ,  Y")  to  the  fundamental  equation  will  be  called  nonsingular 
if  X'X'  *  pY'Y'. 

Theorem  3.5.1:  The  problem  of  factorization  of  *  admits  a  solution  if  and  only 
if  there  .  exists  a  nonsingular  solution  {Xr ,  Y',  Y"}  to  the  fundamental 
equation. 

Proof:  Necessity  obviously  follows  from  (3.2.8a)  and  that  W^W'WVO.  If  {X£,  Y£, 
Y£)  is  a  nonsingular  solution  to  the  fundamental  equation  then  due  to  Lemma 
3.4.4,  X'  -  oX£  +  0pY'z-\  Y'  -  «Y'  +  0X'z-' ,  YM  -  «Y£  -  0p«^*X£  is  a  solution 
to  the  fundamental  equation.  Straightforward  manipulation  then  yields: 

(X'X'  -  pY'Y' )/(W'W' )  -  ( |a| 2  -  p(e| 2)K1  (3.5.1) 

K1  *  (X1X1  ~  P*{Yp/(W'W')  (3.5.2) 

Since  due  to  Lemma  3.4.3  and  nonsingularity  of  (X.',  Y.' ,  YIM ,  K.  is  a  nonzero 
constant,  if  a  and  0  are  chosen  to  satisfy  ( |a|  -  p|0|  )  -  K.  we  have 
(X'X'-pY'Y' )  -  W'W'.  Furthermore,  there  exists  X"  such  that  (pY'z-  ,  X'z-', 
- pa £  X"),  by  the  virtue  of  Lemma  3.4.1,  satisfies  the  fundamental  equation. 

We  next  show  that  X',  X",  Y'  and  Y"  so  obtained  constitute  a  solution  to  the 
algebraic  equation.  Equation  (3.2.12b)  is  obtained  by  adding  the  product  of 
(3.2.17)  and  X'z—  to  the  product  of  second  equality  of  (3.4.1a)  via  (3.4.1c) 
and  (-Y' )  and  subsequently  by  using  (3.2.8a).  Likewise,  (3.2.12a)  is  obtained 
by  adding  the  product  of  (3.2.17)  and  pY'£—  to  the  product  of  second  equality 
of  (3.4.1a)  via  (3.4.1c)  and  (-X')  and  subsequently  by  using  (3.2.8a). 
Finally,  we  obtain  (3.2.8b)  by  substitituting  (3.2.12a)  and  (3.2.12b)  in  (P2) 
and  then  using  (3.2.11)  and  (3.2.8a).  Thus,  the  pair  of  two-tuples  {Xr,  Y'}  and 
{X",  Y"}  satisfies  the  algebraic  equation  and  via  Theorem  3.1  is  a  solution  to 
the  problems  of  factorization  of  *. 

Two  polynomials  triples  (X£,  Y£,  Y£}  and  {X£,  Y£,  Y^}  each  satisfying  the 
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fundamental  equation  will  be  said  to  be  linearly  dependent  if  there  exists 

constants  a  and  0  not  simultaneously  zero  such  that  aX£  +  0X£  m  aY£  +  0Y£  *  aY£ 

+  0Y£  *  0. 

Theorem  3.5.2:  The  problem  of  factorization  *  admits  a  solution  if  and  only  if 
there  exists  two  linearly  independent  polynomial  triples  {xr,  Yr,  YV},  i-1,2 
each  of  which  satisfy  the  fundamental  equation. 

Proof:  Necessity:  Let  the  polynomials  X' ,  Y' ,  X"  and  Y"  constitute  a  solution 

to  the  factorization  problem.  Clearly,  {X',  Y',  Y"}  is  a  solution  to  the 

~  n'  ""  n'  * 

fundamental  equation.  Due  to  Lemma  3.4.1,  therefore,  {pY'z—  ,  X'z—  ,  -pa'  X"} 

is  also  a  solution  to  the  fundamental  equation,  we  claim  that  these  two 

solutions  are  linearly  independent,  because  otherwise  there  would  exist 

~  n 9 

constants  (3. ,  02  not  simultaneously  zero,  such  that  d^X'  +  02pY'z-  *  0jY'  + 

02X'z-  *  0.  Thus,  (X'X'-pY'Y' )  *  0,  which  in  view  of  (3.2.8a),  would  imply 

that  W'  »  0,  i.e.,  due  to  (3.2.11)  that  W  *  0,  which  is  inpossible. 

Sufficiency:  If  one  of  the  solutions  {xr,  Yr,  Y?},  i-1,2  is  nonsingular  then 
sufficiency  follows  from  Theorem  3.5.1.  If  both  solutions  are  singular  then  the 
triple  {X'f  Y',  Y"}  obtained  as:  X'  -  pX£  +  qX£,  Y'  -  pY'  +  qY£,  Y”  -  pY£  + 
qY£,  where  p  and  q  are  complex  numbers,  satisfy  the  fundamental  equation. 
Straightforward  algebraic  manipulation  via  the  singularity  of  the  triples  {xr, 
Yr,  Y7},  i-1,2  then  yields: 

(X'X'-pY'Y' )/(W'W' )  -  pq*L  +  p*qL,  L  -  (X£X£  -  pY'Y£)/(W'W' )  (3.5.3a,b) 

**  n  >  **  it 

Since  due  to  Lemma  3.4.1  {pY£z-  *  X£z-  '  ”p0  *2^  a^s0  a  so^ut^on*  b Y 

invoking  Lemma  3.4.2  on  the  triples  {X£,  Y£,  Yj}  and  {pYJ£—  >  X^z-  •  -P^Xjj}  it 
follows  that  L  in  (3.5.3b)  is  a  constant  i.e.,  L-L-L  .  Thus,  the  right  hand 
side  of  (3.5.3a)  is  2Re(pq*L),  which,  if  L*0,  can  be  made  equal  to  1  by  proper 
choice  of  p  and  q.  With  p,  q  so  chosen  {X',  Y' ,  Y"}  would  thus  be  a  nonsingular 
solution  to  the  fundamental  equation,  and  by  invoking  Theorem  3.5.1,  it  then 
follows  that  a  solution  to  the  problem  of  factorization  of  *  exists. 

The  proof  of  the  present  theorem  is  next  completed  by  showing  that  L*0.  For 
this,  consider  the  following  cases. 
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(i)  p— 1:  Assume  for  contradiction  that  L  -  0,  which  due  to  (3.5.3a)  implies 
that  X'X'  +  Y'Y'  -  0  and  thus,  |X'|2  +  I Y' 1 2  -  0  for  |z|«l.  Consequently,  X'  a 
0,  Y'  *  0,  and  via  (3.2.17)  Y"  ■  0,  which  contradicts  linear  independence  of 
{Xf,  Y!,  Y?},  i-1,2. 

(ii)  p  ■  1:  Assume  for  contradiction  that  L  -  0,  which  due  to  (3.5.3b)  implies 
X'X'  -  Y£Y£.  Since  (X',  Y',  YJJ  is  singular,  we  have  X£X'  -  Y'Y'.  The  last  two 
equations  together  imply  X^/X£  -  Yi>/Y£  *  H2//Hl'  where  Hi  H2  are  C0Prime 
polynomials.  Clearly,  there  exists  polynomials  XQ',  Y  '  such  that 

X1  "  HlXo'  Y1  "  HlYo'  X2  "  h2Xq'  Y2  *■  H2y;  (3.5.4a,b,c,d) 

Considering  the  fundamental  equation  for  the  triple  {X£,  Y£,  Y£)  and  {X£,  Y^, 
Y:;),  we  obtain  (3.5.5a,b)  where  Y"  is  defined  via  (3.5.5c). 

Y1  *  HlYo'  Y2  ■  H2Yo'  (YXi  ■  ^  ■  PO'Z-  W'W'Y^  (3.5.5a,b,c) 

Clearly,  Yo"  is  a  polynomial,  since  otherwise  its  least  denominator  would 
divide  both  and  H2,  i.e.,  and  H2  would  not  be  coprime. 

Furthermore,  it  follows  from  (3.5.4a),  (3.5.4b)  and  (3.5.5a)  that  the  degrees 

of  the  polynomials  in  the  triple  {X^,  Y£,  Y£}  cannot  exceed  the  degrees  of  the 

corresponding  polynomials  in  the  triple  {X£,  Y£ ,  Y£).  Thus,  in  view  of  (3.5.5c) 

{X',  Y',  Y"},  and  consequently  due  to  Lemma  3.4.1,  {pY'z-  ,  X'z-  ,  -p<j*'X"}  is 

a  solution  to  the  fundamental  equation  for  some  X".  This  last  mentioned 

0  -  n'  -  n' 

equation  along  with  (3.5.4)  and  the  fundamental  equations  for  {pYJz:—  ,  X£z—  , 
-pa*XJ}  (cf.  Lemma  3.4.1)  and  {pY£z-  ,  X£z-  ,  -pcr'*X2}  yield: 

X1  “  HlXo'  X2  "  H2Xo  ( 3 . 5 . 6a , b ) 

Next,  if  we  define  FQ  -  X^X^/X,  then  by  eliminating  X£  between  Fq  and  the 
fundamental  equation  for  the  triple  {pY'z-  ,  X'z-  ,  -pc*X"}  one  obtains 

O—  O—  o 

(3.5.7). 

Fq  -  (X'XyW'W')[l  -(  pYY^/X^X) )  (3.5.7) 
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From  Property  (P2)  of  *  in  (3.2.5)  we  have  | Y/X }  <  1  on  |z|=*l,  whenever  X  t  0. 
Since  {X£,  Y£,  Y£}  has  been  assumed  to  be  a  singular  solution  to  the 
fundamental  equation,  we  have  X£X£  *  Y^Y^,  ^hus  v*a  (3.5.4a,  b),  X^X^  =  Y^Y^, 
implying  that  -  1  on  |z|-l,  whenever  Y^  t  0.  Then  (3.5.7)  yields  that 
ReFo  >  0  for  |z|-l,  wherever  FQ  is  well  defined.  Using  an  argument  analogous  to 
that  used  in  the  proof  of  Theorem  3.3.1  it  then  follows  (via  discrete  positive 
nature  of  F  )  that  X^,  X£  are  widest  sense  Schur  polynomials,  and  thus,  cannot 
contain  z ^  as  a  factor  for  any  i*l  to  k.  Since  in  (3.5.6)  X£  and  x£  are 
polynomials,  it  then  follows  that  and  Hj  are  constants.  This  latter 
conclusion  violates,  due  to  (3.5.4)  and  (3.5.5a,  b),  the  linear  independence  of 
(X£ ,  YJ,  Y£}  and  {X^,  Y£,  Y"}. 

The  above  result  can,  in  fact,  be  further  sharpened  as  follows.  Note  that  a 
corresponding  strong  result  for  the  continuous  case,  although  true,  was  not 
given  in  [ 5) . 

Theorem  3.5.3:  The  problem  of  factorization  of  *  admits  a  solution  if  and  only 
if  there  exists  exactly  two  linearly  independent  polynomial  triples  {X?,  Y? , 
YV},  i-1,2  each  of  which  satisfy  the  fundamental  equation. 

Proof:  Due  to  Theorem  3.5.2,  it  is  only  needed  to  show  that  one  can  have  at 
most  two  linearly  independent  solutions.  Assume  for  contradiction  that  {XJ,  Y£, 
YV},  i»l  to  3  be  three  linearly  independent  solution  to  the  fundamental 
equation  and  X'  *  cx^XJ  +  +  o^,  Y*  -  ot^  +  +  a^,  Y"  -  O^YJ  +  o^Y^ 
+  ot^Yj.  Then  {X',  Y' ,  Y"}  and,  due  to  lemma  3.4.1,  {py'z-  ,  x'z-  ,  -pa'*X”)  is 
also  a  solution  to  the  fundamental  equation.  From  Lemma  3.4.3,  K  as  defined  in 
(3.5.8)  is  a  constant . 

K  -  (X'X'-pY'Y'  )/(W'W'  )  *  (a'Y'z-V  +  X'X")/X  (3.5.8) 
Consider  next  the  following  cases: 

(i)  p-1:  Clearly,  we  may  choose  0^,  i*l  to  3  such  that  X' (£)  -  Y"(0)  *  0.  Since 
X(0)*0  due  to  the  scattering  Schur  property  of  X,  it  follows  from  (3.5.8)  that 
K-0  and  thus,  again  from  (3.5.8)  that  X'X'  *  pY'Y' .  Next,  define  F  -  (X'X")/X. 
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From  the  second  equality  of  (3.4.1a)  and  (3.4.1c),  one  obtains  X"  -  (XX'  - 
pYY' )/{W'W' ) .  Eliminating  X"  from  the  last  two  equations  it  follows  that  F  - 
(X'X'/W'W')(1  -  ( YY'/XX' ) ) .  From  Property  (P2)  of  *  and  X'X'  -  pY'  Y' ,  one 
respectively  obtains  for  \z\  -  1  that  |Y/X|  <  1,  whenever  X*0  and  |Y'/X'|  -  1, 
whenever  X'#0.  Using  an  argument  similar  to  that  used  in  proving  Theorem  3.3.1, 
it  is  then  possible  to  show  that  F  is  a  discrete  positive  function,  and 
finally,  X'X"  is  widest  sense  Schur,  which  is  in  contradiction  with  our 
construction  that  X'(0)«0. 

(ii)  p— 1:  We  claim  that  there  exists  Zq  on  |z|  -  1,  such  that  W'W'  *  0  [9]  and 
c^,  i-1  to  3  cam  be  chosen  such  that  X'(z^)  -  Y'(Zq)  *  0.  It  then  follows  from 
(3.5.8)  that  for  z  -  z.,  K  «  0,  i.e.,  |X' |2+|Y' |2  *0  implying  X'  *  0,  Y'  *  0, 
and  via  (3.2.17),  Y"  *  0,  which  contradicts  linear  independence  of  {X?,  Y£, 
Y7},  i-1  to  3. 

To  substamtiate  the  claim  we  show  that  there  exists  z-z^  on  |z|  -  1  with  W'W'  * 
0  such  that  XjY^  -  X^Y£  *  0  and  thus,  it  is  possible  for  any  nonzero  to 
solve  the  following  linear  siraul taneus  equations  for  and  c^. 

X'^)  -  OjX'fz^  +  o^fz^  +  «3X3<zo)  *  0  (3.5.9a) 

Y'  (Zq)  -  o^z^  +  c*^^)  +  c^Y't^)  -  0  (3.5.9b) 

For  this,  consider  a  solution  X'  -  X',  Y.'  -  Y',  X."  -  X",  YV  -  Y"  to  the 
factorization  problem.  Then  from  the  necessity  part  of  proof  of  Theorem  3.5.2 
it  follows  that  X£  -  pY^z-  ,  Y^  -  X^z-  •  Y£  *  -p<x'*X£  is  also  a  solution  to  the 
fundamental  equation  and  that  (Xf,  Yf,  Y? } ,  i-1, 2  are  linearly  independent. 
Next,  for  some  z^  on  |z|  -  1  if  X'Y^-X^Y'  -  0  then  |X^|2  +  |Y^|2  -  0  i.e., 
X^  -  0,  Y£  -  0,  amd  thus,  due  to  the  equation  corresponding  to  (3.2.8a) 
satisfied  by  X^,  Y^  we  have  W'  «  W'  -  0.  Since  W'  is  scattering  Schur,  there 
exists  z^  on  |z|-l  with  W'*0  (9]  amd,  consequently,  with  X£Y^  -  X£Y£  t  0. 

The  fundamental  equation  (3.2.17),  when  considered  as  a  set  of  linear 
simultaneous  equations  involving  the  coefficients  of  the  polynomials  X' ,  Y' ,  Y" 
along  with  the  upper  bounds  on  their  degrees,  turns  out  to  be  overdetermined 
in  general  (except  when  k-1).  More  explicitly,  we  note  that  the  unknown 


polynomials  X',  Y'  and  Y"  contain  a  total  of  u  unknown  coefficients,  whereas 
the  total  number  of  linear  simultaneous  equations  can  easily  be  found  to  be 
equal  to  e,  u  and  e  being  as  given  in  (3.5.10a,b)  below. 

k  k  k 

u  -  n  (nV+l)  +  2  n  (n;+l),  e  -  n  (nr  +  2n{  +  1)  (3.5.10a,b) 

i-1  1  i-1  1  i-1  1  1 

Since  for  k>l  we  have  e>u  in  a  generic  situation  a  solution  to  the  problem  of 
factoring  #  into  two  matrices  of  identical  kind  may  not  exist. 

Delay  free  loop:  In  order  for  the  digital  network  synthesized  via  the 
factorization  of  #  to  be  'computable'  it  may  not  contain  delay  free  loops 
arising  from  interconnection  of  two  sections.  It  is  known  [1]  that  this  problem 
can  always  be  circumvented,  at  least  in  the  one-dimensional  case,  by 
incorporating  digital  equivalents  of  unit  elements.  The  structures  resulting 
from  factorizations  0  -  0'0",  r  «  PP,  E  -  E'E”  are  shown  in  Figures  3.1,  3.2 
and  3.3  respectively.  An  examination  of  directions  of  signal  flows  in  Figure 
3.3  shows  that  the  topological  structure  arising  from  the  factorization  of  E 
as  E  «  E'E"  cannot  contain  any  delay  free  loop  at  the  junction  of  the  two-ports 
E'  and  E".  On  the  otherhand,  Figures  3.1  and  3.2  clearly  show  that  the 
topological  structures  arising  from  the  factorization  of  ©  as  0  »  0'©"  and  r  as 
T  -  FT"  may  contain  delay  free  loops  unless  special  attention  is  paid  to  this 
issue  (note  that  both  of  these  cases  correspond  to  the  choice  p-1).  However,  as 
shown  in  the  following,  delay  free  loops  at  the  junction  of  the  two-ports 
associated  with  #'  and  *"  may  always  be  avoided  by  extracting  an  appropriate 
constant  (generalized)  lossless  two-port  matrix  from  *"  and  subsequently 
combining  it  with  *'  (the  obvious  alternative  of  extracting  a  constant  matrix 
from  ♦'  and  combining  it  with  ♦"  also  apply). 

Fact  3.5.4:  Any  generalized  lossless  two-port  matrix  *"  with  p-1  can  be 
factored  into  the  product  of  two  matrices  and  of  the  same  type  i.e., 
where  -  constant,  and  i”  is  such  that  the  Y-polynomial  associated 
with  it  assumes  a  zero  value  for  £-0. 

Proof:  Let  k  -  Y"(£)/X"(0),  where  X"  and  Y"  are  corresponding  polynomials 
associated  with  and  be  defined  as  follows. 
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1 


K 


(3.5.11) 


Since  p  *  1,  it  respectively  follows  from  Properties  (P2)  and  (P4)  that  |Y"/X"| 

<  1,  wherever  X"*0  on  |z|-l  and  that  X"  is  scattering  Schur.  The  last  two 
properties,  due  to  an  extended  multidimensional  version  of  the  maximum  modulus 
theorem  proved  in  (93,  imply  that  either  Y"/X"  -  unimodular  constant  or  |Y"/X"| 

<  1  for  |z|  <  1.  In  the  former  case  we  would  have,  due  to  Property  (P2) 

satisfied  by  that  W"*Q,  which  is  impossible.  In  the  latter  case,  we  have, 
in  particular,  |k|  *  |Y"(  0)/X"(  0)  j  <  1.  Thus,  is  a  generalized  lossless 
two-port  matrix  with  p-1.  Since  it  is  easily  verified  that  (#J1)  also 
satisfies  this  property  it  follows  from  Fact  3.2.1  that  *"  -  (*£)  *"  is  also  a 

generalized  lossless  two-port  matrix.  Finally,  the  fact  that  Y-polynomial 
associated  with  *"  say  Y",  satisfies  Y"(0)  «  0  follows  from  Y"  - 

it  2  ^  ^  ^ 

(Y"-K  X" )// ( 1—  |  k |  )  (obtained  by  considering  the  (2,1)  entry  of  the  matrix 
equation  *”  -  (*")~1*n)  and  k  -  Y"(0)/XM( 0) . 


Next,  if  *  is  factorable  as  *  »  *'*"  then  due  to  Fact  5.4  we  may  also  write  *  » 
♦£*£,  where  »*'*£,  due  to  Fact  3.2.1,  is  a  generalized  lossless  two-port 
matrix.  Further,  it  is  trivially  verified  that  if  #'  and  #"  satisfies  the 
requirements  imposed  in  the  'Main  problem'  of  Section  3.2  then  and  *£  also 
satisfies  the  same  requirements.  Thus,  *  ■  is  a  valid  solution  to  the 
factorization  problem. 


If  ♦"  is  viewed  as  a  chain  matrix  0",  then  (assuming  that  the  operation  of 
shifting  the  factor  from  *"  into  *'  has  been  carried  out)  we  have  that  Y"(0) 
*  ®21^— ^  ®  the  corresponding  transfer  function  matrix  would  in  view  of 
(3.2.1)  satisfy  ^J^(O)  ■  0.  Consequently,  there  vnuld  be  no  direct  path  from 
'a'  to  'b'  via  *"  »  0"  in  Figure  3.2.  Similarly,  if  ♦"  is  viewed  as  a  hybrid 
matrix  T"  then  the  corresponding  trasnfer  function  matrix  would  satisfy  ^(.O) 
»  0,  thus  guaranteeing  no  direct  path  from  'a'  to  'b'  via  ♦"  -  P  in  Figure 
3.3.  In  either  case,  no  delay  free  loop  exists  at  the  junction  of  the  two 
two-ports.  Note  further  that  ♦"  as  in  (3.5.10)  correspond  to  the  chain  matrix 
or  the  hybrid  matrix  of  the  well  known  Gray-Markel  section. 
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Furthermore,  when  p-1,  X"  is  scattering  Schur  and  thus  X"(0)  *  0.  It  then 
follows  from  (3.2.12b)  that  if  Y(0)  -  0  and  Y"(0)  -  0  then  Y'(0)  -  0.  This  fact 
guarantees  that  the  prescribed  generalized  lossless  two-port  matrix  *  can  be 
successively  factored  into  product  of  generalized  lossless  two-port  matrices  of 
progressively  lower  complexity  in  such  a  way  that  the  Y-polynomial  associated 
with  each  of  the  factors  of  #  except  possibly  the  one  at  the  extreme  left  when 
*  is  interpreted  as  chain  matrix  6  is  equal  to  zero  for  z-0.  Similar 
considerations  apply  when  *  is  a  hybrid  matrix  r.  Absence  of  delay  free  loops 
at  each  junction  of  the  constituent  two-ports,  when  a  given  two-port  is 
fragmented  into  am  interconnection  of  more  elementary  two-ports  via  the  method 
of  factoring  *,  is  thus  guaranteed. 

The  algorithm  for  factoring  *  can  then  be  summarized  as  follows: 

Step  1:  If  the  prescribed  *  be  such  that  associated  XXz-  and  WWz-  are  coprime 
then  proceed  to  Step  2.  Otherwise,  factor  ♦’■*1*0*2  as  described  in  the  proof  of 
Proposition  3.3.2.  Replace  +  by 

Step  2:  Find,  if  possible,  two  linearly  independent  solution  {X!,  yr,  Y? ) ,  i«l, 
2  to  the  fundamental  equation  (3.2.17).  In  the  1-D  case  such  a  solution  always 
exist.  Factorization  of  #  is  inpossible  if  such  solutions  are  nonexistent. 

Step  2:  If  at  least  one  of  the  two  linearly  independent  solution  is  nonsingular 
i.e.,  xrxr  *  pY!Yi  for  any  i  then  proceed  to  Step  4.  Otherwise,  proceed  to  Step 
Step  5. 

Step  j4 :  Assuming  that  {XJ,  Y£,  YJ)  is  a  nonsingular  solution,  find  XJ  from  the 
second  equality  (3.4.1a)  and  (3.4.1c)  where  Y'  and  X'  are  replaced  by  Y!  and  XI 

2  1  2  1 

respectively.  Also,  find  from  (3.5.2)  and  a,  0  such  that  |a|  -  p|0|  - 
1/Kr  Finally,  form  X'-  «X£  +  /3pY'z-\  Y'  «  oY|  +  ex^z-' ,  Y"-  aY”  -  0p<r'*X£  and 
proceed  to  Step  6. 

Step  5>:  Find  the  constant  L  as  in  (3.5.3b)  and  p,  q  such  that  2Re(pq*L)  -  1. 
Form  X'  -  pX£  +  qX£,  Y'  -  pY£  +  qY£,  Y"  -  pY£  +  qY". 

Step  6:  Find  X"  from  (3.2.12a).  Thus,  {X',  Y'}  and  {X",  Y")  i.e.,  ♦'  and  ♦"  are 
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obtained. 


Step  7:  If  p-1  then  from  X",  Y"  associated  with  find  «»Y"(0)/X"(0) ,  ♦£  as 
in  (3.5.11),  *"  -  (♦j)-1*"  and  let  -  *'*£.  Thus,  *  -  without  delay  free 
loop  at  the  junction. 

Remark;  Since  Kx  in  (3.5.2)  is  a  real  constant  it  is  possible  to  choose  real 
values  of  a  and  0  such  that  the  right  hand  side  of  (3.5.1)  is  equal  to  1.  If  * 
is  real  (i.e.,  X,  Y,  W  have  real  coefficients  and  c-+l),  then  X£,  Y',  Y£  as  a 
solution  to  the  fundamental  equation,  and  thus,  X',  Y',  Y"  must  also  have  real 
coefficients  if  a  and  0  are  chosen  to  be  real.  Since  this  implies  that  the 
coefficients  of  X"  are  real,  the  factors  and  *"  would  then  also  be  real. 
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3.6.  CNE-DIHQJSICNAL  SYNTHESIS  AS  A  SPECIAL  CASE: 


In  the  one-dimensional  case  i.e.,  if  k  -  1,  a  closer  examination  of  (3.5.10a,b) 
reveals  that  u  -  e  -  2,  and,  therefore,  there  are  two  more  unknown 
coefficients  them  the  number  of  linear  equations  in  the  set  of  linear 
simultaneous  equations  which  determine  the  solution  to  the  fundamental 
equation.  Thus,  there  are  (at  least)  two  linearly  independent  solutions  of  the 
fundamental  equation,  and  in  view  of  Theorem  3.5.2,  the  problem  of 
factorization  of  ♦  always  admits  of  a  solution.  Consequently,  structurally 
passive  synthesis  for  #  is  achieved  by  performing  a  sequence  of  further 
factorizations  of  ♦'  and  ♦"  into  the  same  kind  of  matrices  of  progressively 
lower  complexity  i.e.,  n£  <  n^,  n£  <  n^  (since  n£  +  n£  -  n^),  until  a  stage  is 
reached  when  each  of  the  resulting  matrices  cannot  be  factorized  any  further. 
This  latter  situation  corresponds  to  the  case  that  each  of  the  two-ports 
resulting  from  the  decomposition  satisfy  n^  «  1,  i.e.,  degj  C  <  1  and  deg^  B  < 
1  and  deg^^  A  <  1.  However,  if  the  prescribed  *  is  such  that  X,  Y,  W  have  real 
coefficients,  o-+l  and  W  has  complex  roots  then  it  is  necessary  to  allow 
two-port  sections  with  n^-2,  deg^W-2  if  realization  involving  only  real 
multipliers  are  sought.  In  order  to  avoid  delay  free  loops  at  the  junction  of 
the  two-ports  we  further  require  that  two-port  sections  satisfy  Y(0)«  0  when  p 
■  1  or  equivalently,  B  -  0  when  ♦  -  6  and  C  -  0  when  ♦  -  r  for  Zj  -  0. 
Two-port  sections  of  the  above  types  will  be  called  elementary  sections  and  can 
in  turn  be  realized  in  structures  possibly  other  than  those  considered  here  by 
exploiting  synthesis  techniques  as  discussed,  for  example,  in  (4).  Thus,  the 
following  elementary  sections  are  obtained. 

An  arbitrary  lossless  chain  matrix  9  with  -  1,  B(0)  *  0  can  be  synthesized 
by  using  the  procedure  described  in  (4]  in  a  structure  given  in  Figure  4, 
whereas  an  abitrary  lossless  chain  matrix  9  with  n^  -  2,  B(0)  «  0  can,  by 
following  the  same  procedure,  be  synthesized  in  the  structure  of  Figure  3.5. 

A  lossless  hybrid  matrix  with  n^  -1  and  -  2  (assuming  Y(0)  *  C(0)  *  0)  can 
respectively  be  realized  by  the  same  elementary  sections  as  described  in  the 
Figures  3.4  and  3.5,  but  after  a  clockwise  rotation  of  the  corresponding 
diagrams  by  an  angle  of  90  degrees. 


On  the  other  hand,  an  arbitrary  lossless  transfer  function  matrix  l  (in  this 
case,  we  may  not  assume  Y(0)  to  be  zero)  with  ■  1  or  n^  ■  2  can  be  realized 
by  using  the  elementary  sections  described  above  only  after  a  Gray  -Markel 
section  has  been  extracted  from  the  corresponding  chain  matrix  (or  hybrid 
matrix)  so  as  to  effect  a  zero  value  for  B(0). 

Thus,  an  arbitrary  lossless  two-port  can  indeed  be  synthesized  as  an 
interconnection  of  Gary-Markel  sections  and  the  sections  depicted  in  Figures 
3.4  and  3.5  only.  Note  that  sections  of  Figure  3.4  and  3.5  were  introduced 
Dewilde  and  Deprettere  in  the  context  of  cascade  synthesis  (19J,  and  can 
viewed  as  scaled  versions  of  interconnections  of  wave  digital  filter  adapters 
[181. 
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3.7.  DISCUSSIONS  AND  ILLUSTRATIVE  EXAMPLES: 

The  purpose  of  this  section  is  many  fold.  First,  we  examine  the  structure  of 
the  fundamental  equation  in  somewhat  more  detail  to  facilitate  the  method  of 
solution  both  for  the  1-D  and  the  k-D  case.  Although  it  has  been  remarked  in 
Section  3.5  that  the  ^-matrix  is  generically  not  factorable  in  multidimensions, 
the  possibility  of  synthesis  for  special  classes  of  *-matrices  may  not  be  ruled 
out.  Furthermore,  in  k-D,  nonfactorability  of  any  one  of  the  three  matrices,  2, 
©,  or  T  associated  with  a  lossless  two-port  does  not  rule  out  the  factorability 
of  other  two  matrices.  This  fact  is  next  substantiated  via  examples,  thus 
justifying  the  need  to  study  factorization  of  all  three  types  of  matrices  (in  a 
unified  manner).  Finally,  in  practice,  all  multidimensional  frequency  filtering 
problems  require  some  form  of  symmetry  in  the  k-D  frequency  response,  and  it  is 
known  that  such  syrmnetries  dictate  that  the  two-port  be  either  symmetric  or 
(quasi)  antimetric  in  the  sense  of  classical  network  theory  (to  be  made  precise 
later  in  this  section).  This  is  indeed  the  case,  for  example,  in  the  design  of 
2-D  fan  [201  and  k-D  circularly  symmetric  [211  wave  digital  filters  based  on 
transformations  from  analog  prototypes.  Therefore,  the  factorabilty  of  the 
E-matrix  associated  with  these  subclasses  of  discrete  lossless  two-ports  is 
also  undertaken  in  the  present  section. 

A.  Computational  considerations: 

For  the  purpose  of  the  present  discussion,  the  following  notations  will  be 
adopted. 

X(z)-2  P.(z')zJ,  Y(z)-E  Q.(z')zJ,  poW'W'z-'-E  R.(z' )zj  (3.7.1  a, b,c) 

X'(z)-E  P?(z')zJ,  Y'(z)-E  Qf(z')zj,  Y"(z)— E  R!(z')zJ  (3.7.2  a,b,c) 

where  z'  is  the  k-tuple  of  integers  (z2,  z^.-.z^)  and  P.-P^z'),  Qi-Qi(z')» 
Ri*Ri(z')  and  Pf*PJ(z'),  Q["Q[(z')»  Rl-Rf(z')  are  polynomials  in  z' .  Then  the 
fundamental  equation  (3.2.17)  can  be  written  in  the  form  of  (3.7.3): 

V(z')T(z')«0,  T-T(z')-(tJ  I  T*  I  T*]fc  (3.7.3  a,b) 
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V(z')  - 


-P' 
*  *0 


I  #  •  •  •  Rq  I 


(3.7.3  c) 


where  the  superscript  't'  denotes  matrix  transposition,  is  a 

(n£+l)x(n^+n£+l)  lower  shift  matrix  whose  first  row  is  the  (k-1)  variable 

polynomial  matrix:  [P  , . ..,Pn,  0,...0]  and  subsequent  rows  of  which  are 

n1  u 

obtained  by  shifting  the  previous  rows  by  one  step  towards  the  right.  The 
(n£+l)x(nj+nj+l)  matrix  T2  and  the  (n£+l)x(nj+n£+l)  matrix  are  similarly 
obtained  from  the  polynomial  row  vectors:  , ...Qq,  0,...0]  and  tR2nJ  — R0- 

0,...0]  respectively.  For  a  given  T(z' )*  a  solution  V(z' )  to  f3.7.3a) 
corresponds  to  a  solution  to  the  fundamental  equation  if  in  (3.7.3c)  the 
following  degree  restrictions  for  all  v  are  satisfied  for  each  i«2  to  k. 


deg^  <  nr,  degj-^  <  nr,  deg.E"  <  n?  (3.7.4) 

As  remarked  earlier,  if  k>l  for  a  given  T(z)  a  solution  V(z)  satisfying 
(3.7.3a)  and  (3.7.4)  may  not,  in  general,  exist.  However,  the  following 
approach  may  be  adopted  in  attempting  a  viable  solution.  First,  find  the 
Hermite  reduced  form  H( z' )  of  T( z' )  via  the  pseudo-division  algorithm  as 
described,  for  example,  in  (22],  [23]  i.e.,  find  a  unimodular  matrix  U(z' )  such 
that  U(z' )T(z' )  «  H(z' )  is  in  Hermite  form.  Since  T(z/),  and  thus,  H(z' )  has 
(n^+n^+1)  columns,  but  (n^+n^+3)  rows,  the  last  k  rows  (k>2)  of  H(z' )  are 
identically  zero.  Consequently,  each  of  the  last  k  rows  of  U(£' )  belong  to  the 
left  null  space  of  T(z').  However,  in  order  forany  vector  belonging  to  this 
space  to  correspond  to  a  solution  of  the  fundamental  equation,  (3.7.4)  must  be 
satisfied,  which,  generically,  may  fail.  To  elaborate  further  on  this  it  may  be 
remarked  that  (proof  oranitted  for  brevity)  in  the  special  case  of  2-D  i.e., 
when  k-2  and  if  n^-n^  a  necessary  and  sufficient  condition  for  the  existence  of 
two  linearly  independent  solutions  to  the  funamental  equation  i.e.,  that  of 
factorability  of  ♦  is  that  the  dimensionality  of  left  null  space  of  Tfz'J-T^) 
be  exactly  two  and  the  two  left  Kronecker  indices  (23]  of  the  polynomial  matrix 
T(z')*T(z2)  be  each  equal  to  n^-nij.  We  have  so  far  been  unable  to  establish  an 
analogous  characterizaion  of  factorability  when  n^/n !j. 


In  the  1-D  case  (k-1),  however,  bo-Ji  V-V(z')  and  T-T(z')  in  equation  (3.7.3a) 
are  constant  matrices.  Furthermore,  since  it  is  known  that  the  lower  shift 
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matrices  T^,  and  are  closely  related  to  Toeplitz  as  well  as 

resultant-like  matrices,  the  linear  simultaneous  equation  (3.7.3a)  can  be 
potentially  solved  by  exploiting  recently  developed  fast  algorithms  for  solving 
such  equations  124).  Furthermore,  it  can  be  shown  by  pursuing  the  proof 
technique  for  Theorem  3.5.3  that  in  obtaining  two  linearly  independent 
solutions  to  the  fundamental  equation,  one  of  the  zeros  of  the  polynomial 
Y"(z^)  may  be  chosen  arbitrarily.  Once  such  a  choice  is  made,  the  solution  to 
the  fundamental  equation  becomes  essentially  unique  except  for  a  constant  scale 
factor  multiplying  each  of  the  polynomials  [X',  Y',  Y")  in  the  solution.  It  can 
be  shown  that  each  of  these  three  polynomials  in  this  solution  can  in  turn  be 
expressed  via  closed  form  determinantal  formulas  as  discussed  in  [25].  From  a 
ccmputational  standpoint  this  latter  method,  as  opposed  to  the  Toeplitz-like 
method  mentioned  above  may  not,  however,  be  the  most  inexpensive  when  the 
integers  nj,  n£,  and  thus,  n  are  large. 

B.  Examples  on  factorability  of  E,  0  and  T: 

We  next  illustrate  by  three  examples  that  in  multidimensions  by  viewing  the 
♦-matrix  as  three  different  types  of  matrices  associated  with  a  two-port, 
namely  the  transfer  function  matrix  E,  the  chain  matrix  ©  and  the  hybrid  matrix 
T,  the  synthesis  of  a  larger  class  of  discrete  lossless  two-ports  can  be 
attempted  than  by  considering  the  factorization  of  a  matrices  of  only  one  of 
the  above  kinds. 

(I)  Consider  a  discrete  lossless  two-port  described  by  A  »  PQ,  B,  and  C  -  2RS, 
n  -  (3,3)  and  o»- 1  as  in  (3.2.1),  where 

2  2 

,  Q-z^Z2~2z^Z2+Z2+4 ' 

B»4z^Z2-4z^Z2-5z^Z2-6z^Z2+7z^Z2-2z^Z2+4z^Z2-6z^Z2-4z2-z^-2z^-5z^+4, 

2  2  2  2 

R-Zj^-l >  S-Z^Z2-Z2+Z2Z2+^Z1Z2~Z2-Z1~Z1-^ 

In  attempting  a  nondegenerate  factorization  of  the  corresponding  ©-matrix  (p*l) 
into  nonconstant  0'  and  0"  we  encounter  the  following  distinct  possibilities: 
(i)  n'-(2,2),  n"-(l,l),  Y'-  C'«  S,  Y"-  C"-  2R  (ii)  n'-(l,l),  n"-(2,2),  Y'»  C'  - 
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R,  Y"»  C"  -  2S.  Neither  in  case  (i)  nor  in  case  (ii)  we  have  two  independent 
solutions  to  the  fundamental  equation  (3.2.17).  Via  Theorem  3.5.3,  we  thus 
conclude  that  0  cannot  be  factored  as  ©'©".  However,  the  corresponding  L  can 
factored  as  Z  -  E'E"  where  Z'  and  Z"  are  described  as  (clearly,  p«p'«p"— 1  in 
this  case) : 

W'«  A'  *  P,  X'=  B'»  2z1Z2-2z1z2-2z2+2,  Y'  =  C'=  z2+2z2+1,  <t'*-1 

W"»  A"  -  Q,  X"-  B"=  2z^z2-2z1z2-2z1+2,  Y"=*  C"=  z\+2z^l,  a"»-l 

Since  the  above  ©  can  also  be  viewed  as  a  hybrid  matrix  r  (with  slight 
modifiation  in  the  sign  of  a),  the  example  also  demonstrates  that  there  exists 
discrete  lossless  two-ports  for  which  the  associated  transfer  function  matrix 
can  be  factored  but  the  associated  hybrid  matrix  may  not  be  factorable. 

(II)  Consider  the  discrete  lossless  two-port  given  by  A,  B  -  2PQ,  C  -  RS, 
n-(2,2),  cr»l,  where 

2  2  2  2 

A»3-z1-z2-z1z2,  P-l-z^,  O-l+z^Zj,,  R«z1z2+z2+z1+1,  S-z^-z^z^l 

An  attempt  to  factor  r  into  nondegenerate  nonconstant  factors  F  and  r"  gives 
rise  to  the  following  two  distinct  possibilities  with  n'-n"«(l,l)  in  both 
cases:  (i)  Y'«  B'«  2P,  Y”*  B"»  Q  (ii)  Y'-  B'-  Q,  Y”«  B"=  2P.  In  neither  of  the 

above  two  cases  the  fundamental  equation  is  found  to  have  two  linearly 

independent  solutions,  thus  proving,  in  view  of  Theorem  3.5.3,  the 
impossibility  of  the  intended  factorizaion.  However,  the  corresponding  ©-matrix 
can  be  factored  as  ^>©'0",  where  ©'  and  0"  are  described  withy  p~ p'«p"»l  as: 

X'-  A'-(z1z2+2z1+2z2+3)//3,  Y'«  B'-(z1+z2+2)//3,  W'»  C'«  R,  <t'«1 

X"-  A"-(z1z2-2z1-2z2+3)//3,  Y”«  B”-(-2z1z2+z2+z1)//3,  W"-  C"-  S,  e"«l 

Note  that  Y"(0)»  B"(j))  -  0;  thus,  there  is  no  delay  free  loop  at  the  junction 
of  the  discrete  two-ports.  By  interchanging  the  roles  of  B-polynomial  and  the 
C-polynomial  in  the  above  example  we  can  similarly  demonstrate  the  existence  of 
a  discrete  lossless  two-port  for  which  the  T-matrix  can  be  factored  as  T-FT", 
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but  the  associated  ©  is  not  factorable  as  ^©'©". 

(Ill)  Consider  next  a  discrete  lossless  two-port  described  by  A,  B,  C-RS, 
n-(2,2),  <r-l,  where 

2  2  2  2  2  2  2 

A—  ( 7+3Zj-2z2-2z^Z2-Z2+3ziz2 ) >  B=ziz2'‘^2lz2+zl+zl-^zlz2~^zl22-^ ' 

2 

R-z^+1,  S—  2z^Z2“ 2z^ Z2”" 2z2^*2 

A  detailed  examination  of  the  degrees  of  the  polynomials  A,  B  and  C  reveals 
that  the  only  possible  way  of  factoring  the  transfer  function  matrix  I 
associated  with  the  two-port  is  to  attempt  either  (i)  n'-(l,0),  n"-(l,2),  A'-l, 
A"-A  or  (ii)  n'-(l,2),  n"-(l,0),  A'-A,  A"-  1.  In  both  cases,  however,  the 
fundamental  equation  (3.2.17)  fails  to  yield  two  linearly  independent 
solutions.  Thus,  I-matrix  associated  with  the  discrete  two-port  under 
consideration  cannot  be  factored.  On  the  otherhand,  the  associated  chain  matrix 
0  can  be  factored  as  ©■©'©",  where  ©'  and  ©"  are  described  as: 

a"-l,  a'-l,  W'«  C'-  R,  W"-C"«S,  X'-  A'— (z1+7)//15,  Y'-  B'-(2z1~4)//l5, 

X"-  A"-(4z^Z2-8z^Z2+4z^-Z2-2z2+15)//15,  Y"-  B"-(-z^Z2+4z2+2z^Z2+8z2-z^)//15 

Note  again  that  Y"(0)-B"(0)  ensures  that  there  is  no  delay  free  loop  at  the 
junction.  By  viewing  the  chain  matrix  as  a  hybrid  matrix  the  same  example  with 
minor  modofications  can  be  used  to  show  the  existence  of  a  discrete  lossless 
two-port  for  which  the  transfer  function  matrix  I  cannot  be  factored,  although 
it  is  possible  to  factor  the  associated  hybrid  matrix  r. 

C.  Symmetric  and  (quasi)  antimetric  two-ports: 

A  discrete  lossless  two-port  will  be  called  symmetric  or  (quasi)  antimetric  if 
(3.7.5)  (or  (3.7.5'))  holds  true.  Note  that  the  former  case  corresponds  to  fan 
type  symmetry  (20],  whereas  the  latter  case  correponds  to  circular  symmetry 
(21]  in  frequency  response. 

OoCz-,  B— aBz-,  (O-cCz-,  B—  aBz-)  (3.7.5)  ((3.7.5')) 
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Let  the  rational  functions  and  Lj  be  defined  as  in  (3.7.6)  (or  (3.7.6'))  for 
symmetric  (or  (quasi)  antimetric)  filters. 

Lj-tB+O/A  (or  {B+jC)/A  )  (3.7.6a)  ((3.7.6'a)) 

L2-(B-C)/A  (or  (B-jC)/A  )  (3.7.6b)  ((3.7.6'b)) 

Then  it  clearly  follows  from  (3.2.4c)  that  L^L^«L2L2*1,  and  consequently, 

| Li | * | L2 ) “l  for  |z|*l,  wherever  or  L2  are  well  defined.  Since  both  L^  and  L2 
have  scattering  Schur  denominators  [9],  it  follows  that  and  L2  are  both 

multidimensional  discrete  all-pass  functions.  By  making  use  of  (3.7.6)  (or 

(3.7.6')),  E  as  in  (3.2.1)  can  be  expressed  as  in  (3.7.7a)  and  (3.7.7b)  in  the 
symmetric  case  and  as  in  (3.7.7'a)  and  (3.7.7'b)  in  the  (quasi)  antimetric 
case. 


b^+L2 )/2 ,  E^2“^2l”^ Ll~^2 (3.7.7  a,b) 

b^+L2 ^ '  ^21m-^12=^ ^ ^2~^1  ^  (3.7.7'  a,b) 

Conversly,  for  any  two  discrete  all-pass  functions  L^  and  L2  the  matrix  E 
obtained  by  using  (3.7.7)  (or (3. 7. 7'))  is  a  valid  transfer  function  matrix  of  a 
discrete  lossless  two-port.  Thus,  any  multidimensional  discrete  lossless 
two-port  can  be  equivalently  described  by  means  of  two  multidimensional 
all-pass  functions  L.,  L2<  We  then  have  the  following  important  result. 

Theorem  3.7.1:  Let  { L^ ,  L2),  {L£,  L^} ,  {L£,  L£}  be  the  all-pass  functions 
associated  with  symmetric  (or  (quasi)  antimetric)  discrete  lossless  two-port 
transfer  function  matrices  E,  E'  and  E"  respectively.  Then  E  -  E'E"  if  and  only 
if  Lj  «  LjL£  and  L2  *  L^L£  hold  true. 

Proof:  Expressing  E  «  E'E"  in  terms  of  the  corresponding  L^  and  L2  via  (3.7.7) 
(or  (3.7.7'))  and  its  counterparts  for  E'  and  E",  it  follows  that  the 
factorability  condition  E  -  E'E"  is  equivalent  to  L^  *  L£L£,  Lj  *  L2L2' 

Since  it  can  be  easily  shown  by  pursuing  methods  outlined  in  [9]  that  any 
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rational  all-pass  function  L  can,  in  fact,  be  expressed  as  L-(P/P)z-  in 
irreducible  form,  where  P  is  a  scattering  Schur  polynomial,  Theorem  3.7.1 
conveys  the  important  fact  that  the  factorability  of  multidimensional  symmteric 
or  (quasi)  antimetric  discrete  lossless  two-port  transfer  function  matrices  can 
be  simply  expressed  in  terms  of  factorability  of  two  scattering  Schur 
polynomials. 


3. 8.  SUMMARY: 


The  present  work  has  been  motivated  by  the  possibility  of  designing 
structurally  passive  multidimensional  digital  filters.  A  simple  algorithm 
involving  the  examination  of  solution  ot  a  set  of  linear  simultaneous  equations 
for  studying  the  synthesizability  of  an  arbitrary  multidimensional  discrete 
lossless  two-port  has  been  derived  via  factorization  of  the  associated  chain 
matrix  0,  hybrid  matrix  T  and  transfer  function  E  by  introducing  a  generalized 
lossless  two-port  matrix  *,  which  in  turn  can  be  considered  as  a 
multidimensional  version  of  the  sigma-lossless  transfer  functions  discussed  in 
the  1-D  literature.  It  turns  out  that  under  a  generic  situation,  synthesis  not 
be  feasible.  In  the  special  case  of  one-dimension  our  algorithm  provides  new 
methods  of  realizing  structurally  passive  filters  directly  in  the  digital 
domain.  Although  in  the  multidimensional  (k>l)  case  synthesis  may  not  be 
feasible  for  an  arbitrary  discrete  lossless  0,  T,  and  E,  the  possibility  of 
synthesis  for  special  classes  of  discrete  lossless  two-ports  is  by  no  means 
ruled  out.  Examples  of  such  subclasses  of  two-ports  such  as  the  symmetric  or 
the  (quasi)  antimetric  discrete  lossless  two-ports  have  been  discussed. 
Existence  of  other  classes  of  discrete  lossless  two-ports  admitting  synthesis, 
albeit  in  special  topological  structures,  seems  feasible,  but  remains  to  be 
identified.  This  is  especially  true  in  view  of  cascade  synthesizability  of 
certain  classes  of  two-dimensional  continuous  time  systems  arising  in  studies 
of  lumped-distributed  netwoks  (22).  It  may  be  noted  that  the  cascade 
synthesizability  of  lumped-distributed  networks  can  be  characterized  in  terms 
of  properties  of  certain  polynomial  matrices  having  the  structure  of 
bigradients  (otherwise  called  resultants).  The  occurrence  of  polynomial 
matrices  of  similar  type  has  been  noted  in  our  study  in  the  context  of 
computing  a  solution  to  the  fundamental  equation  (cf.  equation  (3.7.3)). 
However,  further  investigation  is  needed  to  explore  this  connection  in 
successfully  utilizing  the  results  of  lumped-distributed  network  theory  in 
multidimensional  digital  filter  design. 
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QPlPTER  4 

SYNTHESIS  AtO  DESICM  OF  STRUCTURALLY  PASSIVE 
FOLLY  RECURSIVE  2-D  DIGITAL  FILTERS 

4.1.  INTRODUCTION 

Various  recursive  schemes  have  been  proposed  in  the  multidimensional  (m-D) 
digital  filter  literature.  Among  these  the  most  widely  studied  are  the  quarter 
plane,  the  asymmetric  and  the  symmetric  half-plane  recursive  scheme.  More 
recently,  motivated  by  needs  for  parallel  processing  of  2-D  signals  a  scheme 
known  as  the  fully  recursive  half-plane  scheme  has  been  proposed  in  [15],  and  a 
method  of  designing  transfer  functions  of  filters  having  this  recursive 
structure  has  been  outlined  in  [5].  The  impulse  responses  of  the  class  of 
filters  just  mentioned  satisfies  the  characteristic  property  that  the  region  of 
support  is  a  half-plane  and  the  filter  is  recursive  in  both  horizontal  and 
vertical  direction.  More  specifically,  the  recursion  equation  describing  the 
relation  between  the  input  x  and  output  y  of  a  filter  of  this  type  is  given  by: 

M{yn(n»))]«-  Lb[{y  .(m)}]  +  ^  L*  [  [x.  , .  (m) } )  (4.1.1) 

0  11  i«i  i  n  j.0  1  m-i; 

where  xn(m),  yn(m)  denote  the  n-th  row  of  the  input  and  the  output  signal? 
thus,  for  example,  x  (m)  -  x(m,n)  for  m  -  0,  +  1,  +  2,... etc;  and  the  (row) 

.  **  i  —  — 

operations  L?[ .  ]  and  L“[ . ]  respectively  denote  1-D  linear  shift  invariant 
convolution  operations  with  the  1-D  sequences  a.(m)  and  b. (m) .  Considering  the 

*  ^  a  k 

2-D  Z-transform  of  (4.1.1),  and  assuming  that  the  operations  L*[.)  and  l“( .  ] 
are  all  rational  we  then  have  (4.1.2)  for  the  transfer  function  of  the  filter, 
where  the  rational  functions  A^(z^),  B^(z^)  representing  the  row  operations 
just  mentioned  are  expressed  in  irreducible  rational  form  (i.e.,  as  the  ratio 
of  two  relatively  prime  polynomials)  as  in  (4.1.3). 


Y(z. ,z-)/X(z. ,z~)  -  ^  A. (z.  )zi/  r0  B. (z, ) 

1  1  1  *  i-0  1  1  1  i-Q  1  1 


(4.1.2) 


Ai(zl)  "  i  -  0,1,...!^ 


(4.1.3a) 


Bi(z1)  -  nJ(z1)/cJ(z1),  i  -  0,1,.. .Lp 


(4.1.3b) 
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On  the  other hand  it  is  now  well  known  that  an  input-output  description  such 
as  the  one  expressed  in  (4.1.2),  (4.1.3)  is  not  enough  for  the  successful 
operation  of  a  digital  filter  but  structural  considerations  need  to  be  taken 
into  account.  The  class  of  structurally  passive  filters  variously  known  as  the 
wave  digital  filters  (16),  orthogonal  filters  [17]  or  the  lossless  bounded  real 
filters  [18],  when  properly  designed,  are  known  to  satisfy  the  properties  of 
insensitivity  to  coefficient  perturbation  and  non-linear  arithmatic  conditions 
resulting  from  overflow,  finite  precision  arithmatic  etc.  Although  much 
progress  has  been  documented  in  the  synthesis  and  design  of  1-D  structurally 
passive  filters,  methods  for  two  and  higher  dimensions  are  still  evolving. 
Synthesis  methods  for  two  and  multi-dimensional  wave  digital  filters,  which  are 
quarter  plane  type  filters  have  been  reported  in  [16],  [7],  Quarter  plane  and 
asymmetric  half-plane  generalizations  of  1-D  lattice  filters  which  are,  in 
fact,  structurally  passive,  have  been  discussed  recently  in  the  context  of 
random  field  modeling  in  [11], [19]. 

Following  1-D,  in  the  present  paper  (pseudo)  passive  or  (pseudo)  lossless 
fully  recursive  half-plane  2-D  digital  filters  are  introduced  and  a  method  of 
their  structurally  passive  synthesis  and  subsequently  that  of  their  design  is 
discussed  for  the  first  time.  The  problem  of  synthesis  of  quarter  plane  causal 
(thus,  including  filters  causal  in  a  convex  cone  [20])  structurally  passive 
multidimensional  filters  of  the  type  mentioned  above  can  be  equivalently  viewed 
as  the  classical  network  theoretic  problem  of  synthesizing  a  lossless  but 
otherwise  arbitrarily  prescribed  multidimensional  transfer  function  as  an 
interconnection  of  elementary  building  blocks  such  as  capacitors  and  inductors. 

This  latter  problem  is  completely  unresolved  in  multidimensions  (m>2),  whereas 
in  2-D  synthesis  is  feasible  only  in  an  unconstrained  topological  structure 
[20].  On  the  otherhand,  it  has  been  shown  that  if  certain  ladder-like 
constraints  are  imposed  on  the  structure  in  which  the  filter  is  to  be 
synthesized  then  the  prescribed  2-D  transfer  function  must  satisfy  further 
restrictions  in  addition  to  input -output  losslessness  [21],  [22],  [23]. 
Related  other  synthesis  results  [7],  [16]  in  this  context  deal  with  important 
special  cases  when  the  multidimensional  frequency  response  of  the  filter 
possesses  certain  symmetries.  In  contrast,  the  present  work  provides  us  with  a 
synthesis  of  arbitrary  lossless  fully  recursive  half-plane  2-D  filters. 
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Additionally,  unlike  the  quarter  plane  case  referred  to  earlier  the  synthesis 
is  obtained  in  a  fixed  predetermined  structure  potentially  useful  for  practical 
i  nplemen ta t i on . 

As  in  most  passive  or  lossless  filter  design  techniques  our  synthesis 
method  proceeds  by  viewing  the  prescribed  passive  transfer  function  as  being 
embedded  into  the  transfer  function  of  a  lossless  two-port.  The  synthesis  of 
this  fully  recursive  half-plane  lossless  two-port  takes  advantage  of  a  recent 
algorithm  for  the  design  of  structurally  passive  1-D  filters  advanced  by  Rao 
and  Kailath  [6]  as  an  extension  of  the  celebrated  Schur  algorithm  [9).  Unlike 
all  other  methods  known  for  the  synthesis  of  1-D  continuous  as  well  as  discrete 
lossless  two-ports  including  those  available  in  the  classical  circuit  theoretic 
literature,  the  algorithm  of  (6]  enjoys  the  unique  feature  that  given  a 
transfer  function  associated  with  the  lossless  two-port  the  synthesis  algorithm 
makes  use  of  rational  arithmatic  operations  only  (i.e.,  nonrational  arithmatic 
operations  such  as  polynomial  factorization  is  not  required)  [10].  The 
synthesis  method  for  fully  recursive  half-  plane  filters  to  be  presently 
described  fully  exploits  this  rational  character  of  the  1-D  algorithm  in  [6]. 
Although  the  details  of  the  method  differ  nontrivially  form  1-D  due  to 
considerations  characteristic  of  multidimensional  problems  (e.g.,  those 
utilizing  techniques  from  elementary  algebraic  curve  theory  [3],  [12]),  the 
synthesis  to  be  outlined  can  be  considered,  at  least  at  a  conceptual  level,  to 
be  a  generalization  of  the  result  in  [6]  to  two-port  transfer  functions  the 
coefficients  of  numerator  and  denominator  polynomials  of  which  belong  to  a 
field  of  rational  functions  (instead  of  the  field  of  rational  numbers).  From  a 
different  perspective  the  present  work  can  also  be  viewed  as  a  generalization 
of  1-D  Schur  algorithm  to  2-D  fully  recursive  half-plane  filters,  thus  making 
it  possible  to  cast  the  present  discussion  in  the  closely  related  framework  of 
modeling  of  stationary  random  fields  and  scattering  theory  [9]. 

A  note  regarding  the  stability  of  the  filter  is  in  order.  The  region  of 
analyticity  of  the  transfer  function  of  our  filter  will  be  found  to  marginally 
differ  from  those  previously  considered  in  the  2-D  half-plane  literature  [4], 
[5],  This  is  primarily  due  to  the  fact  that  the  results  such  as  those  in  [4], 
[5]  are  motivated  by  bounded- input-bounded-  output  considerations,  whereas,  in 
contrast,  our  results  are  driven  by  passivity  considerations.  The  fact  that 
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this  difference  in  consideration  does  indeed  lead  to  diverging  formulations  of 
stability  in  multidimensions  (m>l),  but  not  in  1-D,  is  now  known  [1],  [2]. 
Thus,  there  is  no  contradiction  between  our  stability  results  and  those 
existing  in  the  half-plane  literature  so  far. 

In  Section  4.2  the  fully  recursive  half-plane  passive  one-ports  are 
characterized  in  terms  of  their  transfer  function.  Similar  considerations  in 
the  context  of  two-port  transfer  functions  form  the  context  of  Section  4.3.  A 
representation  theorem  for  fully  recursive  half-plane  lossless  two-ports 
analogous  to  that  of  the  Belevitch  canonical  form  (8)  of  representation  for 
lossless  1-D  continuous  two-ports  of  classical  network  theory  is  developed 
here.  In  Section  4.4  the  synthesis  method  based  on  this  representation  theorem 
is  described,  and  in  Section  4.5a  design  methodology  is  proposed  by  taking 
into  account  the  symmetry  requirements  [14]  on  the  frequency  response  imposed 
by  many  practical  multidimensional  processing  tasks. 
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4.2.  FULLY  RECURSIVE  SYMMETRIC  HALF-PLANE  PASSIVE  SYSTEMS: 

By  associating  the  (pseudo)  energy  EE|x(nlfn2)  |2  to  the  input  x(n1#  n2)  of  the 
system  we  first  develop  conditions  necessary  for  the  transfer  function  of  a 
fully  recursive  half-plane  filter  to  be  passive  in  the  sense  that: 

EE  |y(nr  n2)|2  <  EE  |x(nr  n2)|2  (4.2.1) 

for  any  choice  of  square  sunmable  input  sequence  x(n^,  n2). 

By  choosing  x(n1#  n2)  -  ^{n^,  n2)  i.e.,  the  2-D  impulse  function,  the  impulse 
response  h(n1,  n2)  of  the  filter  can  be  obtained  as  the  corresponding  output. 
We  then  have: 

CD 

y(n1#  n2)  -  h(n1,  n2)  -  E  hj^nj)  S(n2~k)  (4.2.2) 

where  hk(n^),  k»0,l, . . .etc.  are  certain  1-D  sequences  and  6(.)  is  the  1-D 

impulse  sequence. 


Considering  the  z-transform  of  (4.2.2)  we  obtain: 

*  nl  k 

H(z,,  z,)  -  E  (E  h,  (n.  )z.  ]  z~  (4.2.3) 

1  1  k-0  nx  *  1  1  1 

Using  the  Schwartz  inequality  it  follows  from  (4.2.3)  that 

OB 

|H(z.,  z7)|2  <  k(z-)  .  E  | E  h.  (n. Jz.1!2  (4.2.4) 

1  Z  1  k-0  n:  *  1  1 

0  A 

where  k(z2)  -  1  +  |z2P+  \z2>  +  -  etc. 

If  we  consider  the  special  case  -  exp(jw1)  then  we  have  (4.2.5)  from 
(4.2.4). 


0«1  2  *2 

|H(e  \  z2)T  <  k(z2)  .  E  I  Hk  ( )  r 

k— 0 


(4.2.5) 


where  H^^)  is  the  Fourier  transform  of  \(n1)  for  each  k  -  0,1,2, ..  .etc. 
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On  the  otherhand,  substituting  z ^  -  expij^),  i  -  1,2  in  (4.2.3)  we  obtain  that 
the  Fourier  transform  Hfo^,^)  of  hfn^n^  is  given  by: 

H(w1,  u^)  -  I  Hk(«x)  «  (4.2.6) 

which  in  turn  yields  that: 


*  2  :Kn 

£  j H,  ( ox. )  |  +  terms  involving  e  with  n*0 

k-0 


(4.2.7) 


Assuming  termwise  integrability  of  the  right  hand  side  of  (4.2.7)  we  then  have: 

+rt+n  j  n  •  _ 

S  I  IHtw^M^rdw^  -  2n  J  {  £  H^cOn  d «.  (4.2.8) 

-n-n  -n  k*0 

However,  by  the  2-D  Parseval's  formula  the  left  hand  side  of  (4.2.8)  is  equal 
2  2 

to  (4n  ) £ | h ( n1 ,  n2)|  i.e.,  the  total  (pseudo)  energy  in  the  signal  hfn^,  n2), 

whereas  the  right  hand  side  can  be  similarly  interpreted  as  the  sum  of  the 

(pseudo)  energies  contained  in  the  row  outputs  hg(n^),  h^(n^) , . . .etc. 

Furthermore,  if  we  assume  the  fully  recursive  filter  under  consideration  to  be 

passive  then  from  (4.2.1),  first  part  of  (4.2.2),  and  the  fact  that  x(n.,o,)  « 

2  i  i 

&(n^,n2),  it  follows  that  E£  |h(n^,n2)|  <  •.  Thus,  the  integral  in  the  left 

hand  side  of  (4.2.8)  is  finite,  and  consequently,  except  possibly  isolated 
values  of  in  the  interval  [ — n, n ) ,  the  integrand  in  the  right  hand  side  is 
bounded  i.e.,  we  have : 

£  |H.(«.)|2  <  -  (4.2.9) 

k-0  *  1 

In  view  of  (4.2.5)  arid  the  fact  that  k(z2)  <  1  for  |z2|  <  1  we  then  conclude 
that  for  all  (z^(  -  1  except  possibly  isolated  values  on  the  unit  Circle  and 
for  all  |z2|  <  1,  H(z^,z2)  is  bounded;  and  furthermore,  if  H(e  z2Q)  is 

unbounded  for  some  |z-Q|  <  1  and  real  w^,  then  £  | Hk ( Q )  ] 2  must  be 

k*0 

unbounded,  and  thus  H(e  ,z2)  must  also  be  so  for  all  z2* 


We  digress  temporarily  to  examine  a  consequence  of  passivity  reflected  on  the 


transfer  function  H(z^,  z2)  of  the  filter.  By  using  the  2-D  Parseval's 

theorem,  (4.2.1)  and  the  fact  that  Yfw^u^)  *  Hfaj^tt^Xfo^,^) ,  where 
Xtca^,^)  are  the  respective  Fourier  transforms  of  y(n^,n2),  xfn^,^)  we  have 
that 

n  n  2  2 

/  /  (l-|H(o»1,w2)r)  dc^  d«2  >  0  (4.2.10) 

— n—  it 


Since  (4.2.10)  is  true  for  any  input  Xju^,^),  we  have  that  |H(w^, )  I  <1  for 
all  real  two-tuples  (o^,^)  except  possibly  for  finitely  many  of  them. 


3<o10 

Consequently,  if  H(e  ,z-QL  is 
then  as  shown  previously  H(e  ,z 
for  all  z2  on  |z2|  -  1  in  particular, 
the  conclusion  of  the  last  paragraph. 
1  and  lz-,1  <  1. 


unbounded  for  some  |z2q|  <  1  and  real  oo^, 
2)  would  be  unbounded  for  all  z2,  and  thus 
This  latter  situation  would  then  violate 
Thus,  H(zirz2)  is  bounded  for  all  |z^|  - 


We  next  make  the  simplifying  assumption  that  for  each  k  -  0,1,... etc. 

( z^ )  -  L  h^n^z^1  i.e.,  the  z-transform  of  h^n^)  are  rational  functions  in 

zr  Since  the  recursion  equation  for  the  fully  recursive  symmetric  half-plane 
filter  is  given  by  (4.1.1)  it  is  clear  that  the  transfer  function  H(zlfz2)  is  a 
rational  function  of  z2>  Under  the  present  assumption  however,  H(zirz2),  in 
view  of  (4.2.3),  becomes  as  in  (4.1.2)  and  (4.1.3)  a  rational  function  of  both 
z^  and  z2,  and  can  be  expressed  as  the  ratio  of  two  relatively  prime 
polynomials  b(z^,z2)  and  a(zlfz2)  as: 

b(z1,z2) 

H(z^,z2)  ■  — — — 

a(z1(z2)  (4.2.11) 

We  now  claim  that  for  passive  systems  presently  under  consideration,  the 
polynomial  a(z^,z2)  in  (4.2.11)  cannot  have  infinitely  many  zeros  on  the 
distinguished  boundary  |z^|  -  |z2)  «  1  of  the  unit  bi-disc.  For,  if  a(z^Q,z2Q) 
■  0  for  some  lz^g!  *  lz20^  “  *  t^,en  *n  v*ew  (4.2.11),  in  order  for 
| H( ,a>2 )  1  <  1  to  be  satisfied  we  would  need  bU^z^)  -  0  i.e.,  a(z1,z2)  and 
b(z^,z2)  would  have  a  common  zero  on  |z^|-|z2|«l.  However,  the  presence  of 
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infinitely  many  such  zeros  would,  in  view  of  Bezout's  theorem  in  algebraic 
curve  theory  [3],  require  that  b(z^,z2>  and  a(z^,z2)  have  a  common  factor, 
which  has  been  hypothesized  to  be  absent  in  (4.2.11). 


A  further  consequence  of  passivity  is  that  when  the  transfer  function  H(z^,z2), 
is  expressed  in  terms  of  ratio  of  two  relatively  prime  polynomials  b(zlfz2)  and 
a(zlfz2)  as  in  (4.2.11),  it  must  be  true  that  a(z^,z2)  is  relatively  prime  with 
the  polynomial  a(z^,z2)  defined  as: 


dl  d2 


*  *“1  *“1 


a(z1,z2)  m  a(z^,z2)z^  z2  \  a(z^,z2)  *  a  (z^  #^2  ^ 


(4.2.12  a,b) 


where  d^,  d2  are  the  partial  degrees  of  a  in  z^  and  z2,  and  *  denotes  complex 
conjugation.  To  prove  this  let  g(z^,z2)  be  the  gcd  between  a(z^,z2)  and 
a(z1,z2).  Then  as  shown  in  (2)  we  must  have  g(z^,z2)  -  yg(Zj.,z2)  where  r  - 
constant,  |r|“l;  and  g(zlfz2)  is  either  a  constant  or  must  have  infinitely  many 
zeros  on  |z^|  -  |z2|  -  1.  In  the  latter  case,  a(z1#z2)  must  have  infinitely 
many  zeros  on  | (  -  |z2|  -  1,  which  has  been  shown  to  be  impossible  in  the 
last  paragraph.  Thus,  g(Zj,,z2)  -  constant,  and  a(z^,z2)  is  relatively  prime 
with  a(zlfz2) . 


The  essential  features  of  the  above  discussion  are  summerized  in  the  following 
result. 


Property  4.2.1:  The  transfer  function  of  a  passive  fully  recursive  symmetric 
half-plane  filter,  when  expressed  in  irreducible  rational  form  as  in  (4.2.11), 
satisfies  the  following  two  conditions:  (i)  a(z2,z2)  *  0  for  l z^ |  -  1  and 
|z2l  <  1  i.e.,  H(zlfz2)  is  analytic  in  (z^l  -  1  and  |z2l  <  1.  ( i i )  a (z^,z2) 

and  a(Zj,z2)  as  defined  in  (4.2.12),  do  not  have  a  common  factor. 

The  condition  (ii)  in  the  above  property  can,  in  fact,  be  replaced  by  any  one 
of  the  conditions  expressed  in  the  following. 

Assertion  4.2.1:  Let  a(z^,z2)  be  a  polynomial  in  z^  and  z2  such  that  a(z^,z2) 
*  0  for  |ZjJ  -  1  and  |z2|  <  1.  Then  the  following  conditions  are  all 
equivalent. 
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(a)  a(z1#z2)  does  not  have  infinitely  many  zeros  on  the  distinguished  boundary 
|z1l  -  \z2\  -  1. 

(b)  a(z.,z2)  and  a(z^,z2)  are  relatively  prime  polynomials. 

(c)  Each  irreducible  factor  of  a(z^,z2)  has  at  least  one  zero  in  the  domain 

I z^ I  *  1 t  I z2 I  ^  1 • 

Proof:  It  has  already  been  proved  that  (a)  implies  (b).  To  show  that  (b) 
implies  (a)  observe  that  if  for  some  Iz^gl  *  |z2q|  “  a^z10'z20*  *  0  then 
from  (4.2.12)  a(zio,z20^  *  0.  Consequently,  if  a(zlfz2)  has  infinitely  many 
zeros  on  | z1 1  «  1  then  a(z^,z2)  and  a(zlfz2)  would  have  infinitely  many  common 
zeros  (on  | z^ |  «  |z2|  «  1).  Therefore,  due  to  Bezout's  theorem  (],  a(z^,?2) 

and  a(z. ,z2)  would  not  then  be  relatively  prime  polynomials.  Thus,  (a)  and  (o) 
are  equivalent. 

Next,  if  a1(z1,z2)  is  any  irreducible  factor  of  a(zlfz2)  then  obviously 
al*zl'z2*  *  0  or  'zl'  "  | z2 1  <  1.  Furthermore,  if  a1(z1,z2)  does  not 

contain  any  zero  in  |ZjJ  -  1,  |z2l  >1  then  for  any  z^  on  Iz^J  »  1,  a^z^z^  * 
0  in  |z2l  <  1  as  well  as  in  |z2l  >  1,  and  thus  the  values  of  z2  such  that 

al*zl'z2^  "  0  must  be  on  | z2 1  ■  1.  Consequently,  a^z^z^,  and  thus  a(z^,z2) 
would  have  infinitely  many  zeros  on  |z^j  «  |z2)  -  1.  Therefore,  (a)  (or 
equivalently  (b))  implies  (c). 

To  prove  that  (c)  implies  (b)  let  g  -  gcd(a,a)  i.e.,  a  -  g.e,  a  *  gf,  where  e 
and  f  are  relatively  prime  polynomials.  Then,  as  shown  in  (2]  g  -  yg  where  r 
is  a  constant.  Assuming  g  to  be  a  nonconstant  polynomial,  if  each  irreducible 
factor  of  a  contains  at  least  one  zero  in  |z1|  -  1,  | z2  j  >  1  then  g  and  thus  g 
-  yg  must  have  a  zero  in  |z^|  »  1,  |z2|  >  1.  However,  this  implies  that  the 
polynomial  g  and  thus,  in  view  of  a  ■  ge,  the  polynomial  a  must  have  a  zero  in 
|ZjJ  ■  1,  |z2l  <1,  which  is  a  contradiction.  Thus,  g  ■  constant  and  a  and  a 
are  relatively  prime. 

We  also  have  the  following  important  result. 

Property  4.2.2a:  If  a  rational  function  H-H(z^,z2)  as  in  (4.2.11)  is  such  that 


| H |  <1  on  |z^|-|z2(«l  excePt  possibly  at  finite  number  of  points  where  H  is  not 
well  defined  and  if  H  satisfies  the  conditions  expressed  in  property  4.2.1(i) 
(thus,  if  H(z2»z2)  is  tranfer  function  of  a  passive  fully  recursive  half  plane 
filter)  then  | H |  <1  for  all  lz^l-1  and  |z2|<l.  Furthermore,  for  some  ( z1o»z2o^ 
with  |z10|-l,|z20l<l.  Then  H(z1Q,z2)  is  a  constant  independent  of  z2. 
Assuming  H  to  be  nonconstant,  the  latter  situation  can  arise  for  at  most 
finitely  many  values  of  z^g  (with  jz^g(-l). 

Proof:  Due  to  property  4.2.1  (i)  the  denominator  polynomial  of  H(z^,z2)  cannot 
be  zero  for  some  fixed  lz^j-1  and  arbitrary  values  of  z2<  Thus,  if  for  any 

z1Q  with  1 2io I “1  we  define  hi“hi^z2^"h^z10'z2^  t^)en  c*ue  to  our  hypothesis  is 
well  defined,  analytic  in  |z2)<l  |H^|<1  for  all  )  z2 1  —1  (except  at  the 
possible  poles).  Thus,  by  maximum  modulus  theorem  | | <1  for  all  |z2l<l. 
Since  this  is  true  for  arbitrary  z^  on  |z1l-l  the  first  part  follows. 

To  show  the  second  part  assume  for  contradiction  that  for  some  |z^g|«l, 

I *20 1 <1 '  we  have  H^z10'z20*“**  Then  as  shown  above  the  maximum  modulus  theorem 
applies  to  Hj-H^ZjJ-H^z^.Zj)  and  thus  [H1(z20)|-1  with  |z2g|-l  implies  that 
Hi-H^Zj) -Oconstant  i.e.,  in  view  of  (4.2.1)  b(z1Q,z2)  -  C  a(z1Q,z2). 
However,  if  H  is  nonconstant  then  b(z^,z2>  and  a(zlfz2)  are  relatively  prime 
nonconstant  polynomials.  Consequently,  b(z1Q,z2)  and  a(z1Q,z2)  treated  as 
one-variable  polynomials  may  fail  to  be  relatively  prime  for  at  most  finitely 
many  values  of  z^g  [1).  The  second  part  is  thus  established. 

In  fact,  the  following  result  in  Projperty  4.2.2b  can  also  be  proved.  This 
result  shows  that  the  polynomials  of  the  type  described  in  properties  4.2.1  (i) 
and  4.2.1  (ii)  characterize  denomination  of  irreducible  rational  functions 
satisfying  the  (half-plane)  boundedness  property:  | H | <1  for  ) z^ | —1 ,  |z2|<l. 

Property  4.2.2b:  If  H  is  a  noncons taunt  irreducible  rational  function  as 
expressed  in  (4.2.11)  and  is  such  that  | H | <1  for  |z^|-l, |z2|<|  then  either  a  is 
a  constant  or  satisfies  properties  4.2.1  (i)  and  4.2.1  (ii). 

Proof:  Obviously  it  is  impossible  to  have  a-0  and  b*0  for  any  |z^|-l,  |z2|<l, 
because  otherwise  |H|  would  be  unbounded  there.  If  a-b-0  for  some  |z^g|-l, 

I z2q I <1  then  consider  an  arbitrary  small  arc  of  |z^|-l  issuing  from  z^g. 
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Let  I*2  be  the  continuous  [12]  arc  traced  out  by  z2  (beginning  fro®  z2Q)  such 
that  a(z1,z2)-0  is  satisfied.  Note  that  since  is  assumed  arbitrarily  small, 
due  to  the  continuity  property  of  zeros  of  a  polynomial  as  a  function  of  its 
coefficients,  I*2  must  lie  completely  within  |z2|<l.  We  next  claim  that  for 
Zjerj  there  must  exist  values  of  z2cr2  such  that  a“0,  b®0308®  otherwise  a 
and  b  would  have  infinitely  many  common  zeros,  which  due  to  Bezout' s  theorem 
[3]  violate  the  fact  that  a  and  b  are  relatively  prime  polynomials.  However, 
since  r1c{z1, |z^|-l]  and  r2c{z2; |z2l<l)  this  latter  conclusion  has  already  been 
shown  to  be  inpossible.  Thus  a*0  for  J z^ ) —1 , | z2 1 <1 . 

Finally,  if  a(zlQ,z2Q)-0  for  some  jz10l-|z20|-l  then  b(z10,z2Q)-0  because 
otherwise  |H|<1  would  be  violated  in  (z^-1,  | z2 1 <1  at  the  vicinity  of 
( zio,z20^ '  Thus  existence  of  infinitely  many  such  (z10'z20*  would  again 
violate  the  relative  primeness  of  a  and  b. 

We  next  assume  the  filter  to  be  (pseudo)  lossless  in  the  sense  described 
earlier  i.e.,  equations  (4.2.1)  and  (4.2.10)  are  satisfied  with  equality. 
Consequently,  from  (4.2.10)  we  then  have  that  for  all  2-tuples  (c^,^)  with  the 
possible  exception  of  finitely  many  (4.2.13)  holds  true. 

IHfWj^)  |  ■  1  (4.2.13) 

we  first  note  that  the  rational  transfer  function  H(zlfz2)  of  a  (pseudo) 
lossless  fully  recursive  half-plane  transfer  function  satisfies  the  property 
that 


H(z1,z2)  H(zlfz2)  -  1  (4.2.14) 

To  substantiate  this  result  we  observe  from  the  definition  of  the  operation 
that  H(z1,z2)  -  H*(z1#z2)  for  |z1(-|z2|-l.  Consequently,  from  (4.2.13)  it 

follows  that  H(z1,z2)-H-1(z1,z2)  for  all  2-tuples  (z^Zj)  on  |z1|-|z2l«l  with 
possible  exception  of  at  most  finitely  many  values.  Thus  the  two  variable 
rational  function  H(zlfz2)  and  H_1(zlfz2)  assume  equal  values  at  infinitely 
many  distinct  points  (z^Zj),  and  consequently,  due  to  analytic  continuation 
are  identically  same,  i.e.,  H(z^,z2)  -  H-1(z^,z2)  for  all  z^^  and  z2. 
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For  convenience  of  further  exposition  the  following  terminology  will  be 
introduced.  Any  rational  function  H(z^,z2)  as  expressed  in  (4.2.11)  will  be 
said  to  be  a  fully  recursive  half-plane  all-pass  function  if  H(zlfz2)  satisfies 
the  conditions  stated  in  property  4.2.1  and  inequation  (4.2.14).  Thus, 
transfer  functions  of  (pseudo)  lossless  fully  recursive  half-plane  filters  are 
fully  recursive  half-plane  all-pass  functions. 


A  function  A(z^,z2)  of  two  variables  z^,  z2,  when  expressible  as  a  polynomial 
in  z2  with  coefficients  as  rational  functions  in  will  be  said  to  be  a 
pseudopolynomi al  (in  z2).  Thus,  if  A(z^,z2)  is  a  pseudopolynomial  then 

N2 

A(z1,z2)  -  +  a1(z1)z2  +...+  oN^(z1)z2  (4.2.15) 

where  a^z^'s  are  rational  functions  in  z^.  With  a(z2,z2)  as  given  in 
(4.2.15),  where  (z^)  is  not  identically  zero,  the  integer  n2  will  also  be 
denoted  by  deg^.  ^  Furthermore,  the  notation  A(z1,z2)  will  be  used  to  denote 
the  pseudo-  polynomial  obtained  from  A(z^,z2)  as: 


M zl'z2^  m  A^zl'z2^z2 


(4.2.16) 


Two  pseudopolynomials  Bfz^Zj)  and  C(zlfz2)  are  said  to  be  coprime  if  there  is 
no  pseudopolynomial  D(z2,z2)  actually  involving  z2  such  that  B(z^,z2)  - 
D(z^,z2)  B1(z1,z2)  and  C(z2,z2)  -  D(z2,z2)  C^(z2,z2)  for  some  pseudopolynomials 
Bl*zl'z2*  301:1  Cl*zl'z2*'  The  ^ollowin9  property  then  holds  true. 

Property  4.2.3:  Any  fully  recursive  half-plane  all-pass  function  H(z^,z2) 
(thus,  rational  transfer  function  of  (pseudo)  lossless  fully  recursive 
half-plane  filter)  can  be  expressed  as  follows: 


H(z^,Z2)  *  — D(z^)  ( A( z^ , z2 )/A( z^ , Z2 ) ] 


(4.2.17) 


where 


i) 


A(z2,z2)  is  a  pseudopolynomial 


N 


ii)  D(z^)  «  z^y  (d(z1)/d(z1) ]  ,  where  d(z^)  is  a  polynomial  in  z^,  y 

is  a  constant  of  unit  modulus  and  N  -  integer. 


iii)  the  pseudopolynomials  A(z^,z2)  and  A(z^,z2)  are  coprime. 


iv)  A(z1,22)  *  0  for  all  |z^|-l,  |z2|<l. 

Conversely,  any  rational  function  expressible  as  in  (4.2.17)  with  (i),  (ii), 
(iii)  and  (iv)  in  force  is  a  fully  recursive  half-plane  all-pass  function. 


Proof  of  property  4.2.3:  Let  H(zlfz2)  -  A(z1,z2)/B(z1,z2) ,  where  A  -  A(z.,z2) 
and  B  -  b(z1,z2)  are  pseudopolynomials  expressible  as  A  *  a^/Ap  and  B  "  b^^/bp, 
where  in  turn  -  a^ Zj,z2 ) ,  bN  -  bN(z^,z2)  are  polynomials  in  both  z^  and  z2, 
whereas  -  a^z^)  and  bD  -  bD(z^)  are  polynomials  in  z^  only. 

We  further  assume  that  H  -  H(z^,z2)  expressed  as  in  (4.2.18)  is  in  irreducible 
rational  form  i.e.,  the  pairs  of  polynomials  (a^a^,  (b^bp),  (a^b^  and 
(bD,j^)  are  relatively  prime. 

H  *  (bDaN)/(aDbN)  (4.2.18) 


Then  from  (4.2.18)  equations  (4.2.19),  (4.2.20)  and  (4.2.21)  follows,  where  the 
generic  notation  n^  for  denoting  the  degree  of  the  polynomial  p  in  the  i-th 
variable  has  been  used. 


H  -  H'1-  (^,1/^,1 


V  -V 


(4.2.19) 


N 


n 


lb. 


+  n 


laN 


~  °'  vd  "  ni%  +  "u^  *  0 


(4.2.20  a,b) 


and  N.  -  n,_  - 

Since  H  is  analytic  in  IZjJ-l,  |z2|<l  30(1  neither  a^  nor  b^  can  have  a  factor 
z2,  it  clearly  follows  that  >  0.  Also,  since  H  in  (4.2.18)  is  in 
irreducible  rational  form,  it  follows  by  comparing  (4.2.18)  and  2.19)  that 


n 


2bN 


(4.2.21) 


°Vn  ■ 


“V’n  ■  Vm  H 


{4.2.22  a,b) 


where  a  «  a(z^,z2)  is  a  polynomial  in  z^  and  z2.  By  inserting  (4.2.22b)  into 
(4.2.18)  and  subsequently  making  use  of  the  relations  between  aj^  and  a^, 
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between  and  and  finally  by  using  A  -  a^a^  (4.2.20)  and  (4.2.21)  we 
obtain  the  following 


-  -(nUM  +  nUV.)  "2aN 

H  -  aI(aDbD)/(aDbD)HA>^].z1  N  z2 

2bN 


(4.2.23) 


n, 


By  defining  d  -  apbp  and  noting  the  fact  that  A  -  A  z2  we  then  have: 


H  -  «( d/d)  (A/A)  z1 


-(n,.  +  n1K  ) 
l3N  lbN 


(4.2.24) 


Since  H  in  (4.2.18)  is  irreducible  and  analytic  in  |z^ |-1,  | z2 1 <1  we  note  that 
a^bjj  cannot  have  a  factor  z2-  Invoking  this  fact  and  considering  the  *  of 
(4.2.22a)  we  then  have: 

“Vd  ■  SVi*  (4-2-251 


where  in  (4.2.25)  k  is  the  total  multiplicity  of  z^  in  (a^b^).  By  substituting 
(4.2.25)  into  (4.2.22b)  we  obtain  ao^z^D_k  (note  that  since  from  (4.2.20b)  vD  « 
degree  of  (apfc^)  in  Zj  it  obviously  follows  that  vD~k>0).  Consequently,  it  must 
be  true  that  a  is  a  monomial  involving  z^  only  i.e.,  is  of  the  form 


v  -k 

a-yz^  for  some  constant  y. 


Then  aa  ”  yy*  m  z 


V* 


(4.2.24)  yields  (4.2.17)  with  N 
4.2(ii)  are  thus  established.  To  show  tha' 


V'ni 


+  "i 
2(  iii 


'VS1 


.  Thus,  | r 1—1 •  Therefore, 
+  k).  Properties  4.2(i)  and 
holds  true  note  that 

(4.2.26) 


Consequently,  if  A  and  A  has  a  pseudopolynomial  cannon  factor  then  it  follows 
from  A^  -  aN/^aD  and  (4.2.26)  that  aj^  and  a^  must  have  a  common  factor  involving 
Zj.  In  view  of  (4.2.22a,b)  then  and  bD&^  would  not  be  relatively  prime, 
thus  violating  the  irreducibility  of  H  in  (4.2.18).  Finally,  to  prove  (iv) 
note  that  it  follows  from  property  4.2.1,  (4.2.18)  and  (4.2.22b)  that  and 
thus  a^  is  nonzero  for  |z^|  •  1,  and 

l*2l  <  1. 


The  converse  proposition  follows  trivially  from  the  fact  that  any  H  -  H(z^,z2) 
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satisfying  (4.2.17)  along  with  (i)  through  (iv)  is  necessarily  analytic  in  | z1 | 
■  1,  Izjl  <1  and  has  the  property  of  HH  ■  1  on  |z^|  -  Iz^l  -  1. 
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4.3.  FULLY  RECURSIVE  SYMMETRIC  HALF-PLANE  LOSSLESS  TWO-FORTS: 

Characterizations  of  fully  recursive  symmetric  half-  plane  passive  as  well 
as  lossless  one-port  filters  have  been  established  in  the  previous  section  in 
terns  of  the  transfer  function  of  the  filter.  In  this  section  we  make  use  of 
the  results  of  the  previous  section  to  characterize  fully  recursive  symmetric 
half-plane  lossless  two-ports.  In  particular,  a  convenient  representation  for 
such  two-ports  analogous  to  the  Belevitch  canonical  representation  of 
continuous  time  1-D  lossless  circuits  of  classical  network  theory  [8]  is 
developed.  This  representation  is  then  subsequently  used  in  section  4.4  to 
synthesize  the  filter  in  a  specific  structure. 

A  two-port  system  of  the  type  under  consideration  is  lossless  if  (4.3.1) 
holds  true  for  any  finite  (pseudo)  energy  inputs  x.(n.,n2)  and  x2(nlfn2). 

2  2  2  2 
££ | y^ ( n^ , n2 ) |  +  £E|y2(n^,n2) |  *  E£|x^(n^,n2) |  +  EE|x2(n^,n2) | 


(4.3.1) 

Considering  x2(n^,n2)  -  0  for  all  n^,  n2  we  then  have  from  (4.3.1)  that  for 
any  finite  (pseudo)  energy  x^(n^,n2): 

EE|y1(nl,n2)|2  <  EElx^n^)  |2,  EE|y2(nrn2)  1 2  <  EElx^n^)  |2 


(4.3.2a,b) 


On  the  other  hand,  if  S  -  S(z^,z2)  -  tS^j(z^,z2)]  is  the  transfer  function  of 
the  two-port  then  for  XjUjyZ^  ■  0  we  have  Y^z^Zj)  -  S11(z1,z2)X1(z^,z2  and 
Y2(Zi,z2)  -  s2i( zi rZ2)X1(z1,z2) ,  where  as  in  last  section  capital  letters  are 
used  to  denote  2-D  z-transforms.  Thus,  due  to  (4.3.2a)  and  (4.3.2b)  the 
transfer  functions  -  S^(z^,z2)  and  ■  S21*zl'z2*  are  (pseudo)  passive, 
and  thus  satisfy  property  4.2.1.  Similary,  by  considering  x^(n^,n^)  -  0  for 


all  n 


1' 


n2  it  can  be  shown  that  S^2 


S12^zl'z2^  arK*  S22  *  S22^zl'z2^  sat^s^Y 


Property  4.2.1. 


Furthermore,  by  considering  2-D  Parseval's  theorem  (4.3.1)  can  be  made  to  yield 


98 


(4.3.3)  where  the  column  vector  ■  (X^u^,^)  X^o^,^))*",  and  * 

denotes  the  combined  operation  of  complex  conjugation  and  matrix  transposition. 

®  ^  *  * 

J  /  X  (wl/w2)(I2  -  S  -  0  (4.3.3) 

— n —  rt 


Since  (4.3.3)  holds  for  any  Xfw^.o^)  it  follows  that  for  any  2-tuple  (o^,^) 
except  possibly  finitely  many,  we  have  S  (w^w^JStMj^M,)  -  l2-  This  latter 
equation,  by  exploiting  arguments  similar  to  those  used  in  the  paragraph 
following  (4.2.14)  yields  that  for  all  z^,  z2: 


S(z1,z2)  S(zlfz2)  -  I2  (4.3.4) 

A  (2  x  2)  rational  matrix  S  ■  S(z^,z2)  is  said  to  be  fully  recursive  half-plane 
lossless  bounded  if:  (i)  each  entry  of  S  satisfies  the  conditions  expressed  in 
properties  4.2.1  (i),  (ii)  and  equation  (4.3.4)  holds  true. 


Note  that  the  transfer  function  of  a  fully  recursive  half-  plane  lossless 
two-port  is  necessarily  of  the  above  type.  As  a  consequence  of  property 
4.2.2,  we  then  have  the  following  important  conclusion. 


Proposition  4.3.1:  Each  entry  of  a  fully  recursive  half-plane  bounded  matrix  S 
■  (S^jJ  satisfies  |S^j|  <  1  for  all  |z^|-l  and  |z2|<l. 

Proof:  Since  on  |z^|  -  |z2l-l,  we  have  S*S  -  I2,  and  thus  |Sj^|^+  |S2^|^  *  1» 
,S22,2+  |S12|2  "  l'  which  in  tum  imply  (S^l  <  1  for  all  i,j.  The  result  then 
follows  from  property  4.2.2. 


Consider  next  a  fully  recursive  half-plane  lossless  bounded  matrix  S.  Since 
each  entry  of  S  satisfies  property  4.2.1,  the  rational  fnction  det  H  also 
satisfies  property  4.2.1.  Also,  it  follows  from  (4.3.4)  that  (det  S)(det)  »  1. 

Thus,  S  is  a  fully  recursive  symmetic  half-plane  all-pass  function  as  defined 
in  section  4.2  and  admits  of  the  representation  (4.2.17)  described  in  property 
4.2.3,  i.e.,  (4.3.5)  holds. 


det  S  -  -D.  (A/A) 


(4.3.5) 
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Consequently/  it  follows  from  (4.3.4)  and  (4.3.5)  that 

Sn  -  -  (1/D)(A/A)S22;  S21  -  (1/D)(A/A)S12  (4.3.6a,b) 

S12  -  (1/D)(A/A)S21;  s22  -  -  (1/D)(VA)S11  (4.3.7a,b) 

It  may  be  shown  via  the  operation  ~  that  (4.3.6a)  is  identical  with  (4.3.7b), 
whe re  as  (4.3.6b)  is  identical  with  (4.3.7a).  Now  since  in  (4.3.7b) 
satisfies  property  4.2.1,  it  follows  that  when  S11  is  expressed  in  irreducible 
rational  form,  each  of  the  irreducible  factors  of  the  denominator  polynomial 
must  have  zeros  in  (z^l-l,  |z2|>l.  However,  due  to  property  4.2.1,  S22  is 

analytic  in  |z.|-l,  | z2 I <1  and  thus  S22  is  analytic  in  ( | —1 ,  |z2|>l.  This  is 

possible,  however,  only  if  the  denominator  polynomial  of  in  irreducible 

rational  form  is  completely  cancelled  by  the  numerator  of  (A/A) .  Thus, 
(4.3.8a)  follows,  where  B  is  a  pseudopolynomial.  By  inserting  (4.3.8a)  in 
(4.3.7b)  and  performing  the  operation  “,  (4.3.8b)  follows  via  the  use  of  the 
identity  DD  ■  1.  Equations  (4.3.9a,b)  follow  from  (4.3.7a)  in  a  similar 
manner,  where  again  C  is  a  pseudo-polynomial. 

Sn  -  B/A,  S22  -  -  D(B/A)  (4.3.8a,b) 

Sn  -  C/A,  S12  -  D(C/A)  (4.3.9a,b) 

Inserting  (4.3.8)  and  (4.3.9)  in  the  expression  for  (det  S)  in  (4.3.5)  we  then 
have 

M  -  BB  +  CC  (4.3.10) 

Also,  since  S--  and  S15  are  analytic  in  |z. |*1,  |z,|<l  we  have  from  (4.3.8b) 
and  (4.3.9b)  that: 

deg2B  <  deg^  ;  deg2C  <  deg^  (4.3.11) 

The  above  discussion  can  be  succinctly  expressed  in  the  following 
representation  of  a  fully  recursive  symmetric  half-  plane  lossless  bounded 
matrix. 

Property  4.3.1:  Any  fully  recursive  syranetric  half-plane  lossless  bounded 
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matrix  can  be  represented  in  terms  of  three  pseudopolynomials  A,  B  and  C  as  in 
(4.3.8)  and  (4.3.9),  where  A  is  nonzero  in  (z^J-l,  |z2l<l;  the 
pseudopolynomials  A  and  A  are  coprime;  and  furthermore  (4.3.10),  (4.3.11)  hold 
true.  Conversely,  any  matrix,  which  admits  of  the  above  representation  is  fully 
recursive  symmetric  half -plane  bounded. 


Since  property  4.2.1  is  automatically  satisfied  by  each  entry  of  a  matrix  S 
expressed  via  (4.3.8),  (4.3.9)  and  (4.3.10),  to  establish  the  above  converse 
proposition  we  only  need  to  observe  that  straightforward  algebraic  manipulation 
with  (4.3.8)  through  (4.3.10)  yield  (4.3.4).  For  convenience  of  exposition  any 
S  expressed  as  in  (4.3.8)  and  (4.3.9)  will  be  referred  to  as  in  standard  form. 

A  fully  recursive  half-plane  lossless  two-port  can  be  alternatively  described 
by  means  of  a  chain  matrix  T  -  T(z^,z2)  defined  as  in  (4.3.12).,  where  X2 
and  Y^,  Y 2  ace  respective  inputs  and  output  signals  from  ports  1  and  2. 


(4.3.12) 


It  can  be  easily  shown  from  (4.3.8)  through  (4.3.12)  that  the  following 
property  characterizes  the  chain  matrices  of  the  type  described  above. 

Property  4.3.1' :  The  chain  matrix  T  -  { j 1  associated  with  a  fully  recursive 
half-plane  two-port  is  lossless  if  and  only  if  it  can  be  expressed  as 


Tn  -  DA/C  ;  T12  -  B/C 


(4. 3. 13a, b) 


T,.  -  DBz/'A/C  ;  T,,  -  A/C 


(4. 3. 14a, b) 


where  nA  -  deg^,  and  A,  B,  C  and  D  satisfies  the  same  restrictions  described 
in  property  4.3.1.  Also,  any  T  as  in  (4.3.13),  (4.3.14)  is  said  to  be  in 
standard  form. 


4.4. 


SYNTHESIS  OF  FULLY  RECURSIVE  HALF-PUNE  LOSSLESS 


TOD-PORTS: 

A  procedure  for  synthesizing  fully  recursive  half-plane  lossless  two-ports  as 
an  interconnection  of  more  elementary  building  blocks  of  the  same  type  will  be 
developed  in  this  section.  The  synthesis  algorithm  can  be  viewed  as  a 
generalization  of  the  algorithm  for  synthesizing  1-D  discrete  lossless 
two-ports  as  described  by  Rao  and  Kailath  in  (6).  Our  synthesis  procedure 
exploits  the  unique  feature  of  the  algorithm  described  in  [6]  that  (in  1-D) 
given  (polynomials)  A,  B,  C  the  arithmetic  operations  needed  to  be  performed  on 
the  coefficients  of  A,  B  and  C  in  each  cycle  of  the  algorithm  requires  rational 
operations  only.  To  the  best  of  our  knowledge  this  is  the  only  algorithm  of 
the  above  mentioned  type  available  for  synthesis  of  1-D  discrete  as  well  as 
continuous  domain  lossless  two-ports  including  those  in  classical  network 
theory  (8].  (All  other  algorithms  known  prior  to  (6]  required  nonrational 
operations  e.g.,  polynomial  factorization).  The  basic  structure  of  the  filter 
to  be  presently  synthesized  would  thus  be  the  same  as  in  (6],  whereas  the 
elementary  building  blocks  are  certain  1-D  two  port  sections  to  be  referred  to 
as  the  generalized  Gray-Plarkel  sections  (GGM  section)  and  z2~type  delays,  each 
of  which  are  fully  recursive  half-plane  lossless. 

A  generalized  Gray-Markel  section  is  a  1-D  two  port  as  shown  in  figure  4.4.1 
where  the  1-D  transfer  functions  (assumed  rational)  k1  and  k2  satisfy  the 
relations?,  *•: 


klkl  +  k2k2  “  1  (4.4.1) 

~  ~  2 

and  are  such  that  k^  (and  thus  k2  in  view  of  (4.4.1))  satisfies  kjk^  -  ( |  < 

1  everywhere  on  |z.  |  -  1  with  the  possible  exception  of  finite  number  of  values 

i  -  "2 

of  z ^  for  which  we  may  have  kjk^  -  jk^|  -  1. 

We  first  note  that  given  any  rational  function  k^  of  z^  satisfying  the  above 
conditions  it  is  always  possible  to  find  a  rational  function  k2  satisfying  the 
same  conditions  as  that  of  k^  along  with  (4.4.1).  (The  role  of  kj^  and  k2  can 
obviously  be  interchanged  in  the  present  considerations).  To  show  this  let  k, 
-  n^/d^  where  n^  are  polynomials  in  Zj.  Then  (1  -  k^)  -  Nj/td^),  where 
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-  djd^-n^n^.  Thus,  N^  »  N^f  and  for  all  2^  on  |z^|  -  1,  N^lz^)  is  real  and 
we  have  that  N^z^)  >  0  as  a  consequence  of  <  1.  Therefore,  the 
(spectral)  factorization  N^  -  n2n2'  n2  a  P01!™0®*-®!  in  zi  holds. 
Also,  by  (possibly)  rearranging  the  irreducible  factors  of  (d^d^)  to  write  d^d^ 
-  d2d2'  where  d2  “  polynomial,  we  can  have  k2  -  n2/d2  such  that  (4.4.1)  is 
satisfied.  Note  that  since  the  factorizations  N2  -  n2n2  and  djd^  -  d2d2  are 
not  unique  the  k2  so  obtained  is  not  unique  unless  further  restrictions  are 
imposed. 

The  transfer  function  matrix  SG  -  S^Zj)  associated  with  such  a  GGM  section  can 
be  expressed  as  in  (4.4.2a),  whereas  the  corresponding  chain  matrix  T  is  given 
in  (4.4.2b). 


Lk2  -kiJ 


TG“  <1/k2> 


1 


( 4 . 4 . 2a , b ) 


Since  SG  in  (4.4.2a)  satisfies  the  representation  described  in  property  4.3.1 
with  A  -  1,  B  »  k^,  C-k2  and  D  »  1  the  GGM  section  is  indeed  a  fully 

recursive  half-plane  lossless  two-port. 

To  proceed  with  the  synthesis  of  a  prescribed  fully  recursive  half-plane 
lossless  bounded  matrix  S  or,  equivalently,  corresponding  chain  matrix  T  as 
described  respectively  in  property  4.3.1  or  4.3.1',  we  first  note  that  in  view 
of  porposition  3.1  the  rational  function:  k1  -  k^tZj)  -  S11(z1,0)  satisfies 
|kjJ  <  1  for  all  |z^|  -  1  with  the  possible  exception  of  finite  number  of 

values  of  z ^  where  |k^|-l.  Therefore,  in  view  of  the  preceeding  discussion 
defines  a  GGM  section  i.e.,  a  k2  can  be  found  such  that  |k2|<l  everywhere  on 
|z^|*l  with  the  possible  exception  of  finite  number  of  values  of  z^,  for  which 
| k2 | — 1  and  that  (4.4.1)  is  satisfied. 

Step  1:  The  first  step  is  to  extract  a  GGM  section  with  k^  -  sn(zi'0)  froin 

prescribed  S  or  T  as  shown  in  figure  4.4.2.  Since  a  cascade  connection  of  two 
two-ports  amounts  to  multiplication  of  the  corresponding  chain  matrices,  the 
chain  matrix  of  the  remaining  two-port  is  then  T'  -  T „-1T.  From  (4.3.13), 
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(4.3.14)  and  (4.4.2b)  we  can  write: 


T'  -  (1/Ck2) 


-  "A 
D(A-k1Bz2  *) 

~  nA  “ 

_D(Bz2  -k^A) 


B-kjA 

A-k1B 


(4.4.3) 


We  next  define  the  pseudopolynomials  A',  B',  C'  and  the  1-D  rational  function 
D'  as  in  (4.4.4)  and  (4.4.5)  below,  where  p-p(z^)  is  the  conjugate  reciprocal1 
polynomial  factor  of  largest  degree  present  in  the  numerator  of  A-k^B,  when 
expressed  in  irreducible  rational  form. 

-  n 

pA'  -  A(l-k1S11)  -  A-k1Bz2  ;  pC'  -  Ck2  (4.4.4a,b) 

pB'  -  A(S11-k1)  -  B-kjA  ;  D'  -  D(p/p)  (4.4.5a,b) 


We  claim  that  degjA'  -  degjA.  To  prove  this,  clearly  deg^'  <  deg^  and  note 
that  (4.4.4a)  yields  pA'/A  -  1-kjSj^,  which  implies  that  if  deg^'  <  deg^  then 
for  arbitrary  Zj  we  would  have  k1(z1)S^1(z^,0)  -  | , 0 )  | 2  -  1.  The 
impossibility  of  this  latter  situation  has  already  been  demonstrated.  As  a 
consequence  of  this  we  can  write  T'  as  in  (4.4.6)  and  (4.4.7),  where  nft,  - 
deg2A'  -  deg2A. 


T'x  -  D'A'/C'  ;  T'2  -  B'/C'  (4.4.6a,b) 

T21  -  D'B'Zj  A'/C'  ;  IJj  .  A’/C-  (4.4.7a,b) 

We  next  claim  that  the  pseudopolynomial  Ap  satisfies  the  properties  that  A' 
*  0  for  |ZjJ-l,  | Z2 1 <1  and  that  A'  is  coprime  with  A'.  To  prove  this  we  write 

A'  -  A^/A£  ln  irre^uc^^e  “ational  form,  and  thus  A'  -  (A^/A^).za,  where 
a  -  integer  and  A^/A£  is  in  irreducible  rational  form.  Thus,  from  the 

1 

A  polynomial  p  is  said  to  be  conjugate  reciprocal  if  p*r-p  for  some  constant 


definition  of  A'  it  follows  that  A^  is  devoid  of  conjugate  reciprocal 
polynomial  factors  in  z^  only.  If  we  assume  for  the  prupose  of  a  proof  by 
contradiction  that  for  some  value  of  zi“zio'z2*z20  I z10 1*1'  lz20*<1  we 
have  A'-O  i.e.,  A^-0  then  since  A^  cannot  have  a  factor  (z^-z^),  by  changing 
the  value  of  z^  from  z^q  along  an  arbitrarily  small  arc  of  the  unit  circle 
|z^|«l  it  would  be  possible  to  find  a  continuous  [12]  set  (z^^)  of  zeros  of 
A^  i.e.,  also  of  A'  with  z^eTf[z^;  |z^  |-1)  and  IZ2I<1. 

Also,  since  it  follows  from  (4.4.4a)  and  deg2B<nA«nA'  that  pA'-Ad-k^S^) 
and  A*0  in  |z^|»l,  | I <1  (cf.  property  4.3.1)  we  would  then  have  that  for  all 
z^ef^  there  exists  some  Z2  in  |z2l<l.  such  that  kjS^-1.  Since  | |  <1 ,  | (  <1 
for  all  |z.|-l, |z2l<l  (cf.  Proposition  4.3.1)  the  last  conclusion  implies  that 
there  exists  z2  in  )z2l<l  511011  tliat  lki  l“|S11(z1,0)  |-1  and  |S11(z1,z2)  !-l  for 
all  z^el^.  However,  this  last  conclusion  leads  to  a  contradiction  in  view  of 
Property  4.2.2a.  Thus,  A'*0  i.e.,  A^*0  for  |z. |-1, | Zj I <1 • 

Also,  since  due  to  Property  4.2.2a  IkjJ-ISj^Uj^O)  |-1  is  possible  for  at 
most  finite  number  of  values  of  Zy  for  |z^|-|z2l-l  we  have  IS^IO.  and  A*0  for 
at  most  finite  number  of  exceptions  on  the  distinguished  boundary  of  the  unit 
bi-disc,  we  conclude  from  pA'-Ad-kjS^)  that  A',  thus  A^  may  have  at  most 
finite  number  of  zeros  on  jz^l-^l-l.  Since  as  shown  earlier  A£*0  in 
|z1l-l,  |z2l<l  it  follows  from  Assertion  4.2.1  that  Aj^  and  A^j  are  relatively 
prime  polynomials.  Consequently,  the  pseudopolynomials  A  and  A'  are  relatively 
prime. 


Finally,  straightforward  algebraic  manipulation  along  with  (4.4.1)  yields 
AA'  -  BB'  +  CC'  ,  whereas  degjA'  <  deg2B',  deg2A'  <  deg2C'  follow  from 
(4.4.4a,b),  (4.4.5),  (4.3.11)  and  nA-nA,.  Since,  clearly  D'  as  in  (4.4.5b) 
possesses  the  requisite  properties  for  T'  to  be  in  standard  form,  in  view  of 
Property  4.3.1'  all  the  conditions  necessary  for  T  -  [T^],  as  given  in 
(4.4.6),  (4.4.7),  to  be  a  fully  recursive  half-plane  lossless  two-port  chain 
matrix  are  satisfied. 

We  further  note  that  as  a  consequence  of  the  choice  k^-S^^(z^,0)  we  have 
from  (4.4.5a)  that  B'(z^,0)  -  0  for  arbitrary  z^  i.e.,  the  pseudopolynomial  B' 
contains  Z2  as  a  factor.  Also,  if  C  contains  a  pseudopolynomial  factor  Z2  then 
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so  does  C' . 


(2) 

Step  2:  In  the  next  step  we  form  a  fully  recursive  naif-plane  two  port  T  by 
interchanging  the  two  output  terminals  in  each  port  of  T'  as  shown  in  figure 


4.4.3.  It  can  be  easily  shown  that  T 

,(2)  „(2)  „(2) 


(2) 


can  then  be  written  in  terms  of 
(2) 


pseudopolynomials  A  ,  Bv  ,  C  and  the  rational  function  D 
form  as  expressed  in  property  4.3.1',  where 


in  standard 


,  (2) 


B 


<2>  -C', 


.(2) 


,(2) 


-D' 


(4.4.8) 


Step  3:  A  GGM  section  is  then  extracted  from  the  two-port  with  chain  matrix 


by  iterating  step  1  on  to  get  a  fully  recursive  half-plane  lossless 

C*3*  and  rational 


two-port  chain  matrix  3  ^ .  If  pseudopolynomials  A*3^,  B^ 3  ^ 

function  D*3*  represent  3 ^  in  standard  form  as  in  property  4.3.1'  then  B^3* 

would  have  a  factor  z2<  Also,  since  from  (4.4.4b)  and  (4.4.8)  we  have  C^3*  - 

C*2*lu*3*  -  B'k,(3\  where  k./3^  defines  the  GGM  section  presently  extracted, 

/  z  z  ,  2 . 

and  B'  has  a  “actor  z-  we  conclude  that  C  '  has  a  factor  z,.  From  this  and  the 


CUM  a  itao  a  auuui.  wo  uictt  uoo  a  Lauuvi  •  r  suiu  u 

fact  that  a(3)a(3)  -  b(3)B(3)  +  C(3)C(3)  it  follows  that  a'3,A(3)  • 


0  and  for  arbitrary  z. .  Since  A  3  *  0  for  |z.  |  -  1  and  |z9|  <  1,  we  conclude 

(3)  1  r3) 

A'  ■  0  for  z,  -  0  and  for  arbitrary  z. .  Consequently,  A'  has  a  factor  z 

•  A  /  A  \  i  A  \  t  A 


0  for  z2  * 


and  it  is  possible  to  write,  for  some  pseudopolynomials  A 
that 


<4>  B<4>andC<4? 


a(4)-Z2a(3), 


B(4)  -z2B(3),  C<4) 


Z2C 


(3) 


(4.4.9) 


Step  4 


The 


?3) 


from 

which  in  standard 
Furthermore, 


last  step  in  the  synthesis  cycle  is  to  extract  a  z.,  type  delay 

.  c -  A  t  A  L. . 1- _  . _ ..  •  .L  .k.J. .-J  m(  4  ) 


as  in  figure  4.4.4  to  produce  a  two-port  with  chain  matrix  T 


form 
a<4>-  a<3> 


can  be 
t  0  for 


expressed  in  terms  of 


(4)  _(4)  „(4) 


most  finitely  many  zeros  on  |z.|  -  |z,  | 

i  \  A  Z 


m  1, 

-  1. 


lz-,1  <  i; 


a<4> 


A'  B 
.  r(3) 


can  have  at 


losslessness  of  T 


(3) 


Also,  it  follows  from  (4.4.9)  and 


that 


and 


;<4>a<4>  .£<4>b<4>  +c(4)c(4> 


deg2B^4^  <  deg^4^,  deg2C^  <  degjA^4^ 


(4.4.10) 
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Thus,  the  two-port  associated  with  T 
Furthermore,  note  that 


is  fully  recursive  half-plane  lossless. 


(4) 

deg^^  ”  degjA^-l  ■  deg2A^-l  «  deg2A'-l  -  deg2A-l, 

where  the  first  equality  follows  from  (4.4.9);  the  second  and  the  fourth  from 
the  fact  that  in  step  1  we  have  nA  ■  nA, ;  and  the  third  from  (4.4.8). 

Consequently,  after  iterating  deg^  times  the  cyclic  algorithm  described  in 
Steps  1  through  4,  we  obtain  a  lossless  chain  matrix  independent  of  z2, 
which  in  standard  form  is  described  by  A^  «  Ag(z1),  «  Bf(z1),  ■  C^fz^) 

and  D  -  D^(z^) . 

Terminal  Step:  In  the  final  step  we  extract  another  GGM  section  as  in  Step  1 
to  produce  a  fully  recursive  half-plane  lossless  two-port  SQ  with  AQ,  Bq,  Cq 
and  Dq  in  standard  form.  Since  Af,  B^  are  functions  of  z ^  only  it  follows  from 
(4.4.5a)  that  BQ  ■  0.  Also,  since  A^Aq  »  BqBq  +  CgCg  and  Ag  -  Ag  *  0  for  all 
|z^|  ■  1  the  1-D  transfer  functions  (Sg)12  ■  DoC(/A0  31x1  ^S0^21  "  Ct/A0  are 

both  well  defined  and  of  unit  modulus  on  |z.  |  -  1  i.e.,  they  are  all-pass 

2  1 
functions.  The  realization  for  such  a  two-port  is  shown  in  figure  4.4.5.  The 

resulting  composite  filter  structure  is  as  shown  in  figure  4.4.6. 


^Note  that  (sq^12  an<*  (Sg)21  are  not  necessar^1y  stable  rational  functions, 
i.e.,  may  have  poles  in  |z^|<l. 


4.5. 


DESICT  OF  2-D  FULLY  RECURSIVE  HALF-PLANE  FILTERS: 


Two-dimensional  filters  with  only  two  kinds  of  symmetries  in  their 
magni tude  responses,  namely  the  fan  type  synraetry  and  the  circular  symmetry  are 
of  practical  interest.  The  locii  of  canstant  gain  in  the  plane  for  the 
fan  filters  are  required  to  be  approximate  straight  lines,  whereas  those  for 
the  circularly  symmetric  filters  are  required  to  be  closed  circles  in  an 
approximate  sense.  In  addition,  we  also  require  the  pass  (or  the  stop)  region 
of  the  fan  filter  to  be  the  region  approximately  lying  within  the  straight 
lines  -  oa*^  and  -  -  au^  for  some  0<a  <1. 


Our  design  proceeds  by  requiring  the  transfer  function  S21  -  C/A  (cf. 
equation  (4.3.9a))  of  the  lossless  two-port  S  to  have  the  desired 
characteristics.  However,  unlike  the  corresponding  problem  in  1-D,  due  to 
nonfactorability  of  m-D  polynomials  it  is  in  general  not  possible  to  find  a 
pseudo-polynomial  B  satisfying  (4.3.10)  from  A  and  C.  To  circumvent  this 
problem  it  will  be  further  assumed  that  the  two-port  is  either  symmetric  i.e., 


,,  S21  -  S22  or  antimetric  i.e.,  Sn--S2 


symmetricnand  in  the  antimetric  case  we  respectively  have  B  -  -DBz2 


S21  "  S12*  nT,lus'  in 
_ »  _  .JL  A 


B  -  DBz, 


whereas  we  also  have  C 


in  both  cases.  We  next  define  two 


o  ■  t  mreivoo  "v  jwvc  v.  —  ia.^2  im  iajui  tasca •  nc  iical  uciiuc  lwu 

rational  functions  and  S2  as  in  (4.5.1)  and  5.2)  respectively  for  symmetric 
or  antimetric  two  ports. 


Sx  -  (B  +  C)/A,  S2  -  (B  -  C)/A  (4.5.1a,b) 

S2  -  (B  +  jC)/A,  S2  -  (B-  jC)/A  (4.5.1'a,b) 

From  (4.5.1)  it  is  easily  verified  that  SjS^  -  S2S2  -  1.  Thus  for  each  i,  | SA | 
■  1  for  all  | 2^ |  -  |z2|  -  1  except  possibly  finitely  many  values  where  it  is 
undefined.  Furthermore,  A  *  0  for  | |  »  1,  |z2|  <  1.  Thus  via  Property  4.2.2a 
it  follows  that  S^,  for  each  i  -  1,2  in  (4.5.1)  must  be  a  fully  recursive 
half-plane  all-pass  function.  Exactly  same  conclusions  hold  for  and  S2  in 
(4.5.1').  Consequently,  S^,  S2  can  be  expressed  as  in  (4.5.3),  where  D^,  D2 
and  A^  Aj  satisfy  properties  analogous  to  D  and  A  in  Property  4.2.3. 

S1  -  -  S2  -  -  d2A2//^2  (4.5.3a,b) 
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Note  that  even  if  A,  B,  C  are  real  rational  functions,  and  S2  are  real  in 
(4. 5. lab)  but  not  in  (4.5.1'a,b).  Thus,  a  syrnnetric  filter  can  be  realized  by 
making  use  of  the  relation  S21  -  C/A  -  (Sj-S2)/2,  where  the  one-ports  and  S2 
are  realized  as  in  Appendix  A.  Although  S21  -  C/A  -  -  j ( S^— S2 )  holds  true  in 
the  antimetric  case,  a  realization  in  terms  of  this  last  mentioned  equation  is 
not  feasible  due  to  the  presence  of  the  factor  j  unless  complex  filter 
realizations  are  called  for.  In  this  case,  the  pseudopolynomials  A,  B,  C  which 
are  real,  can  be  found  from  (4.5.1'a,b)  and  subsequently  S21  can  be  realized  as 
being  embedded  in  a  real  two-port  S  described  by  A,  B,  C  in  standard  form.  The 
design  problem  then  boils  down  to  appropriately  choosing  the  real  1-D  rational 
functions  D^,  D2,  and  real  pseudopolynomials  a^,  A2  so  that  the  frequency 

response  requirements  on  |S21(  are  satisfied.  This  latter  step  may  be  carried 
out  by  using  numerical  optimization  (e.g.  Levenberg-Marquadt) .  For  the  purpose 
of  numerical  optimization,  however,  the  following  symmetry  observations  have 
the  effect  of  reducing  the  number  of  parameters  to  be  optimized. 

Note  that  if  |  possesses  either  the  fan-type  or  the  circular-type 

symmetry  then  I S21  ( )  I  roust  have,  in  particular,  the  so  called  quadrantal 
symmetry  ( 14 J  i.e.,  Is2i^wi,w2^ I  is  unaltered  if  the  signs  of  either  or 
or  both  are  changed.  This  requirement  coupled  with  the  stability  property  of 

and  A2  in 

(4.5.3)  satisfy  a  certain  factorability  property.  It  then  proves  to  be 
expedient  to  carry  out  the  optimization  procedure  on  these  factors  rather  than 
on  the  pseudopolynomials  A^,  A^. 


S21  demands  that  the  denominater  pseudopolynomial  of  S2^,  and  thus  A^ 
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AFPDDIX: 


In  this  appendix  we  prove  that  a  fully  recursive  symmetric  half-plane  all-pass 
function  H-H( z^,z2)  can  ^  synthesized  as  an  interconnection  of  GGM  sections 
(cf.  Section  4.3)  and  z2~type  delays.  This  can  be  considered  to  be  a 
generalized  form  of  Schur's  algorithm  [91. 

Let  k ^  -  k^z^  «  Hfz^O).  Since  H  is  as  in  Property  4.2.3  it  follows  from 
Property  4.2.2a  that  |k^|<l  for  all  |z^j-l  with  the  possible  exception  of 
finite  number  of  values  of  z^,  where  |k^|-l.  Thus,  a  k2  satisfying  (4.4.2)  can 
be  found  i.e.,  k^  and  k2  defines  a  GGM  section.  Consider  next  the  function  Hj 
-  H-fz^Zj)  defined  as  in  (4.A1.1),  which  can  be  interpreted  as  the  residual 
transfer  function  after  extraction  of  the  GGM  section  just  mentioned  from  H^. 

H'  -  (H  -  k^/U  -  kH)  (4.A1.1) 

From  (4.2.17)  it  then  follows  that  H'  «  -  pA'/B' ,  where  pA'  -  DA  +  kjA,  B'  «  A 
+  k.  DA,  p  being  the  conjugate  reciprocal  factor  of  largest  degree  present  in 
the  numerator  of  (DA  +  k.A)  when  expressed  in  irreducible  rational  form.  Next, 
since  we  have  D  -  D-i  it  follows  that  pA'/A  -  D(l-kjH).  Consequently,  if  deg2 
A'  <  deg  A  then  we  would  have  k^HU^O)  —  |H(zlf0)  |2—  1  for  arbitrary  Zy  which 
is  impossible  (cf.  Property  4.2.2a).  Thus,  deg^'-  deg  A.  It  then  clearly 
follows  that  pA'  -  D(A  +  k^DA)  -  DB' ,  thus  H'  -  -  D^A'/A');  D1  -  D(p/p). 
Also,  since  A'  ■  -  A(H-k^)  it  follows  that  A'  -  0  for  z2-0  and  arbitrary  z^. 
Thus,  the  pseudopolynoraial  A'  has  a  factor  z2*  By  defining  Aj  via  zA^  *  A'  we 
can  write  -  z2H' ,  where  -  -  D^(A^/A^).  Note  that  can  be  constructed 
simply  by  extracting  a  z2~type  delay  from  H' .  Clearly,  satisfies  condition 
(ii)  of  property  4.2.3.  Also,  by  following  arguments  similar  to  that  used 
after  (4.4.7a,b)  it  can  be  shown  that  A^  ■  A'  t  0  for  [z^  -  1,  |z2|  <  1  and 
the  pseudopolynomials  A^  and  A^  are  coprime.  Thus,  condition  (iii)  and  (iv)  of 
Property  4.2.3  are  satisfied  by  H^,  which  has  now  been  proved  to  be  fully 
recursive  symmetric  half-plane  lossless.  Since  deg^  -  deg2A'  -  1  -  deg^  -  1 
the  procedure  just  described  when  applied  degjA  times  yields  a  circuit  as  shown 
in  figure  4.A.1,  in  which  the  terminating  section  is  an  all-pass  (not 
necessarily  stable)  in  z^  only. 
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CHAPTER  5^ 

CONCLUSIONS  MD  FURTHER  WORK 


5.1  CONCLUSIONS: 

Due  to  potential  applications  of  multidimensional  passive  filtering  schemes  in 
manifold  areas  of  signal  processing  e.g.,  including  frequency  filtering, 
modelling  of  random  fields  associated  with  detection  and  estimation  of 
parameters  of  multidimensional  signals,  fundamental  issues  relating  to  the 
description  of  passive  m-D  systems  and  their  synthes izability  have  been 
addressed  in  the  present  report.  First,  various  classes  of  stable 
multidimensional  polynomials  essential  to  the  discrete  domain  description  of 
passive  m-D  systems  has  been  delineated.  The  synthes izability  of  quarter  plane 
type  causal  m-D  digital  filters  has  been  investigated  in  a  very  general 
setting.  Finally,  due  to  its  potential  benefit  to  be  derived  from  currently 
emerging  parallel  hardware  architectures  [1],  [2],  a  method  of  passive 
filtering  within  the  framework  of  an  alternative  recursive  scheme,  namely,  the 
fully  recursive  half-plane  recursive  scheme  has  been  introduced  and  studied. 
Filter  synthesis  procedures  within  this  framework  has  also  been  examined. 

A  major  problem  in  processing  such  m-D  discrete  signals  in  real  time  is  the 
large  amount  of  data  rate  involved.  Conventional  digital  filtering  algorithms, 
which  process  data  sequentially,  is  therefore,  inappropriate  for  such 
processing  purposes.  In  spite  of  the  fact  that,  as  in  the  1-D  case,  the  basic 
linear  algebra  operations  such  as  vector  matrix  multiplication  form  the  core  of 
many  m-D  signal  processing  algorithms  available,  a  detail  study  of  their  fast 
VLSI/optical  implementation ,  which  utilizes  the  underlying  structure  of 
multidimensional  problems  is  yet  unavailable.  On  the  other  hand,  a  consistent 
scattering  formalism  for  passive  multidimensional  systems  has  begun  to  emerge 
as  a  result  of  the  presently  reported  work,  it  should  now  be  possible  to 
undertake  an  investigation  into  the  design  of  a  broad  variety  of  concurrent, 
numerically  stable  and  fault  tolerant  m-D  signal  processing  algorithms. 
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5.2.  FUK1HER  WORK: 


Study  of  novel  multidimensional  signal  processing  algorithms  of  the  concurrent 
type  for  implementation  in  VLSI/optical  architectures  via  the  framework  of 
passive  scattering  theory  should  thus  form  the  major  emphasis  of  future 
research  in  this  area.  Two  generic  considerations  ensuing  from  the  present  work 
work  may  potentially  form  the  basis  of  this  investigation.  To  elaborate  on  the 
first,  note  that  traditionally  most  m-D  processing  (filtering)  have  relied  on 
the  quarter-plane  type  recursion  schemes  in  the  causal  order  of  data  points. 
However,  since  unlike  1-D,  in  most  multidimensional  applications  the 
independent  coordinates  describing  the  signal  may  not  have  temporal,  but  only 
spatial  significance,  to  impose  such  a  causality  restriction  is  not  only 
unnecessary,  but  it  may  cause  severe  drawback  in  fully  exploiting  the 
concurrency  offered  by  the  optical  architechtures.  Due  to  this,  the  fully 
recursive  symmetric  half-plane  scheme,  seems  to  be  most  appropriate  for  a  large 
number  of  m-D  signal  processing  problems  of  diverse  nature.  To  justify  this 
remark  in  2-D,  it  may  noted  that  the  computational  model  under  consideration 
have  the  remarkable  property  that  all  data  points  belonging  to  a  row  in  the  2-D 
lattice  space  may  be  simultaneously  (in  parallel)  computed  front  the  data  points 
belonging  to  an  adjacent  row  in  the  2-D  lattice  space.  Furthermore,  as  has  been 
demonstrated  in  the  case  of  frequency  filtering  in  chapter  4,  the  computation 
just  referred  to,  in  fact,  involves  computing  the  convolution  of  the  adjacent 
rows  with  a  fixed  1-D  sequence  determined  by  the  transfer  function  of  the 
filter.  Since  1-D  convolution,  among  other  linear  algebra  operations,  are  known 
to  easily  yield  to  high  speed  optical  architechtures  using  acousto-optic 
devices  [1],  the  time  complexity  of  the  entire  computation  can  be  shown  to  be 
drastically  less  than  filters  using  other  recursive  schemes.  (In  fact,  the 
speed  of  this  2-D  processing  scheme  proves  to  be  of  comparable  order  to  that  of 
a  1-D  filter,  when  an  implementation  of  this  type  is  called  for).  Furthermore, 
since  for  certain  3-D  applications,  (e.g.,  in  the  processing  of  3-D  time 
varying  imagery)  this  type  of  processing  corresponds  to  processing  all  of  the 
data  points  in  a  given  frame  of  the  image  simultaneously,  while  the  direction 
of  filter  recursion  corresponds  to  the  direction  of  flow  of  time,  the 
processing  scheme  under  consideration  is  not  only  adequate  for  fast  optical 
implementation,  but  is  the  most  natural  choice  unlike  the  somewhat  contrived 
quarter  plane  causal  recursive  scheme,  which  has  been  almost  universally 
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adopted  in  the  m-D  scattering  based  signal  processing  techniques  so  far. 

A  second  possibility  of  optical  implementation  of  passive  filtering  algorithms 
arises  from  the  fact  that  1-D  digital  lattice  filters  has  been  shown  to  be 
easily  implementable  via  analog  optical  devices  such  as  the  single  mode  fiber 
and  directional  couplers  [1].  The  obeservation  that  the  digital  lattice  section 
(also  known  as  the  Gray-Markel  section)  and  variations  of  it  can  be  implemented 
in  terms  of  optics,  can  be  potentially  utilized  in  the  implementation  of 
passive  m-D  digital  filters.  This  conclusion  derives  from  the  fact  that,  as 
shown  in  chapters  3  and  4  any  passive  »-D  digital  filter  can,  in  fact,  be 
implemented  as  interconnection  of  modular  building  blocks  each  of  which  can  in 
turn  be  viewed  as  interconnections  of  a  small  number  of  digital  lattice 
sections.  To  carry  this  point  little  further  let  us  note  that  single  mode 
optical  fibers  and  directional  couplers  can  be  interpreted  as  passive  digital 
two-port  networks  [1]  and  their  propagation  characteristics  can,  in  fact,  be 
described  in  terms  of  scattering  parameters,  which  are  exactly  the  tools  in 
deriving  a  large  variety  of  1-D  and  m-D  signal  processing  algorithms  of 
interest  to  us  in  the  present  context.  However,  we  are  unaware  of  any  existing 
work  which  develops  this  connection  further  so  as  to  tailor  specific  m-D 
algorithms  to  fit  into  the  optical  architectures. 
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Figure  4.A.l:synthesis  of  a  fully  recursive 
symmetric  all-pass  one-port 
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